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PREFACE 


The  Materials  Research  Laboratory  at  Brown  University  is  engaged  in  interdisciplinary  research  directed  toward 
the  solution  of  problems  in  various  areas  of  Materials  Science. 

In  this  report.  Sections  l-IV  represent  major  thrust  areas;  the  remaining  sections  are  less  extensive,  since  they 
cover  either  new  areas  of  activity,  emerging  thrusts  with  a  smaller  number  of  investigators,  or  derive  their  principal 
funding  from  sources  other  than  the  NSF-MRL  program. 

Although  the  research  summaries  in  this  report  are  contributions  from  individual  investigators,  introductions  to  the 
four  major  thrust  area  summaries  have  been  provided  by  investigators  who  during  the  reporting  period  assumed 
the  functions  of  coordinators  of  those  thrust  areas.  These  are:  Microscopic  and  Macroscopic  Dynamic  Plasticity, 
R.  J.  Clifton;  Fracture  of  Solids,  J.  R.  Rice;  Inorganic  Glasses,  W.  M.  Risen,  Jr.,  and  Chemisorption  and  Related 
Surface  Interactions,  P.  J.  Estrup. 

Since  many  of  the  research  areas  are  interrelated,  cross  references  have  been  made  when  the  work  might  with 
equal  justice  have  been  placed  in  two  or  more  sections. 

The  largest  single  source  of  support  of  this  Laboratory  is  the  Materials  Science  Branch  of  the  National  Science 
Foundation.  Substantial  support  also  came  from  sixteen  other  agencies  and,  of  course,  Brown  University. 

The  Materials  Research  Laboratory  at  Brown  University  is  administered  technically  by  a  Director,  Associate 
Director,  and  a  committee  composed  of  faculty  members  in  representative  areas  of  Materials  Science. 


The  present  appointees  are: 

Director: 

Associate  Director: 
Advisory  Committee: 


G.  S.  Heller 

C.  Elbaum 

R.  H.  Cole 
L.  N  Cooper 
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SECTION  1 


Microscopic  and  Macroscopic  Dynamic  Plasticity 


Introduction 


Research  in  this  area  is  directed  toward  development  of  sufficient  understanding  of  the  mechanisms  of  plastic 
deformation  to  be  able  to  formulate  constitutive  relations  that  model  plastic  deformation  of  single  crystals  and, 
ultimately,  structural  metals.  Both  analytical  and  experimental  investigations  are  being  conducted.  The  size 
scale  on  which  mechanisms  are  examined  varies  from  the  atomistic  level  in  quantum  rate  theory  calculations, 
through  the  microstructural  level  in  dislocation  mobility  studies  in  single  crystals,  to  the  macroscopic  scale  of 
plastic  deformation  of  polycrystals. 


Re-examination  of  the  analysis  of  the  decay  of  the  elastic  precursor  in  plate  impact  experiments  on  single  crys¬ 
tals  has  drawn  attention  to  the  possible  importance,  and  even  dominance,  of  the  effects  of  surface  damage.  A 
critical  experiment  to  examine  this  possibility  is  being  developed.  In  other  work  on  dislocation  mobility  in  single 
crystals,  a  technique  is  being  developed  for  determining  the  attenuation  and  the  change  in  elastic  wave  speed 
of  ultrasonic  waves  during  and  immediately  following  plastic  deformation  caused  by  a  relatively  large  amplitude 
stress  pulse.  6 


The  relationship  between  strain  rate,  stress,  and  temperature  has  been  examined  for  bcc  metals.  Interesting 
strain  rate  history  effects  have  been  obtained  at  low  temperatures  ( — 190°C)  for  a  hot-rolled  steel.  Possible 
models  for  strain-rate  history  effects  and  temperature  history  effects  have  been  considered. 


On  the  atomistic  level,  quantum  rate  theory  has  been  extended  to  the  case  of  an  asymmetric  double-well  poten¬ 
tial.  In  addition,  the  form  of  the  elastic  moduli  of  centro-symmetric  crystals  has  been  derived  in  terms  of  the  elec¬ 
tron  density  and  its  distortion  due  to  the  crystal  deformation.  Conditions  for  the  moduli  to  satisfy  the  Cauchy  rela¬ 
tions  have  been  established „ 


Other  investigations  have  been  directed  toward  understanding  dynamic  plastic  deformation  at  higher  strain 
rates,  under  non-homogeneous  deformations,  and  for  non-proportional  loading  paths.  Pressure-shear  plate 
impact  experiments  on  thin  specimens  of  commercially  pure  aluminum  sandwiched  between  hard  elastic  plates 
have  shown  a  marked  increase  in  flow  stress  between  strain  rates  of  id’s'  and  2  x  ids  '.  This  increase  in 
flow  stress  suggests  a  change  in  the  rate  controlling  mechanism  for  plastic  deformation  at  these  exceptionally 
high  strain  rates.  The  stability  of  homogeneous  plastic  deformation  at  high  strain  rates  where  inertia  effects  are 
significant  has  been  examined.  Explicit  results  have  been  obtained  for  several  limiting  cases.  Finally,  the  effects 
of  non-proportional  loading  paths  have  been  examined  by  an  analysis  of  plastic  waves  of  two-dimensional 
anti-plane  strain  and  by  pressure-shear  impact  experiments  on  a  -titanium  plates.  Comparison  of  calculated 
final  strains  for  a  smooth  yield  surface  model  and  a  corner  theory  model  indicates  little  difference  in  the  calcu¬ 
lated  final  states  of  anti-plane  strain.  Velocity-time  profiles  for  the  a  -titanium  plates  indicate  that  the  effects  of 
strain-rate  sensitivity  as  well  as  non-proportional  loading  paths  will  have  to  be  included  before  the  material  s 
dynamic  plastic  response  can  be  modeled  satisfactorily. 
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SECTION  1 

Microscopic  and  Macroscopic  Dynamic  Plasticity 


Individual  Contributions 


Mechanics  of  Dislocations  In  Single  Crystals 

The  decay  of  the  elastic  precursor  has  been  re-examined  from  the  viewpoint  of  the  etastodynamics  of  disloca¬ 
tions.  Solutions  for  long  straight  dislocations  set  in  motion  by  an  incident  plane  wave  have  been  used  to  derive  the 
rate  of  decay  of  the  elastic  precursor.  The  resulting  equation  for  precursor  decay  is  essentially  the  same  as  the 
commonly  used  one,  which  is  based  on  an  elastic/visco-plastic  model  of  the  crystal.  The  only  difference  is  a  small 
correction  due  to  the  fact  that  for  the  elastic/visco-plastlc  model  dislocations  of  opposite  sign  moving  in  opposite 
directions  have  the  same  effect  on  precursor  decay,  whereas  according  to  the  etastodynamics  analysis,  disloca¬ 
tions  propagating  toward  the  wavefront  have  a  greater  effect  than  those  propagating  away  from  the  wavefront. 
Based  on  this  confirmation  of  the  validity  of  the  precursor  decay  relation,  an  explanation  for  the  larger-than-pre- 
dicted  decay  observed  in  experiments  has  been  sought  in  terms  of  possible  errors  in  either  the  dislocation  veloc¬ 
ity  or  dislocation  density  that  appear  in  the  precursor  decay  relation.  Analytical  and  experimental  evidence  sug¬ 
gests  that  usual  estimates  for  dislocation  velocity  are  reasonable,  but  that  large  errors  in  estimates  for  dislocation 
density  occur  at  the  impact  face.  Surface  damage  during  specimen  preparation,  and  the  impact  of  rough  sur¬ 
faces,  are  responsible  for  dislocation  densities  at  the  surface  that  are  several  orders  of  magnitude  larger  than  in 
the  interior  of  the  crystal.  The  large  dislocation  densities  near  the  impact  face  are  believed  to  be  responsible  for 
much  of  the  larger-than-predicted  reduction  in  precursor  amplitude.  The  decrease  in  precursor  amplitude  with 
increasing  distance  of  propagation  appears  to  be  compatible  with  dislocation  densities  at  the  wavefront  once  the 
nonlinear  elasticity  effects  of  unloading  waves  from  the  impact  plane  are  taken  into  account.  Overall,  it  appears 
that  precursor  decay  can  be  explained  by  considering  finite  elasticity  effects  and  large  dislocation  densities  near 
the  impact  surface.  In  order  to  test  the  validity  of  this  explanation,  an  experiment  has  been  designed  to  reduce  the 
generation  of  dislocations  at  the  front  and  rear  faces  of  the  specimen  by  placing  the  specimen  plate  between  two 
liquid  layers.  Preliminary  experiments  suggest  the  feasibility  of  conducting  such  critical  experiments. 

Principal  Investigator:  R.  J.  Clifton  (Engineering). 

Personnel:  X.  Markenscoff  (Engineering),  P.  Kumar  (Engineering),  and  L.  Hermann  (Engineering). 

Publications:  "The  Transient  Motion  of  a  Nonuniformly  Moving  Dislocation,”  X.  Markenscoff,  J.  of  Elasticity  10 
(1980),  193. 

“On  the  Kinetics  of  a  Frank-Read  Source,”  P.  S.  Stelf  and  R.  J.  Clifton,  Materials  Sci.  and  Engrg.  41  (1979),  251. 


Dislocation  Bahavlor  under  Stress-Pulse  Loading 

We  have  redesigned  an  apparatus  for  studying  the  dynamic  behavior  of  dislocations  during  and  immediately  after 
the  application  of  a  short  (~ 100  microseconds)  duration  stress  pulse.  This  investigation  is  conducted  by  sampling 
(at  10  microsecond  intervals)  the  attenuation  and  velocity  changes  of  a  high  frequency  (~  100  to  250  MHz)  ultra¬ 
sonic  wave,  while  the  stress  pulse  Is  propagating  through  the  specimen  under  study.  The  above  sampling  pro¬ 
cess  continues  for  a  period  of  several  milliseconds,  thus  providing  a  profile  of  dislocation  behavior  both  during  the 
stressing  and  immediately  after  the  removal  of  the  stress,  i.e.,  during  the  very  early  stages  of  recovery.  By  adjust¬ 
ing  the  stress  pulse  “shape”  it  is  also  possible  to  obtain  Information  on  the  behavior  of  dislocations  while  the  stress 
is  rising  and  falling.  At  the  end  of  the  initial  fast  sampling  period,  such  sampling  can  be  extended,  or  its  rate  can  be 
decreased  and  the  recovery  process  following  the  stress  pulse  application  can  be  monitored  over  an  extended 
period  of  time. 

Principal  Investigator:  C.  Elbaum  (Physics  and  Applied  Mathematics). 
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Personnel:  A.  Hikata  (Applied  Mathematics)  and  S.  F.  Hahn  (Physics). 
Supported  by  the  Materials  Research  LaboratoryINSF. 


Effects  of  Strain  Rate  and  Temperature  on  Plastic  Deformation  of  Steels 

A  comprehensive  experimental  study  involving  static,  dynamic  and  incremental  strain  rate  plasticity  experiments 
on  low  carbon  cold-rolled  and  hot-rolled  steels  at  temperatures  ranging  from  -150°C  to  + 1 10°C  has  now  been 
extended  for  the  hot-rolled  steel  to  -190°C.  An  unexpected  result  of  these  latter  experiments  is  an  inverse  strain 
rate  history  effect.  This  is  seen  in  the  fact  that  the  flow  stress  in  the  dynamic  portion  of  the  Incremental  strain  rate 
experiment  is  significantly  greater  than  the  flow  stress  resulting  from  loading  entirely  at  the  same  constant 
dynamic  strain  rate.  Furthermore,  this  inverse  effect  becomes  more  pronounced  as  the  testing  temperature  is 
decreased.  A  possible  explanation  would  involve  thermal  effects — the  isothermal  character  of  quasi-static  strain¬ 
ing  in  contrast  to  the  nearly  adiabatic  conditions  during  dynamic  straining.  A  technical  report  on  this  subject  is 
being  prepared.  An  oral  presentation  was  made  at  the  ASTM  Symposium  on  Mechanical  Testing  for  Deformation 
Model  Development  held  in  Bal  Harbour,  Florida  in  November  1980  and  publication  will  be  in  an  ASTM  Special 
Technical  Publication. 

Dr.  J.  Klepaczko  of  the  Institute  of  Fundamental  Technological  Research  of  the  Polish  Academy  of  Sciences  in 
Warsaw  has  been  a  visitor  at  Brown  since  April  1 , 1 980.  In  the  time  he  has  been  here  he  has  initiated  an  exten¬ 
sive  analysis  of  the  results  of  experiments  (some  of  which  were  conducted  at  Brown  or  at  the  Polish  Academy  of 
Sciences)  with  the  goal  of  clarifying  the  source  of  strain  rate  sensitivity,  strain  rate  history  effects  and  temperature 
history  effects  in  bcc  metals — especially  ferritic  steels  and  iron.  This  study  has  resulted  or  will  result  in  at  least 
three  technical  reports: 

1)  "Recent  Progress  in  the  Understanding  of  Strain  Rate  and  Temperature  History  Effects  in  bcc  Metals,"  with  J. 
Duffy,  presented  at  the  ASTM  Symposium  on  Mechanical  Testing  for  Deformation  Model  Development  at  Bal 
Harbour,  Florida  in  November  1980.  To  be  published  in  an  ASTM  Special  Technical  Publication. 

2)  "On  Strain  Rate  and  Strain  Rate  History  Effects  in  bcc  Metals  and  Ferritic  Steels." 

3)  "Relation  of  Thermally  Activated  Flow  in  bcc  Metals  to  Strain  Rate  and  Temperature  History  Effects." 

Experimental  results  including  published  results  in  the  literature  are  being  analyzed  on  the  basis  of  thermally  acti¬ 
vated  control  of  strain  rate  effects.  A  literature  survey  includes  a  search  for  evidence  of  strain  rate  history  and 
temperature  history  effects  in  ferritic  steels  and  other  bcc  metals. 

Principal  Investigator:  J.  Duffy  (Engineering). 

Personnel:  J.  Klepaczko  (Engineering),  R.  H.  Hawley  (Engineering),  G.  J.  LaBonte,  Jr.  (Engineering).  P.  Rush 
(Engineering),  A.  Rosakis  (Engineering),  and  A.  Morrone  (Engineering). 

Publications:  ‘On  the  Uniqueness  of  Representation  of  the  Stress  Field  of  Plane  Polygonal  Dislocation  Loops," 
A.  J.  Rosakis,  Scripta  Metallurgies  14  (1980),  1261-65. 

‘Strain  Rate  History  Effects  in  LIF  Single  Crystals  During  Dynamic  Loading  in  Shear,”  C.  Y.  Chiem  and  J.  Duffy, 
Brown  University  Technical  Report  NSF  ENG75- 18532/11,  July  1979.  Accepted  for  publication  in  Material  Sci¬ 
ence  and  Engineering. 

Supported  by  the  Materials  Research  LaboratoryINSF,  the  National  Science  Foundation,  and  the  Army 
Research  Office. 
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Rate  Theory 


Methods  previously  developed  for  the  study  of  quantum  rate  theory  in  a  symmetric  double-well  potential  have 
been  extended  to  the  asymmetric  case.  The  approach  involves  the  decomposition  of  the  real  stationary  states 
into  complex  right-  and  left-moving  states  and  the  determination  of  exact  transmission  coefficients  for  the  latter. 
Numerical  solutions  for  a  selection  of  representative  parameter  values  were  obtained  for  the  time  evolution  of 
states  which  corresponded  initially  to  a  Gaussian  wave  packet  localized  in  the  upper  well.  Both  resonant  and 
nenresonant  cases  were  studied  and  the  significance  of  the  transmission  coefficient  was  established  for  each. 

Principal  Investigator:  J.  H.  Weiner  (Engineering). 

Personnel:  S.  T.  Tse  (Engineering). 

Publication:  “Tunneling  in  Asymmetric  Double-Well  Potentials,"  J.  H.  Weiner  and  S.  T.  Tse,  J.  Chem.  Phys., 
forthcoming  (February  15,  1981). 

Supported  by  the  Materials  Research  LaboratoryiNSF  and  the  Gas  Research  Institute. 


Crystal  Elasticity 

Theoretical  calculations  of  the  elastic  moduli  of  crystals  are  generally  based  on  assumed  interatomic  force  laws 
which  in  turn  are  intended  to  reflect  the  electron  states  in  the  given  system.  In  this  work,  the  Hellmann-Feynman 
theorem  has  been  used  to  derive  a  formula  for  the  elastic  moduli  of  centro-symmetric  crystals  directly  in  terms  of 
the  electron  density  and  its  distortion  due  to  the  crystal  deformation.  One  consequence  of  this  formula  is  the 
observation  that  these  moduli  satisfy  the  Cauchy  relations  if  the  electron  density  deforms  according  to  the  same 
affine  transformation  which  governs  the  crystal  deformation. 

Principal  Investigator:  J.  H.  Weiner  (Engineering). 

Publication:  “Hellmann-Feynman  Theorem,  Elastic  Moduli,  and  the  Cauchy  Relations,”  J.  H.  Weiner,  submitted  to 
Phys.  Rev. 

Supported  by  the  Materials  Research  LaboratoryiNSF  and  the  Gas  Research  Institute. 


Plastic  Shearing  at  Very  High  Strain  Rates 

Thin  plate  specimens  (~0.4  mm  thick)  of  aluminum  sandwiched  between  two  high  strength  steel  plates  have 
been  subjected  to  shear  strain  rates  of  104  -  10*  s 1  by  means  of  pressure-shear  plate  impact  experiments.  In 
these  experiments,  a  flyer  plate,  inclined  to  the  launch  tube  axis,  impacts  a  parallel  target  plate  assembly  consist¬ 
ing  of  the  specimen  backed  by  a  high  strength  steel  plate.  By  using  the  pressure-shear  configuration  instead  of 
the  conventional  normal  impact  configuration,  it  is  possible  to  sustain  high  shear  strain  rates  until  unloading 
waves  arrive  from  boundaries  of  the  plates.  Preliminary  results  indicate  the  feasibility  of  using  such  experiments 
to  measure  the  flow  stress  at  strain  rates  of  105  s'1.  For  commercially  pure  aluminum  these  experiments  show  a 
marked  increase  in  flow  stress  over  the  values  reported  for  strain  rates  of  103  s'1. 

Principal  Investigator:  R.  J.  Clifton  (Engineering). 

Personnel:  C.  H.  Li  (Engineering). 

Supported  by  the  United  States  Army  Research  Office  and  the  Materials  Research  LaboratoryiNSF. 


Shear  Localization  in  Plastic  Deformation 


A  previous  quasi-static  analysis  of  the  stability  of  a  pure  shearing  deformation  has  been  extended  to  include  the 
effects  of  Inertia.  As  in  the  earlier  quasi-static  case,  the  analysis  includes  the  effects  of  strain  hardening,  heat 
conduction,  and  the  sensitivity  of  the  flow  stress  to  strain  rate  and  temperature.  Conditions  are  derived  for  the 
growth  of  infinitesimal  variations  from  a  homogeneous  deformation.  Explicit  stability  conditions  are  obtained  for 
both  the  long  and  short  wavelength  limit  as  well  as  for  the  special  case  in  which  the  material  is  modeled  as  a  fluid 
with  temperature  dependent  viscosity.  Stability  conditions  for  more  general  cases  require  determination  of  the 
roots  of  a  polynomial  of  either  degree  7  or  degree  4,  depending  upon  the  level  of  generality  retained. 

Principal  Investigator:  R.  J.  Clifton  (Engineering). 

Supported  by  the  United  States  Army  Research  Office  and  the  Materials  Research  Laboratory/NSF. 


Elastic-Plastic  Wave  Propagation 

The  situation  of  a  step  shear  load  which  moves  steadily  on  the  surface  of  an  elastic-plastic  half  space  at  a  speed 
exceeding  the  elastic  shear  wave  speed  of  the  material  has  been  examined.  The  orientation  of  the  shear  traction 
is  such  that  the  deformation  is  two-dimensional  antiplane  strain.  Two  different  representations  of  the  rate  inde¬ 
pendent  elastic-plastic  material  response  have  been  considered.  The  first  material  model  is  based  on  the  associ¬ 
ated  flow  rule  and  the  Mises  yield  condition  with  isotropic  hardening,  whereas  the  second  model  is  based  on  a 
particular  flow  theory  of  plasticity  which  represents  incremental  behavior  at  a  corner  of  the  instantaneous  yield 
surface.  Both  models  predict  the  same  response  under  the  same  proportional  loading.  The  stress  history  experi¬ 
enced  by  a  typical  material  particle  during  passage  of  the  load  step  has  been  determined,  and  the  variation  of  final 
strain  with  the  magnitude  of  the  load  step  has  been  calculated.  One  conclusion  resulting  from  comparison  of 
results  for  the  two  material  models  for  this  problem  is  that  yield  surface  vertex  formation  is  not  significant  in 
determining  the  final  states  of  plastic  strain. 

Principal  Investigator:  L.  B.  Freund  (Engineering). 

Personnel:  A.  S.  Douglas  (Engineering). 

Publication:  “Shear  Waves  in  an  Elastic-Plastic  Material  Due  to  a  Supersonic  Surface  Step  Load."  L.  B.  Freund 
and  A.  S.  Douglas,  Wave  Motion  2  (1980),  159-165. 

Supported  by  the  Materials  Research  Laboratory/NSF  and  the  National  Science  Foundation. 


Pressure-Shear  Waves 

Symmetric  impact  of  parallel  plates  inclined  relative  to  the  direction  of  approach  have  been  used  to  study  dynamic 
plastic  deformation  for  loading  along  non-proportional  stress  paths.  Further  experiments  on  6061 -T6  aluminum 
appear  to  confirm  that  slip  does  not  occur  at  the  impact  face  for  a  2:1  ratio  of  normal  to  transverse  components  of 
projectile  velocity,  provided  that  the  impact  surfaces  are  not  highly  polished.  A  series  of  experiments  on  a-tita- 
nium  plates  indicate  that  the  wave  profiles  are  similar  in  general  appearance  to  those  obtained  for  aluminum 
plates.  However,  more  detailed  comparisons  suggest  that  strain  rate  sensitivity  has  a  greater  effect  on  the  wave 
profiles  for  a-tftanium.  This  conclusion  follows  from  the  observation  that  simple  wave  solutions  cannot  describe 
the  observed  wave  profiles  as  well  for  a-titanlum  as  for  6061 -T6  aluminum. 


Principal  Investigator:  R.  J.  Clifton  (Engineering). 

Personnel:  A.  Gilat  (Engineering)  and  L.  Hermann  (Engineering). 

Publication:  “Pressure-Shear  Impact  of  6061 -T6  Aluminum,"  K.  S.  Kim  and  R.  J.  Clifton,  J.  Appl.  Mech.  47  (1980), 
11-16. 

Supported  by  the  Materials  Research  Laboratory  I NSF  and  the  National  Science  Foundation. 
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SECTION  2 


Fracture  of  Solids 


Introduction 


Studies  on  the  fracture  of  solids  encompass  a  wide  range  of  problem  areas  in  solid  mechanics  and  materials 
science.  One  area  given  major  emphasis  has  been  microscale  cavitation  processes  and  the  resulting  prediction 
of  macroscopic  fracture  ductility  in  terms  of  microstructure.  Parts  of  this  work  are  directed  to  elevated  tempera¬ 
ture  diffusive  creep  rupture,  with  consideration  of  grain  deformabiiity  and  non-equilibrium  cavity  shapes  (Nee- 
dleman,  Rice).  Other  work  focuses  on  lower  temperature  cavitation  processes,  e  g.,  by  a  grain  boundary  split¬ 
ting  mode  of  ductile  void  growth  in  spheroidized  steels  (Gurland,  Asaro),  and  on  the  effect  of  hydrogen  on 
accelerating  this  process  (Asaro).  Further  studies  have  used  microscale  models  for  ductile  void  growth  os  a 
basis  for  macroscopic,  dilational  plastic  constitutive  relations  and,  through  finite  element  and  other  numerical 
analysis  techniques,  have  led  to  assessments  of  effects  of  necking,  triaxiality  and  strain  rate  on  fracture  ductility 
(Needleman,  Saje). 

A  second  area  of  major  emphasis  is  that  of  plastic  bifurcation  and  instability,  particularly  in  the  form  of  flow  local¬ 
ization  into  a  shear  band  as  a  precursor  to  ductile  rupture.  Some  work  in  this  area  has  continued  the  association 
established  earlier  between  the  constitutive  description  of  double  slip  in  fee  crystals  and  shear  localization 
(Asaro,  Needleman).  Other  works,  directed  to  polycrystalline  materials,  employ  recently  developed  models  of 
vertex-like  yielding  as  a  basis  for  localization  and  bifurcation  (Needleman).  Studies  appropriate  to  wrinkling 
bifurcations  in  sheet  metal  forming  and  surface  bifurcations  in  non-homogeneously  deformed  bodies  (e  g., 
beams  in  finite  bending)  have  been  done  (Needleman,  Abeyaratne).  Also,  a  major  accomplishment  has  been 
the  development  of  finite  element  analysis  procedures  suitable  for  severely  localized  flow  zones  in  the  form  of 
shear  bands  (Tvergaard,  Needleman,  Lo). 

A  further  area  of  major  emphasis  is  on  the  dynamic  and  non-linear  analysis  of  crack  growth.  Crack  toughness 
under  high  rate  loading  of  steels  has  been  measured  and  compared  to  static  toughness  (Duffy).  Advances  in 
computational  techniques  for  dynamic  crack  propagation  have  been  made,  allowing  also  the  possibility  of  elas¬ 
tic-plastic  material  response  (Freund).  A  method  based  on  caustics  has  been  developed  for  incorporation  of 
small  scale  plasticity  effects  in  observational  studies  of  dynamic  crack  tip  fields  (Freund,  Duffy).  The  field  of 
stress  and  deformation  very  near  a  stably  growing  tensile  crack  in  an  ideally  plastic  solid  has  been  analyzed  in 
further  detail,  and  has  served  as  basis  for  a  predicted  form  of  a  ductile  crack  growth  resistance  curve  shown  to 
have  experimental  support  for  a  high  strength  steel  (Hermann,  Rice).  Further  advances  have  been  made  on  the 
analysis  of  crack  tip  fields  in  time-dependent  materials,  intended  to  model  response  under  high  temperature 
creep  conditions  (Kubo). 

Also,  major  advances  have  been  made  on  crack  analysis  under  non-linear  conditions  in  an  idealized  continuum 
model  of  fiber-reinforced  composite  materials  (Pipkin). 

A  variety  of  other  studies  in  solid  mechanics  and  materials  science  have  been  carried  out,  involving,  e  g., 
hydride  formation  and  failure  in  zircalloy  (Simpson),  fracture  mechanisms  in  two-phase  alloys  such  as  tungsten- 
carbide-cobalt  (Gurland),  oxidation  failures  in  reaction-bonded  nitride  (Richman),  cavitation  in  cobalt  titanium 
alloys  (Avery),  crack  tip  fatigue  processes  in  age-hardenable  aluminum  alloys  (Asaro),  shear  fracture  and  crus¬ 
tal  stress  distribution  processes  in  geomaterials  (Lehner,  Rice),  constitutive  descriptions  of  frictional  slip  (Rice), 
and  other  topics. 
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Fracture  of  Solids 


Individual  Contributions 


Analysis  of  Elastodynamic  Crack  Growth 

Research  has  continued  on  the  analysis  of  crack  growth  in  nominally  elastic  materials  at  rates  large  enough  for 
material  inertia  to  be  important.  Work  has  been  completed  on  the  following  projects: 

a)  The  dynamic  propagation  of  a  crack  in  an  anti-plane  shear  deformation  field  has  been  analyzed  by  second-or¬ 
der-accurate  finite  differences.  The  finite  difference  equations  have  been  obtained  by  integrating  the  dynamic 
linear  elastic  equations  of  motion  along  the  bicharacteristic  strips  in  four  perpendicular  directions  and  the  time 
axis  to  0(At3).  The  singularity  in  stresses  around  the  crack  has  been  calculated  by  performing  a  global  energy 
balance  on  a  small  region  containing  the  crack  tip  and  approximating  the  stresses  and  velocity  in  this  region  by  a 
one  term  asymptotic  expansion  about  the  crack  tip.  Results  for  stresses  and  stress  intensity  factor  have  been 
presented  for  a  semi-infinite  crack  propagating  steadily  in  an  infinite  strip,  from  which  errors  in  the  numerical  cal¬ 
culations  are  identified.  Cases  of  typical  non-steady  crack  propagation  in  an  infinite  strip  following  steady  propa¬ 
gation  have  also  been  considered 

b)  The  dynamic  numerical  analysis  of  linear  elastic  bodies  containing  a  crack  (stationary  or  propagating)  has 
been  carried  out  using  four  noded  isoparametric  finite  elements  with  central  differences  in  time.  The  crack  tip 
region  was  represented  by  the  one  term  asymptotic  expansion  of  the  near  tip  field  where  the  stress  intensity  is 
obtained  by  performing  an  energy  balance  on  a  small  region  containing  the  crack  tip.  Examples  of  both  Modes  I 
and  III  crack  propagation  have  been  compared  with  analytic  solutions. 

c)  The  problem  of  an  edge  crack,  propagating  with  constant  velocity  perpendicular  to  the  traction  free  boundary  of 
a  linear  elastic  half  space  has  been  solved  numerically  for  typical  self-similar  loading  cases.  A  limiting  case  of  uni¬ 
form  pressure  on  the  crack  faces  has  been  shown  to  have  a  stress  intensity  factor  close  to  that  tor  the  equivalent 
static  problem,  which  serves  as  a  check  on  the  method.  The  effect  of  stress  gradients  has  been  demonstrated 
and.  for  impact  loading,  it  is  shown  that  cracks  can  only  initiate  under  the  center  of  the  load  if  friction  is  present, 
and  then  only  at  very  low  crack  tip  speeds. 

d)  The  plane  strain  response  of  an  unbounded  elastic  body  containing  a  semi-infinite  crack  subjected  to  a  pair  of 
concentrated  forces  suddenly  applied  to  the  crack  faces  at  some  distance  from  the  crack  tip  has  been  deter¬ 
mined.  The  forces  act  on  opposite  faces  of  the  crack,  in  the  plane  of  the  crack,  and  in  the  same  direction.  An  exact 
solution  has  been  obtained  within  the  framework  of  linear  elastodynamics  using  a  fundamental  solution  obtained 
from  dynamic  elastic  dislocation  theory.  If  the  loading  is  quasi-staticalfy  applied,  then  the  stress  intensity  factor  is 
zero.  However,  if  the  loads  are  suddenly  applied,  the  stress  intensity  factor  varies  with  time  and,  for  a  short  time,  it 
takes  on  very  large  values.  As  time  becomes  large  compared  to  the  transit  time  of  a  Rayleigh  wave  from  the  load 
point  to  the  crack  tip,  the  stress  intensity  factor  decays  to  zero. 

Principal  Investigator:  L.  B.  Freund  (Engineering). 

Personnel:  P.  Burgers  (Engineering)  and  I.  S.  Abou-Sayed  (Engineering). 

Publications:  “Dynamic  Linear  Elastic  Crack  Propagation  in  Antiplane  Shear  by  Finite  Differences,"  P.  Burgers, 
International  Journal  of  Fracture  16  (1980),  261-74 

‘Numerical  Analysis  of  Dynamic  Linear  Elastic  Fracture  Mechanics  Problems  by  Finite  Elements,"  P.  Burgers  and 
L.  B.  Freund,  Brown  University  Technical  Report,  February  1979. 

“Dynamic  Growth  of  an  Edge  Crack  in  a  Half-Space,"  P.  Burgers  and  L.  B.  Freund,  International  Journal  of  Solids 
and  Structures  16  (1980),  265-74, 


“The  Stress  Intensity  Factor  Due  to  Parallel  Impact  Loading  of  the  Faces  of  a  Crack,'  I.  S.  Abou-Sayed, 
P  Burgers,  and  L.  B.  Freund,  Fracture  Mechanics,  ASTM  STP  700,  American  Society  for  Testing  and  Materials 
(1980),  164-73. 

Supported  by  the  Materials  Research  Laboratory /NSF,  the  National  Science  Foundation,  and  the  Office  of 
Naval  Research. 


Dynamic  Elastic-Plastic  Crack  Growth 

A  numerical  analysis  of  steady-state  crack  growth  in  an  elastic-ideally  plastic  material  under  small  scale  yielding 
conditions  has  been  performed.  Crack  growth  is  in  the  antiplane  shear  mode,  and  inertial  resistance  of  the  mater¬ 
ial  to  motion  has  been  included  explicitly.  The  numerical  procedure  is  based  on  the  finite-element  method.  It  has 
been  found  that,  with  increasing  crack  tip  speed,  the  active  plastic  zone  becomes  smaller  and  less  concentrated 
in  the  direction  of  crack  growth,  and  the  crack  tip  plastic  strain  singularity  becomes  weaker.  It  has  also  been  found 
that  the  rate  of  crack  opening  at  a  small  distance  behind  the  moving  crack  tip  is  insensitive  to  changes  in  crack  tip 
speed. 

A  study  was  also  completed  on  the  analysis  underlying  the  use  of  the  optical  method  of  caustics  in  fracture  pro¬ 
pagation  experiments  in  metals.  The  shadow  spot  data  which  are  obtained  in  using  the  method  of  caustics  to 
determine  experimentally  dynamic  stress  intensity  factors  are  usually  interpreted  on  the  basis  of  a  static  elastic 
crack  model.  In  this  work,  an  attempt  was  made  to  include  both  crack  tip  plasticity  and  inertial  effects  in  the  analy¬ 
sis  underlying  the  use  of  the  method  of  reflection.  For  dynamic  crack  propagation  in  the  two-dimensional  tensile 
mode  which  is  accompanied  by  a  Dugdale-Barenblatt  line  plastic  zone,  the  predicted  caustic  curves  and  corre¬ 
sponding  initial  curves  have  been  studied  within  the  framework  of  plane  stress  and  small  scale  yielding  condi¬ 
tions.  These  curves  have  been  found  to  have  geometrical  features  which  are  quite  different  from  those  for  purely 
elastic  crack  growth.  Estimates  have  been  made  of  the  range  of  system  parameters  for  which  plasticity  and  iner¬ 
tia  effects  should  be  included  in  data  analysis  when  using  the  method  of  caustics.  For  example,  it  is  found  that  the 
error  introduced  through  the  neglect  of  plasticity  effects  in  the  analysis  of  data  will  be  small  as  long  as  the  size  of 
the  region  near  the  crack  tip  which  corresponds  to  the  observed  caustic  is  more  than  about  twice  the  plastic  zone 
size.  Also,  it  is  found  that  the  error  introduced  through  the  neglect  of  inertial  effects  will  be  small  as  long  as  the 
crack  speed  is  less  than  about  20%  of  the  longitudinal  wave  speed. 

Pnncipal  Investigator:  L.  B.  Freund  (Engineering). 

Personnel:  A.  S.  Douglas  (Engineering),  A.  J.  Rosakis  (Engineering),  and  D.  M.  Parks  (Engineering). 

Publications:  “Dynamic  Steady  Antiplane  Shear  Crack  Growth  In  an  Elastic-Plastic  Material,'  A.  S.  Douglas, 
L.  B.  Freund,  and  D.  M.  Parks.  Proceedings  of  the  5th  International  Conference  on  Fracture,  Cannes,  France, 
1981. 

“The  Effect  of  Crack  Tip  Plasticity  on  the  Determination  of  Dynamic  Stress  Intensity  Factors  by  the  Optical 
Method  of  Caustics,'  A.  J.  Rosakis  and  L.  B.  Freund,  Brown  University  Technical  Report  N0001 4-0051/5,  Sep¬ 
tember  1980. 

Supported  by  the  National  Science  Foundation  and  the  Office  of  Naval  Research. 


Fracture  Initiation  In  Steel 

In  previous  fracture  initiation  experiments  with  two  low  carbon  steels  (AIS1 1018  cold-rolled  and  1020  hot-rolled) 
the  temperature  ranges  for  cleavage  fracture  under  quasi-static  and  dynamic  loading  conditions  were  estab¬ 
lished.  In  both  instances  microscopic  examination  of  the  specimens  after  testing  showed  that  the  fractures  on  the 
lower  shelf  initiate  at  carbide  plates  that  precipitate  along  the  ferrite  grain  boundaries,  and  that  this  result  holds 
true  for  both  quasi-static  and  dynamic  loading.  In  order  to  pursue  this  subject  further,  and  in  collaboration  with 
R.  J.  Asaro  of  our  material  science  group,  an  investigation  was  undertaken  to  determine  the  fracture  behavior  of  a 
similar  steel  with  (ewer  and  smaller  carbide  plates.  AISI  1020  hot-rolled  steel  was  modified  by  isothermal 
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transformation  to  significantly  reduce  in  size  and  number  the  carbide  plates  along  the  grain  boundaries.  All  other 
physical  and  mechanical  characteristics  were  unaltered  by  this  process  and  it  was  found  that  the  plane  strain  frac¬ 
ture  toughness,  KIC,  is  raised  about  40%  throughout  the  temperature  region  where  cleavage  fracture  predomi¬ 
nates,  but  is  virtually  unchanged  at  higher  temperatures  where  ductile  fractures  prevail.  A  preliminary  report  is 
being  prepared  while  further  work  is  being  accomplished  which  will  track  the  source  of  cracks  from  the  fracture 
surface  to  determine  if  the  source  is  a  carbide  plate. 

Principal  Investigator:  J.  Duffy  (Engineering). 

Personnel:  See  below  (Dynamic  Crack  Propagation). 

Publications:  See  below  (Dynamic  Crack  Propagation). 


Dynamic  Crack  Propagation 

A  combined  experimental  and  analytical  study  of  crack  propagation  including  development  of  a  means  for  deter¬ 
mination  of  the  dynamic  fracture  toughness  Kw  of  a  running  crack  in  steel  by  the  method  of  caustics  is  underway. 
The  experiment  will  utilize  a  wedge  loaded  double  cantilever  beam  specimen  and  a  high  speed  12-frame 
Crantz-Schardin  type  camera  system.  At  the  present  time,  we  have  photographed  a  number  of  single  caustics 
produced  by  a  running  crack  in  4340  steel;  construction  of  the  12-frame  camera  has  begun.  A  method  of  analyz¬ 
ing  the  caustics  produced  by  a  running  crack  has  been  developed. 

Principal  Investigator:  J.  Duffy  (Engineering). 

Personnel:  R.  H.  Hawley  (Engineering),  G.  J.  LaBonte,  Jr.  (Engineering),  P.  Rush  (Engineering),  A.  Rosakis 
(Engineering),  and  J.  C.  Lambropofis  (Engineering). 

Publications:  “The  Effect  of  Loading  Rate  and  Temperature  on  Fracture  Initiation  in  1020  Hot-Rolled  Steel," 
M.  L.  Wilson,  R.  H.  Hawley,  and  J.  Duffy,  Engineering  Fracture  Mechanics  13  (1980),  371-85. 

“Analysis  of  the  Optical  Method  of  Caustics  for  Dynamic  Crack  Propagation,"  A.  J.  Rosakis,  Engineering  Fracture 
Mechanics  13  ,1980),  331-47. 

“The  Effect  of  Loading  Rate  and  Temperature  on  the  Initiation  of  Fracture  in  a  Mild,  Rate-Sensitive  Steel," 
L.  S.  Costin  and  J.  Duffy,  Trans.  ASME,  Journal  of  Engineering  Materials  and  Technology  101  (1979),  258-64. 

Supported  by  the  Materials  Research  Laboratory! NSF,  the  National  Science  Foundation,  the  Office  of  Naval 
Research,  and  the  Army  Research  Office. 


Localization  and  Fracture  of  Solids 

We  have  investigated  numerically  the  possibility  of  localization  in  the  Plane  Strain  Tensile  Test  using  different 
constitutive  relations.  Localization  was  found  to  occur  at  reasonable  strain  levels,  consistent  with  experimental 
observations.  The  effect  of  constitutive  relation  on  the  stress  fields  near  the  crack  tip  was  also  studied  to  see  if 
classical  plasticity  theory  with  a  smooth  yield  surface  gives  reasonable  results  under  conditions  of  unloading  and 
reloading.  It  was  concluded  that  a  constitutive  relation  which  gives  a  different  response  for  changes  in  loading 
directions  does  not  affect  the  result  significantly  in  the  case  that  was  studied. 

Finally,  a  model  which  accounts  for  the  residual  plastic  deformation  at  the  flanks  of  a  crack  under  cyclic  loading  is 
used  to  offer  an  explanation  of  crack  growth  retardation  under  changes  in  the  level  of  loading. 

Principal  Investigator:  K.  K.  Lo  (Engineering). 

Publications:  “Flow  Localization  in  the  Plane  Strain  Tensile  Test,"  K.  K.  Lo,  V.  Tvergaard,  and  A.  Needleman,  to 
appear  in  JMPS. 
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‘Fatigue  Crack  Closure  Under  a  Step  Increase  Load,"  K.  K.  Lo,  to  appear  in  J.  Appl.  Mech. 

"Effect  of  Yield-Surface  Vertex  or  Crack  Tip  Fields  under  Steady  Mode  III,”  K,  K.  Lo,  to  appear  in  JMPS. 
Supported  by  the  Materials  Research  Laboratory/NSF. 


Failure  of  Fiber-Reinforced  Materials 

Plane  deformations  of  fibrous  and  fiber-reinforced  materials  were  studied  by  using  mathematical  models  in  which 
the  fibers  are  treated  as  inextensible.  Variants  of  the  theory  cover  nets  (1,2),  closely-woven  or  coated  fabrics  (3), 
and  stiff  sheets  with  fibers  in  two  orthogonal  directions  (4).  Traction  boundary-value  problems  lead  to  coupled 
non-linear  integral  equations  that  generally  have  more  than  one  solution  (5),  some  stable  and  others  unstable. 
This  general  theoretical  work  provides  the  background  for  applications  to  tearing  and  crack  problems  (2).  The 
theory  is  simple  enough  that  although  finite  deformations  are  considered,  it  may  be  possible  to  predict  dynamic 
crack  or  tear  trajectories.  Work  on  dynamic  crack  trajectories  in  small  deformations  is  in  progress. 

Principal  Investigator:  A.  C.  Pipkin  (Applied  Mathematics). 

Personnel:  T.  G.  Rogers  (Applied  Mathematics)  and  L.  Mannion  (Applied  Mathematics). 

Publications:  "Some  Developments  in  the  Theory  of  Inextensible  Networks,"  A.  C.  Pipkin,  Q.  Appl.  Math.  38 
(1980),  343-55. 

"Holes  in  Inextensible  Networks,"  T.  G.  Rogers  and  A.  C.  Pipkin,  QJMAM  (forthcoming). 

“Plane  Traction  Problems  for  Inextensible  Networks,"  A.  C.  Pipkin,  QJMAM  (forthcoming). 

“Finite  Plane  Stress  of  Stiff  Fiber-Reinforced  Sheets,”  A.  C.  Pipkin,  JIMA  (forthcoming). 

“Boundary  Integral  Equations  for  Inextensible  Materials,"  A.  C.  Pipkin,  J.  Eiast.  (forthcoming). 

Supported  by  the  Materials  Research  Laboratory/ NSF,  the  National  Science  Foundation,  and  the  Science 
Research  Council  (U.K.). 


Fracture  Initiation  Criteria  at  Brittle,  Second  Phaae  Particles 

The  role  of  zirconium  hydride  as  a  fracture  initiation  site  in  Zirconium-2. 5%  Nb  alloy  has  been  studied.  Experi¬ 
ments  utilizing  acoustic  emission  equipment  were  carried  out  to  identify  the  precise  conditions  of  stress  and/or 
strain  in  the  material,  necessary  to  cause  fracture  of  the  hydride  particles.  Results  showed  that  without  develop¬ 
ment  of  triaxial  stress  conditions,  fracture  of  the  hydride  particles  was  difficult,  and  they  had  little  influence  on  frac¬ 
ture  of  the  alloy  (for  the  low  hydrogen  concentrations  studied,  <35  ug/g).  On  the  other  hand,  hydrides  cracked 
readily  after  some  plastic  deformations  had  occurred  (ep<1%),  if  specimens  were  notched  to  produce  a  triaxial 
stress  state.  These  results  were  explained  in  terms  of  a  mechanism  in  which  cracks  were  nucleated,  in  the 
hydrides,  by  slip  bands  impinging  from  the  zirconium  matrix.  The  stability  of  these  cracks  was  dependent  on  the 
normal  stress  on  the  hydrides  such  that  Instability  was  favored  by  triaxiality. 


Principal  Investigator:  L.  A.  Simpson  (Engineering). 

Publications:  See  below  (Prediction  of  Ductile  Fracture  Conditions  Using  J-Resistance  Curves). 
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Prediction  of  Ductile  Fracture  Conditions  Using  J-Reslstance  Curves 

Raw  data  obtained  prior  to  the  author's  arrival  at  MRL  was  analyzed  to  produce  crack  growth  resistance  curves 
for  Zr-2.5%  Nb,  expressed  in  terms  of  the  J-integral.  These  curves  were  used  to  predict  the  critical  crack  length  in 
pressure  tubes  used  in  the  CANOU  nuclear  reactor.  A  strong  specimen  dependence  of  R-curve  shape  was 
observed  which  was  discussed  in  terms  of  current  theories. 

Initial  experiments  were  also  carried  out  on  HV130  steel  to  determine  J-resistance  curves  for  conditions  of  small 
scale  yielding  through  to  full  plasticity,  to  test  the  model  of  J.  R.  Rice.  Project  was  incomplete  at  the  end  of  the 
author's  appointment,  but  is  being  continued  by  R.  Asaro. 

Principal  Investigator:  L.  A.  Simpson  (Engineering). 

Publications:  ‘Expressions  for  Calculating  J-Resistance  Curves,"  L.  A.  Simpson,  Int.  J.  Fracture  16  (1980), 
R247-R249. 

■Effects  of  Specimen  Geometry  on  Elastic  Plastic  R-Curves  for  Zr-2.5%  Nb,”  L.  A.  Simpson,  Fracture  1981 ,  Pro¬ 
ceedings  of  the  5th  International  Conference  on  Fracture,  Cannes,  France,  April  1981 ,  to  be  published  by  Perga- 
mon  Press. 

"Criteria  for  Fracture  Initiation  at  Hydrides  in  Zr-2.5%  Nb  Alloy,”  L.  A.  Simpson,  submitted  to  Metallurgical  Trans¬ 
actions.  February  1981. 

Supported  by  the  Materials  Research  Laboratory INSF  and  Atomic  Energy  of  Canada  Limited. 


A  Combined  Macroscopic  and  Microscopic  Approval  to  the  Fracture  of  Metals 

A  systematic  study  of  void  nucleation  at  cementite  particles  in  spheroidized  steels  during  tensile  deformation  was 
completed.  The  favored  sites  for  void  nucleation  were  the  interfaces  of  large  particles  situated  on  grain  bounda¬ 
ries. 

A  criterion  for  void  nucleation  based  on  an  analysis  of  the  local  stresses  and  elastic  energy  storage  at  spheroid 
particles  in  a  plastically  flowing  matrix  was  developed.  A  critical  normal  stress  requirement  and  energy  balance 
consideration  provided  a  double  criterion  for  stable  void  nucleation. 

Contributions  to  the  stored  elastic  energy  in  the  vicinity  of  the  particles  due  to  applied  loading  and  local  plastic 
incompatibilities  were  considered  and  models  were  evaluated  for  the  present  alloy  systems  in  accordance  with 
the  proposed  nucleation  criterion. 

Also  completed  was  an  experimental  study  of  hydride  cracking  in  zirconium  alloys,  using  acoustic  emission  as  a 
cracking  indication. 

Principal  Investigator:  J.  Gurfand  (Engineering). 

Personnel:  L.  A.  Simpson  (Engineering),  J.  R.  Fisher  (Engineering),  and  H.  Stanton  (Engineering). 

Publications:  ‘Void  Nucleation  in  Spheroidized  Steels  During  Tensile  Deformation,”  J.  R.  Fisher,  Ph.D.  thesis, 
Brown  University,  June  1980. 

The  Strength  and  Fracture  of  Two-Phase  Alloys  •  A  Comparison  of  Two  Alloy  Systems,"  J.  Gurland,  in  Fracture 
79,  Proceedings  of  the  Conference  on  Fracture,  University  of  the  Witwatersrand,  Johannesburg,  South  Africa. 
1979,  p.  69. 

The  Effect  of  Alloy  Deformation  on  the  Average  Spacing  Parameters  of  Non-Deforming  Particles,’  J.  R.  Fisher 
and  J.  Gurland,  Metallurgical  Transactions,  in  press. 

"Characterization  of  Hydride-Induced  Slow  Crack  Growth  in  CANDU  Nuclear  Reactor  Pressure  Tubes,” 
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L  A.  Simpson,  MRL  Technical  Report,  in  preparation. 

Supported  by  the  Materials  Research  Laboratory  I NSF  and  the  U.S.  Department  of  Energy. 


Strength  and  Fracture  of  Two-Phase  Alloys 

A  study  of  the  micro-mechanics  of  two-phase  alloys  has  been  initiated.  The  work  was  focused  on  the  following 
aspects:  1)  a  quantitative  discussion  of  spatial  continuity  of  particulate  aggregates,  2)  consideration  of  the  role  of 
load  transfer  between  soft  and  hard  constituents,  and  3)  an  experimental  study  of  the  yield  and  fracture  of  dual¬ 
phase  s'eeis 

An  interesting  preliminary  conclusion  is  that  load  transfer  is  of  considerable  importance  as  a  strengthening  mech¬ 
anism  in  this  type  of  alloy  It  is  a  function  of  the  aspect  ratio  of  the  hard  particles  and  of  the  connectivity  of  the  hard 
constituent. 

Principal  Investigator:  J.  Gurtand  (Engineering). 

Personnel:  A.  Szewczyk  (Engineering)  and  H.  Stanton  (Engineering). 

Publication:  ‘An  Approximate  Method  for  the  Estimate  of  the  Contribution  of  Load  Transfer  to  the  Internal  Stress 
in  Dispersed  Particles,"  J.  Gurtand,  Scripta  Metallurgies  13  (1979),  967. 

Supported  by  the  Materials  Research  LaboratoryINSF  and  the  U.S.  Department  of  Energy. 


Microstructure  and  Fracture  Toughness  of  Cemented  Carbides 

Work  continues  on  the  theory  of  fracture  of  sintered  WC-Co  alloys.  The  initial  strain  energy  release  rate  is 
related  to  the  actually  observed  fracture  process,  namely  intra-  and  transgranular  fracture  of  WC.  interfacial 
decohesion  between  carbide  and  binder,  and  plastic  deformation  to  rupture  of  binder  ligaments. 

In  addition,  the  crystallographic  relations  between  contiguous  WC  grains  in  the  sintered  microstructure,  have 
been  studied  experimentally  and  theoretically.  The  theory  is  based  on  the  identification  of  coincidence  lattice  sites 
associated  with  low  energy  boundaries.  The  most  frequent  boundary  configurations  are  being  studied  by  meas¬ 
urement  of  dihedral  contact  angles  and  slip  line  directions  in  contiguous  WC  grains. 

Principal  Investigator:  J.  Gurland  (Engineering). 

Personnel:  J.  Hong  (Engineering),  R.  Deshmukh  (Engineering),  and  H.  Stanton  (Engineering). 

Publication:  “The  Fracture  Toughness  of  WC-Co  Two-Phase  Alloys  -  A  Preliminary  Model,"  M.  Nakamura  and 
J.  Gurland,  Metallurgical  Transactions  11A  (1980),  141. 

Supported  by  the  National  Science  Foundation. 


ShMr  Zoom  and  Fracture 

A  common  observation  in  certain  highly  deformed  ductile  solids  is  the  more  or  less  abrupt  appearance  of  bands  of 
highly  localized  deformation — shear  bands— often  leading  to  fracture  along  these  zones.  Analytical  studies  of  cer¬ 
tain  aspects  of  this  Issue  were  conducted. 

In  the  first  paper  (see  below)  the  appearance  of  localized  shear  zones  at  the  tip  of  a  pre-existing  crack  in  a  nonll- 
nearly  elastic  material  was  examined.  An  unsuccessful  attempt  was  made  to  repeat  this  study  in  the  case  of  a 
blunted  notch.  In  the  third  paper,  the  initiation  of  shear  bands  at  an  Imperfection  site  was  examined  considering 
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both  hyper-elastic  and  hypo-elastic  constitutive  relations. 

The  second  paper  is  unrelated  to  these  studies.  It  presents  the  solution  of  a  rather  elementary  boundary-value 
problem  in  the  theory  of  large  strain  plasticity. 

Principal  Investigator:  R.  C.  Abeyaratne  (Materials  Research  Laboratory). 

Publications:  “Discontinuous  Deformation  Gradients  Away  from  the  Tip  of  a  Crack  in  Anti-Plane  Shear," 
R.  Abeyaratne,  to  appear  in  Journal  of  Elasticity. 

“Finite  Elastic-Plastic  Deformation  of  a  Rotating  Hollow  Cylinder,”  K.  K.  Lo  and  R.  Abeyaratne,  submitted  for  pub¬ 
lication  to  the  Journal  of  Applied  Mechanics. 

“On  the  Emergence  of  Shear  Bands  in  Plane  Strain,'  R.  Abeyaratne  and  N.  Triantafyllidis,  to  appear  in  Interna¬ 
tional  Journal  for  Solids  and  Structures. 

Supported  by  the  Materials  Research  Laboratory/NSF. 


Single  Crystal  Constitutive  Laws 

Pierce.  Asaro,  and  Needleman  are  continuing  work  on  large  strain  single  crystal  deformation.  Assuming  physi¬ 
cally  plausible  forms  for  the  slip-system  hardening  moduli,  we  are  using  Asaro’s  constitutive  formulation  to  calcu¬ 
late  important  macroscopic  features  in  a  tensile  test.  For  the  three-dimensional  crystal,  the  effects  of  latent  hard¬ 
ening  on  lattice  rotations  and  overshoot  of  the  <  1 00-<  1 1 1  >  symmetry  boundary  have  been  quantitatively  analyzed 
(D.  Pierce,  M.Sc.  thesis).  In  addition,  critical  conditions  for  the  onset  of  shear  localization,  viewed  as  a  bifurcation 
from  a  homogeneous  deformation  state,  have  been  determined.  A  significant  result,  as  in  Asaro’s  two-dimen¬ 
sional  analysis,  is  the  admissibility  of  bifurcation  while  the  slip  systems  are  still  strain-hardening.  To  further  our 
study,  we  are  again  concentrating  on  the  two-dimensional  model,  but,  by  means  of  the  finite  element  method,  we 
hope  to  understand  better  the  effects  of  inhomogeneous  deformation  (especially  necking)  on  the  inception  of  a 
localized  mode.  This  more  complete  understanding  of  single  crystal  behavior  will  allow  us  next  to  develop  an 
appropriate  description  of  polycrystalline  deformation  at  finite  strains. 

Principal  Investigators:  R.  J.  Asaro  (Engineering)  and  A.  Needleman  (Engineering). 

Personnel:  D.  Pierce  (Engineering). 

Publication:  ‘Latent  Hardening  and  Overshoot  in  the  Deformation  of  Single  Crystals,"  D.  Pierce.  M.Sc.  thesis. 
Brown  University.  June  1980. 

Supported  by  the  Materials  Research  Laboratory/NSF. 


Flow  Localization  In  Porous  Solids 

Ductile  metals,  when  highly  deformed,  are  prone  to  the  localization  of  deformation  in  the  form  of  shear  bands. 
Often  such  metals  contain  second  phase  particles  which  may  crack  or  debond  from  the  surrounding  matrix  mater¬ 
ial  and  the  voids  thus  nucleated  then  grow  by  plastic  deformation.  A  rate  independent  constitutive  model  of  a 
porous  plastic  solid  developed  by  Gurson  has  been  employed  to  Investigate  the  influence  of  void  nudeation 
occurring  during  the  deformation  history  on  shear  localization.  Both  plastic  strain-controlled  and  stress -control  led 
nudeation  processes  are  simulated.  Two  deformation  histories  are  considered,  one  corresponding  to  uniaxial 
tension  and  the  other  to  plane  strain  tension.  The  enhanced  triaxiality  at  the  center  of  a  neck  is  simulated  by 
application  of  Bridgman  s  solution  for  the  stress  and  deformation  state  at  the  minimum  section  of  a  necked  bar. 

For  the  rate  independent  version  of  the  study,  the  enhanced  triaxiality  due  to  necking  is  found  to  have  great 
influences  on  flow  localization  in  uniaxial  tension,  the  necking  has  reduced  the  localization  strain  by  about  a  factor 
of  two,  but  the  necking  only  has  moderate  effects  on  flow  localization  in  plane  strain  tension.  The  destabilizing 


effect  that  arises  when  void  nucleation  is  stress-controlled  and  nucleation  occurs  over  a  narrow  range  of  stress  is 
illustrated.  Results  are  obtained  by  employing  parameter  values  representative  of  spheroidized  carbon  steels 
employed  in  a  recent  experimental  study  carried  out  by  Fisher  (Ph  D.  thesis.  Brown  University.  1980),  and  the 
predictions  of  the  model  are  in  good  agreement  with  experimental  observations. 

In  a  companion  study,  the  effect  of  plastic  strain-rate  sensitivity  on  localization  has  bean  analyzed.  The  parameter 
values,  representative  of  the  above-mentioned  steels  studied  by  Fisher  and  the  rate  sensitivity  of  1018  cold-rolled 
steel  reported  by  the  experimental  results  of  Costin  et  al.  are  employed  to  characterize  a  porous  solid,  and  the 
numerical  predictions  of  the  localization  strain  in  axisymmetric  tension  with  necking  are  in  good  agreement  with 
the  fracture  strains  obtained  experimentally  and,  therefore,  with  those  given  by  the  rate  independent  analysis. 
However,  the  present  numerical  predictions  of  localization  show  that  the  material  rate  sensitivity  raises  the  locali¬ 
zation  strain  in  plane  strain  tension  more  than  in  axisymmetric  tension.  The  ductility  ratio,  the  ratio  of  the  localiza¬ 
tion  strain  in  plane  strain  tensior  to  that  in  axisymmetric  tension,  is  found  to  increase  significantly  from  the  value 
of  the  corresponding  rate  independent  analysis. 

Principal  Investigators:  A.  Needleman  (Engineering)  and  J.  R.  Rice  (Engineering). 

Personnel:  J  Pan  (Engineering)  and  M.  Saje  (Engineering). 

Publication:  “Void  Nucleation  Effects  on  Shear  Localization  in  Porous  Plastic  Solids,”  M  Saje,  J.  Pan,  and 
A.  Needleman,  submitted  for  publication. 

Supported  by  the  Materials  Research  Laboratory /NSF  and  the  National  Science  Foundation. 


Plastic  Shear  Localization  from  Inhomogeneous  Deformation  States 

Localization  of  plastic  flow  in  narrow  shear  bands  is  frequently  observed  in  ductile  metals  subject  to  tensile  or 
compressive  loading.  The  growth  of  such  shear  bands  has  been  previously  analyzed  in  cases  where  the  stress 
state  outside  the  band  is  homogeneous  and  the  band  is  assumed  to  have  constant  width.  In  an  ongoing  series  of 
studies  the  initiation  and  propagation  of  localized  shear  deformations  in  solids  subject  to  non-uniform  states  of 
deformation  is  being  analyzed 

In  the  plane  strain  tensile  test,  the  deformations  are  observed  to  remain  essentially  homogeneous  up  to  the  max 
imum  load  point  after  which  a  diffuse  neck  develops.  Then,  for  structural  metals  the  deformation  mode  often  shifts 
to  one  involving  a  band  or  bands  of  high  localized  shear.  Previous  numerical  studies  have  employed  the  simplest 
smooth  yield  surface  theory  of  plasticity  and  have  found  continued  growth  of  the  diffuse  neck,  with  no  tendency 
exhibited  for  the  deformation  pattern  to  shift  to  one  involving  Intense  localized  shearing.  It  was  not  surprising  that 
these  numerical  studies  gave  no  indication  of  shear  band  formation  since  studies  of  material  instability  reveal  that 
the  classical  smooth  yield  surface  elastic-plastic  solid  is  quite  resistant  to  localization  of  deformation  into  a  shear 
band  (Rice,  p.  207  in  Proc  14th  IUTAM  Congress,  Delft,  1976).  One  feature  of  the  plastic  response  of  polycrys¬ 
tals  is  that  the  incremental  tensorial  relations  between  strain  and  stress  are  expected,  within  a  rate-independent 
model,  to  exhibit  a  "vertex"  structure.  This  means  that  the  direction  (in  an  appropriate  hyperspace)  of  the  strain 
increment  is  at  least  somewhat  dependent  on  the  direction  of  the  stress  increment,  in  contrast  to  classical  consti¬ 
tutive  models  of  the  Prandtl-Reuss-Mises  type.  This  feature  leads  to  localization  instabilities  at  sufficiently 
reduced  strain  hardening  levels. 

A  finite  element  analysis  of  the  plane  strain  tensile  test  for  a  solid  that  develops  a  vertex  on  Its  yield  surface  has 
been  carried  out.  The  phenomenological  J2  comer  theory  of  Christoffersen  and  Hutchinson  (J.  Mech.  Phys.  Sol¬ 
ids  27,  (1979),  465)  was  employed.  In  this  J2  comer  theory  the  instantaneous  moduli  are  associated  with  those  of 
a  finite  strain  generalization  of  J2  deformation  theory  for  nearly  proportional  loading.  For  increasing  deviation  from 
proportional  loading  the  moduli  stiffen  smoothly  until  they  coincide  with  the  linear  elastic  moduli  for  stress  rates 
directed  along  or  within  the  comer  of  the  yield  surface.  A  variety  of  initial  thickness  imperfections  were  also  con¬ 
sidered.  One  component  of  the  imperfection  was  always  taken  in  the  form  of  the  diffuse  necking  mode.  Additional 
shorter  wavelength  imperfections  were  included  in  some  of  the  computations. 

In  all  cases  the  numerical  calculations  exhibited  a  transition  from  the  diffuse  necking  deformation  pattern  to  one 
having  well  defined  shear  bands.  The  particular  shear  band  pattern  that  emerged  depended  on  the  initial 


imperfection.  Among  the  notable  features  of  the  numerical  results  were:  (i)  that  the  shear  bands  tended  to  propa¬ 
gate  along  the  orientation  revealed  by  a  ‘material  instability"  analysis  to  be  most  critical  for  bifurcation  into  the 
shear  band  mode  and  (if)  that  secondary  shear  band  patterns  emerged  after  a  well  developed  primary  shear  band 
pattern  had  formed. 

We  continued  our  study  of  the  development  of  flow  localization  from  nonhomogeneous  deformation  states  by 
considering  shear  oand  development  in  a  plate  subject  to  pure  bending.  One  motivation  for  focusing  on  this  prob¬ 
lem  stems  from  applications  to  various  metal  forming  processes  which  involve  bending  of  sheet  metal.  Bending 
beyond  a  certain  minimum  radius  of  curvature  relative  to  the  plate  thickness  can  lead  to  shear  fracture  at  the  sur¬ 
face.  The  possibility  of  bifurcation  away  from  the  cylindrically  symmetric  fundamental  state  of  pure  bending  has 
been  studied  by  Triantafyllidis  for  an  incompressible  nonlinear  elastic  solid.  The  critical  bifurcation  mode  is  a  sur¬ 
face  wave  mode,  which  occurs  first  in  the  compressive  region,  with  the  shortest  possible  wavelength  being  criti¬ 
cal.  Subsequently,  if  the  fundamental,  unbifurcated  state  is  thought  to  be  continued  beyond  the  first  critical  bifur¬ 
cation  point,  ellipticity  of  the  governing  equations  is  lost  at  the  surfaces.  Thus,  the  bifurcation  behavior  in  pure 
bending,  where  a  short  wavelength  mode  is  most  critical,  contrasts  with  that  in  plane  strain  tension,  where  the 
mode  with  the  longest  possible  wavelength  is  most  critical.  Triantafyllidis'  bifurcation  results  do  not  hold  directly 
for  an  elastic-plastic  solid  since  in  pure  bending  strongly  nonproportional  loading  takes  place  prior  to  bifurcation. 
However,  near  the  surface  regions  that  are  most  highly  stressed  and  thus  most  important  for  bifurcation  propor¬ 
tional  loading  does  occur. 

We  carried  out  a  numerical  study  of  shear  band  development  in  a  plane  strain  strip  subject  to  pure  bending 
employing  the  same  two  material  models  as  used  in  aur  plane  strain  tension  study.  In  some  of  the  calculations  the 
material  parameters  employed  in  this  study  were  identical  to  those  used  in  the  plane  strain  tension  analysis. 

The  shear  band  patterns  that  developed  for  the  J2  comer  theory  solid,  with  strain  hardening,  resembled  those 
associated  with  the  formation  of  a  plastic  hinge  in  an  ideally  plastic  solid.  The  location  and  number  of  these 
“hinges'  depended  on  the  specified  initial  thickness  inhomogeneity  but  the  pattern  of  shear  bands  forming  the 
hinge  was  similar  in  each  case  considered. 

In  pure  bending,  the  shear  bands  must  propagate  inward  against  an  adverse  de formation  gradient  so  that  the 
peak  straining  occurred  at  the  free  surface.  In  contrast,  in  plane  strain  tension,  due  to  diffuse  necking,  the  most 
intense  shearing  can  occur  in  the  interior  of  the  specimen.  In  the  pure  bending  case  we  found  that  intense  surface 
oscillations  occurred  where  the  shear  bands  met  the  free  surface  on  the  compressive  side  of  the  plate.  Our 
results  indicate  that  in  pure  bending  of  a  plate  made  of  a  material  prone  to  shear  band  development,  the  protu¬ 
berances  that  develop  from  these  oscillations  would  serve  as  likely  failure  initiation  sites. 

Principal  Investigator:  A.  Needleman  (Engineering). 

Personnel:  N.  Triantafyllidis  (Engineering),  K.  K.  Lo  (Engineering),  and  V.  Tvergaard  (Engineering). 

Publications:  ‘Bifurcation  Phenomena  in  Pure  Bending,”  N.  Triantafyllidis,  J.  Mech.  Phys.  Solids  28,  (1980),  221 . 

"Flow  Localization  in  the  Plane  Strain  Tensile  Test,"  V.  Tvergaard,  A.  Needleman,  and  K.  K.  Lo,  J.  Mech.  Phys. 
Solids,  in  press. 

“On  the  Development  of  Shear  Bands  in  Pure  Bending,"  N.  Triantafyllidis.  A.  Needleman,  and  V.  Tvergaard,  sub¬ 
mitted  for  publication. 

Supported  by  the  Materials  Research  Laboratory/ NSF  and  the  National  Science  Foundation. 


Wrinkling  In  th«  Swift  Cup  Twt 

We  analyzed  the  onset  of  wrinkling  in  the  Swift  cup  test,  which  is  used  to  determine  the  drawabillty  of  sheet  met¬ 
als.  In  this  test  a  flat  sheet,  held  at  its  periphery  between  a  dfebtock  and  a  blankholder,  is  drawn  into  a  cup  by  a  cyl¬ 
indrical  punch.  The  term  wrinkling  refers  to  a  failure  mode  involving  the  formation  of  buckles  (or  wrinkles)  in  the 
flange,  which  Is  the  portion  of  the  sheet  remaining  between  the  die  and  the  blankholder.  We  modelled  the  flange 
as  an  annular  plate  subject  to  axisymmetrtc  radial  tension  along  Its  Inner  edge.  A  numerical  solution  to  the 
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bifurcation  problem  governing  the  onset  of  wrinkling  was  obtained  by  employing  separation  of  variables  in  con¬ 
junction  with  a  finite  element  method.  The  sheet  material  was  characterized  both  by  the  simplest  flow  theory  of 
plasticity  and  by  the  simplest  deformation  theory  of  plasticity.  In  each  case  isotropy  was  presumed  to  hold  in  the 
plane  of  the  sheet  but  norma!  plastic  anisotropy  was  accounted  for.  The  restraint  provided  by  the  blankholder  was 
modelled  by  resting  the  annular  plate  on  a  linear  elastic  foundation. 

We  found  that  the  critical  drawing  displacement  for  the  onset  of  wrinkling,  with  a  blankholder  present,  increased 
with  increasing  plastic  anisotropy.  These  results  which  were  obtained  employing  deformation  theory  were  consis¬ 
tent  with  experimental  observations  of  Naziri  and  Pearce  (La  Metall.  Italiana,  No.  8  (1968),  727)  who  found  that 
the  blankholder  load  required  to  suppress  wrinkling  decreased  with  increasing  plastic  anisotropy.  Our  results  also 
indicated  that  if  an  attempt  were  made  to  investigate  the  effect  of  plastic  anisotropy  on  wrinkling  without  a  blank- 
holder.  the  effect  would  not  be  observed.  Previous  analyses  which  were  carried  out  within  a  framework  of  beam 
theory  or  which  restricted  attention  to  the  simplest  flow  theory  of  plasticity  did  not  reveal  the  dependence  of  the 
onset  of  wrinkling  on  the  plastic  anisotropy. 

Principal  Investigator:  A.  Needleman  (Engineering). 

Personnel:  N.  Triantafyllidis  (Engineering). 

Publication:  “An  Analysis  of  Wrinkling  in  the  Swift  Cup  Test,"  N.  Triantafyllidis  and  A.  Needleman,  J.  Eng.  Mat. 
Tech.  102  (1980),  241-48. 

Supported  by  the  National  Science  Foundation. 


Effects  of  Hydrogen  on  Fracture 

Cialone  and  Asaro  have  recently  completed  an  experimental  study  Of  the  effects  of  hydrogen  on  ductile  fracture  in 
plain  carbon  steels.  Results  on  two  steels,  viz.  a  1 01 7  and  a  0.1 5  wt  percent  carbon  steel  have  shown  that  hydro¬ 
gen  affects  both  the  process  of  void  initiation  at  cementite  particles  and  void  growth  leading  to  a  significant  loss  in 
ductility.  The  results  are  consistent  with  the  view  that  hydrogen  reduces  the  cohesive  strength  of  the  particle-ma¬ 
trix  interface  thereby  facilitating  decohesion.  Void  growth,  which  appears  to  occur  to  a  large  extent  along  grain 
boundaries,  is  also  accelerated  by  hydrogen  and,  again,  the  results  suggest  that  hydrogen  has  the  effect  of  reduc¬ 
ing  the  cohesion  of  the  interface,  in  this  case  grain  boundaries.  These  results  are  consistent  with  our  earlier 
results  on  a  1045  plain  carbon  steel.  In  these  lower  carbon  steels,  however,  we  observed  a  transition  in  the  frac¬ 
ture  mode  with  increasing  hydrogen  content  from  ductile  void  initiation  and  growth  to  transgranular  cleavage.  The 
tendency  for  cleavage  was  found  to  be  dependent  upon  microstructure.  In  fact,  a  good  correlation  was  found 
between  the  total  amount  of  interface  (grain  boundary,  particle-matrfx,  etc.)  and  the  susceptibility  to  cleavage; 
less  interlace  was  associated  with  an  increased  occurrence  of  cleavage. 

Principal  Investigator:  R.  J.  Asaro  (Engineering). 

Personnel:  H.  Cialone  (Engineering). 

Publications:  “Hydrogen  Assisted  Fracture  of  Spheroidized  Plain  Carbon  Steels,”  H.  Cialone  and  R.  J.  Asaro,  to 
appear  in  Met.  Trans.,  1981. 

‘Hydrogen  Effects  in  the  Fracture  of  Spheroidized  Plain  Carbon  Steels,”  H.  Cialone  and  R.  J.  Asaro,  to  appear  in 
the  Proceedings  of  the  3rd  International  Conference  on  the  Effects  of  Hydrogen  on  Material  Behavior,  ed. 
A.  W.  Thompson  and  I.  M.  Bernstein,  ASM,  1981. 

Supported  by  the  U.S.  Department  of  Energy. 


n 


Oxidation  and  Fracture  Characteristics  of  Reaction  Bonded  Silicon  Nitride  (RBSN) 


The  application  of  high  strength  ceramic  materials  for  high  temperature  turbine  engines  and  components  not  only 
requires  that  these  normally  brittle  materials  have  strength,  toughness,  and  creep  resistance  at  elevated  tempera¬ 
tures,  but  also  that  these  materials  stand  up  to  the  hot  corrosive  and  erosive  atmospheres  within  the  turbine 
engines  and  that  any  oxidation  suffered  by  the  ceramic  materials  does  not  result  in  a  degradation  of  these  proper¬ 
ties.  Oxidation  studies  of  RBSN  have  been  conducted  to  determine  microstructurally  the  path  of  oxidation  in  the 
various  phases  and  microconstituents  present  in  the  material  when  used  in  turbine  engine  applications.  The  oxi¬ 
dation  of  unreacted  silicon,  the  oxidation  of  silicon  nitride,  and  the  penetration  of  grain  boundaries  and  surface 
flaws  by  these  brittle  oxides  leads  to  the  creation  of  Griffith  flaws  which  were  not  present  when  the  part  was  fabri¬ 
cated.  The  continued  utilization  of  ceramic  turbine  engine  components  of  RBSN  is  dependent  on  understanding 
the  effect  of  microstructural  oxidation  and  crack  formation  and  growth  and  also  upon  the  development  of  techni¬ 
ques  and  processes  for  rendering  the  surfaces  of  these  materials  oxidation  resistant  or  retardant.  Such  treat¬ 
ments  are  being  investigated  at  the  present  time  with  the  view  that  a  new  generation  of  RBSN  will  be  fabricated 
which  will  be  resistant  to  oxidation  and  hot  corrosion  and  which  will  allow  greater  turbine  engine  life. 

Principal  Investigator:  M.  H.  Richman  (Engineering). 

Personnel:  M.  Magida  (Engineering),  O.  Gregory  (Engineering),  and  J.  Fogarty  (Engineering). 

Publication:  "Effect  of  Processing  Parameters  on  Reaction  Bonding  of  Silicon  Nitride,"  M.  H.  Richman,  to  appear 
as  NASA-Lewis  Technical  Report. 

Supported  by  the  Materials  Research  Laboratory /NSF  and  the  National  Aeronautics  and  Space  Administration 
(Lewis). 


Shear  Rupture  of  Geological  Materials 

Work  on  brittle-elastic  stress-strain  behavior  of  compressed  rock  by  Kachanov  was  based  on  a  model  in  which 
frictional  sliding  on  microcracks,  accompanied  by  the  growth  of  secondary  (tensile)  cracks,  is  the  dominant 
mechanism  of  inelasticity.  An  approximate  theoretical  description  of  branched  frictional  sliding  cracks,  which  took 
into  account  uplifting  at  asperities,  was  constructed.  It  was  found  that  most  of  the  features  of  pre-peak  macro¬ 
scopic  behavior,  before  extensive  crack  interaction,  could  be  reasonably  well  described  by  the  model.  The  frac¬ 
ture  toughness  KIC  inferred  from  the  model,  based  on  observations  of  the  onset  of  dilatant  deformation  in  various 
rock  tests,  is  relatively  constant  and  has  values  in  the  range  of  those  determined  direc'iy  from  macro-scale  crack 
propagation  fracture  tests.  Time-dependent  crack  growth  (stress  corrosion  cracking)  was  also  considered  in  the 
framework  of  the  model,  as  a  source  for  creep  of  compressed  rock  in  the  elastic-brittle  range. 

Experimental  studies  of  sliding  friction  on  smooth  rock  surfaces  by  Ruina,  in  collaboration  with  Dieterich  at  the  lat¬ 
ter’s  laboratory  at  U.S.G.S.,  Menlo  Park,  provided  a  number  of  significant  results.  It  was  observed  that  for  steady 
slip  (at  fixed  normal  stress)  of  sandwich-type  quartzite  specimens  the  friction  coefficient  p.  approached  (with  ongo¬ 
ing  slip)  a  decreasing  function  of  sliding  velocity,  V*.  which  was  independent  of  previous  slip  history.  However,  if 
the  slip  velocity  was  suddenly  changed  from  V*  (at  which  steady  state  has  been  achieved)  to  another  value  V,  the 
sudden  change  Ap.  in  the  friction  coefficient  was  positive  for  V>V*  and  negative  for  V<V*.  That  is,  the  sudden 
change  in  frictional  resistance  occurs  in  a  sense  opposite  to  the  long-term  change.  On  a  plot  of  p  versus  slip  dis¬ 
tance,  the  result  for  a  transition  from  speed  V,  to  V2  is  very  nearly  the  mirror  image  of  that  for  a  transition  from  V2 
to  V,. 

The  total  friction  coefficient  p  can  be  written  as  p  =  p(  +  pt,  where  p(  and  p§  are  long-term  and  short-term 
components  respectively.  The  long-term  component  pf  of  frictional  resistance  was  found  to  approach  the 
steady  state  resistance,  given  by,  say,  p  =  g(V),  at  a  rate  that  depends  on  the  difference  between  p(  and  g(V). 
This  result  can  be  fit  wel)  by. 

d^/dt  =  -( V/df)(ji,-g(V)l, 

where  d,  is  a  characteristic  slip  distance  (d(~5.2pm  for  the  quartzite  surface  examined).  Short-term  decay 
which  is  responsible  for  the  sudden  increase  in  friction  associated  with  Ap  is  also  fit  tolerably  by  an  exponential 
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decay  of  the  form: 


d[ns-f(V)]/dt  =  -<v/d>, 

where  the  characteristic  slip  distance  ds=0.32p.m  is  much  smaller  than  d(,  and  f(V)  represents  a  direct  velocity 
dependence  associated  with  the  sudden  increases  in  friction  with  rapid  velocity  increases  mentioned  above.  The 
two  differential  equations  for  and  ns,  together  with  the  steady  state  resistance  function  g(V)  and  the  direct 
velocity  dependence  function  f(V)  seem  to  provide  an  adequate  constitutive  model  for  the  results  so  far  observed. 
Over  the  range  of  velocities  tested  (approximately  0.01  to  1  p.m/sec)  these  two  functions  are  well  described  as 
linear  functions  of  log  V. 

Work  has  also  been  completed  on  the  mechanics  of  elastic  lithosphere/viscoelastic  asthenosphere  coupling  in 
time-dependent  stress  redistribution  along  rupturing  plate  boundaries.  A  generalized  Elsasser-type  plate  founda¬ 
tion  model  allowing  for  a  Maxwellian  viscoelastic  asthenosphere.  formulated  earlier  by  Rice,  has  been  analyzed 
by  Lehner.  Li.  and  Rice  for  spatial  and  temporal  stress  alterations  associated  with  suddenly  introduced  or  migrat¬ 
ing  (crack-like  or  dislocation-like)  ruptures  along  finite  segments  of  transform  or  subduction-type  plate  bounda¬ 
ries.  The  stress  which  is  shed  by  a  great  earthquake  onto  the  asthenosphere  is  gradually  transferred  back  to  the 
lithosphere  by  a  relaxation  process.  Accordingly,  it  was  found  that  the  post-seismic  increases  in  stress  outside  the 
rupture  zone  might  exceed  the  coseismic  elastostatic  stress  alterations  several  times  and  amount  in  magnitude  to 
an  appreciable  fraction  of  the  stress  drop  on  the  fault.  For  example,  in  the  distributed  dislocation  model  with  fixed 
slip  distributions  due  to  a  uniform  stress  drop  of  100  bars  on  the  ruptured  fault,  the  stress  alteration  at  one  lithos¬ 
pheric  distance  ahead  gradually  increases  to  about  1 8  bars  at  a  rate  comparable,  for  several  years  after  the  rup¬ 
ture.  to  that  due  to  plate  motion  corresponding  to  repeat  times  of  100  years.  Such  stress  transfer  may  contribute 
significantly  to  the  triggering  of  large  earthquakes  in  adjacent  segments  of  the  plate  boundary.  When  post-seismic 
slippage  along  the  fault  is  unimportant,  the  stress  on  the  fault  will  build  up  in  time  due  to  the  same  mechanism  of 
asthenosphere  relaxation.  The  spatial  and  temporal  characteristics  of  this  post-seismic  stressing  were  found  to 
be  quite  distinct  for  points  located  on  the  rupture  and  beyond  it,  suggesting  characteristic  patterns  of  prolonged 
aftershock  activity  for  locations  inside  and  outside  a  new  rupture  zone. 

Migration  of  seismic  activity  along  a  plate  boundary  has  been  modelled  by  a  steadily  advancing  zone  of  stress 
drop.  The  altered  stress  state  propagating  ahead  of  this  zone  is  found  to  correspond  very  closely  to  an  essentially 
unrelaxed  state  of  the  asthenosphere  at  typical  propagation  speeds  of  100  km/yr.  When  arrested  by  a  barrier,  a 
disturbance  of  this  kind  can  transfer  significant  stresses  to  the  barrier  from  a  then  relaxing  asthenosphere  as 
described  above,  and  could  thus  act  as  a  delayed  trigger  to  a  barrier  breaking  event. 

Principal  Investigator:  J.  R.  Rice  (Engineering). 

Personnel:  F.  K.  Lehner  (Engineering),  R.  Dmowska  (Geological  Sciences),  V.  C.  Li  (Engineering),  A.  Ruina 
(Engineering),  and  M.  Kachanov  (Engineering). 

Publications:  “The  Mechanics  of  Earthquake  Rupture,”  J.  R.  Rice,  in  Physics  of  the  Earth  s  Interior  (Proc.  Int. 
School  of  Physics  “Enrico  Fermi,"  Course  78),  ed.  A.  M.  Dziewonski  and  E.  Boschi,  Italian  Physical  Soc./North 
Holland  (1980),  pp.  555-649. 

Contribution  to  discussion  on  “Outstanding  Problems  of  Geodynamics:  Mechanisms  of  Faulting,”  J.  R.  Rice,  ibid., 
pp.  713-16 


“Microcrack  Model  for  Rock  Inelasticity,”  M.  L.  Kachanov,  Ph.D.  thesis,  Brown  University,  1980. 
Supported  by  the  National  Science  Foundation  (Geophysics  Program)  and  the  U.S.  Geological  Service. 


Elastic-Plastic  Analysis  of  Ductils  Crack  Growth 

Work  on  crack  growth  in  ductile  metals  has  been  based  on  a  combination  of  analytical  work  on  the  structure  of 
elastic-plastic  strain  field  singularities,  numerical  finite  element  studies,  and  experimentation  on  plane-strain 
mode  growth  in  steels. 
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An  experimental  investigation  carried  out  by  Hermann  and  Rice  assessed  the  accuracy  of  a  stable  crack  growth 
criterion  proposed  on  the  basis  of  earlier  theoretical  work  by  Rice,  Drugan,  and  Sham  These  theoretical  results 
combined  an  approximate  asymptotic  analysis  with  numerical  finite  element  studies  to  determine  the  stress  and 
deformation  fields  very  near  the  tip  of  a  growing  crack  in  an  elastic-ideally  plastic  material.  Knowledge  of  these 
near-tip  fields  allowed  statement  of  explicit  conditions  for  variation  of  J  with  crack  length,  at  small  scale  yielding,  to 
meet  a  fracture  criterion  which  hypothesized  a  self-similar  deformation  state  near  the  crack  tip  during  growth.  The 
experiments  were  an  extension  of  work  by  Odegaard  and  Asaro  on  high  strength  4140  steel,  now  heat  treated  to 
a  relatively  brittle  condition  to  facilitate  significant  stable  crack  growth  before  conditions  of  general  yielding  were 
met  This  was  done  because  the  theory  is  fully  developed  only  for  crack  growth  with  well-contained  yielding  The 
work  consisted  of  testing  four  pre-fatigued  compact  tension  specimens,  continuously  recording  the  load  P  and  the 
opening  displacement  1  along  the  load  line,  for  crack  growth  from  the  small  scale  yielding  range  up  through  the 
general  yielding  configuration  (representing  a  substantial  amount  of  stable  crack  growth,  eg.  4  to  9  mm  before 
general  yield)  A  differential  compliance  technique  was  used  to  measure  amounts  of  stable  crack  growth;  total 
compliance-estimated  growth  differed  from  visually  observed  growth  at  the  end  of  each  test  by  less  than  7%  in  the 
worst  case 

Curves  showing  the  experimentally-determined  dependence  of  the  "deformation  theory”  value  of  the  J-integral, 
J  on  crack  length  a  were  plotted  for  each  specimen  and  compared  with  the  theoretical  curves  of  Jd  versus  a 
obtained  by  integrating  the  expression  for  dj  da  given  by  the  small-scale  yielding  theory  of  Rice,  Drugan,  and 
Sham: 

dJ  da  -  (p  <»)(<r2  E)  in  (f><eR). 

Here  E  is  Young  s  modulus,  <rQ  is  the  tensile  yield  strength,  J  is  the  far-field  value  of  the  J-integral,  a  is  the  crack 
length,  e  is  the  natural  logarithm  base.  0  =  5.1  for  v  -  0.3;  correlation  with  the  finite  element  results  for  small  scale 
yielding  suggested  u  0  65  and  R  0  23  EJ Ui2Q.  The  parameter  p  governs  the  intensity  of  near-tip  crack  opening, 
and  is  a  material-dependent  constant  tor  the  hypothesized  criterion  of  growth  with  a  constant  near  tip  profile  The 
constant  was  chosen  to  provide  the  best  fit  of  the  theoretical  curves  with  the  experimental  ones  Proceeding  in 
this  manner,  good  theoretical  to  experimental  correlation  was  found  over  the  entire  range  of  stable  crack  growth 
prior  to  general  yielding. 

A  recently  completed  analytical  investigation  by  Rice  develops  a  procedure  for  exact  elastic-ideally  plastic  asymp¬ 
totic  analysis  of  fields  at  tips  of  growing  cracks  under  various  conditions.  Previous  plane  strain  solutions  have 
been  based  on  the  approximation  that  the  yield  condition  reduces  to  a  planar  form  of  constant  resolved  shear 
stress  in  all  plastically  active  sectors  near  the  crack  tip,  arguing  that  since  this  is  known  to  be  true  in  plastic  sec¬ 
tors  of  singular  strain,  it  is  approximately  correct  for  all  plastic  sectors.  The  new  formulation  treats  the  yield  condi¬ 
tion  as  fully  three-dimensional  in  non-singular  plastic  sectors,  showing  that  for  an  exact  solution  the  stress  fields 
in  such  plastic  sectors  are  not  statically  determinate.  A  compatibility  condition  is  required,  and  formulated,  to 
determine  the  full  set  of  equations  which  govern  the  asymptotic  field.  This  study  also  extends  the  concepts  devel¬ 
oped  previously,  for  asymptotic  analysis  of  crack  problems  in  specific  limited  cases,  to  crack  growth  with  arbitrary 
combinations  of  tensile  and  shear  loading  modes  in  materials  with  arbitrary  anisotropy  of  elastic  moduli  and 
(ideally)  plastic  yield  condition.  Among  the  results  of  the  study,  centered  fan  sectors  of  singular  straining  at  the 
crack  tip  are  shown  to  be  a  general  feature  of  the  large  class  of  material  response  treated,  as  are  also  the  plastic 
strain  singularity  types  associated  with  stationary  and  growing  cracks. 

A  final  finite  element  solution  for  a  stably  growing  crack  under  small-scale  yielding  and  plane  strain  conditions  in 
an  elastic-ideally  plastic  solid  was  initiated  by  T-L.  Sham  in  an  effort  to  improve  the  accuracy  with  which  the 
undetermined  parameters  arising  in  the  asymptotic  crack  tip  analysis  are  known.  The  procedure  is  based  on  the 
minimization  of  a  modified  energy  functional,  written  in  the  spirit  of  Hilton  and  Hutchinson  (Engrg.  Frac.  Mech., 
1971)  for  incremental  plasticity  under  small  displacement  gradient  assumptions.  In  this  procedure,  the  field  quan¬ 
tities  in  the  elastic  region  away  from  the  crack  tip  are  represented  by  series  convergent  at  large  distances  from  the 
tip  in  terms  of  unknown  coefficients.  The  leading  term  in  the  series  is  characterized  by  a  specified  stress  intensity 
K;  this  term  is  the  one  for  which  stresses  decay  as  r  ’ 2  The  near-tip  region,  where  the  material  response  is  elas¬ 
tic-plastic,  is  discretized  into  finite  elements  in  the  usual  manner.  The  nodal  displacement  increments  along  the 
boundary  of  these  two  regions  are  taken  to  be  those  given  by  the  series  representation,  and  the  discrete,  incre¬ 
mental  equilibrium  equations  are  obtained  through  the  minimization  of  the  modified  energy  functional.  The  advan¬ 
tage  of  this  procedure  is  that  the  small  scale  yielding  conditions  are  always  met  irrespective  of  the  size  of  the  plas¬ 
tic  zone  inside  the  finite  element  mesh,  and  we  can  concentrate  all  the  available  finite  element  degrees  of 
freedom  (which  are  limited  by  the  computer  capability)  in  the  plastic  zone  and  region  nearby. 


A  parallel  finite  element  investigation  was  begun  by  Sham,  Tse,  and  Pan  on  the  generation  of  a  stable  crack 
growth  solution  which  progresses  from  small  scale  yielding  to  general  yield  conditions  for  the  deeply  cracked 
bend  specimen  geometry.  This  geometry  and  loading  mode  correspond  closely  to  our  experimental  specimens, 
and  the  work  has  the  dual  purposes  of  showing  how  the  parameters  which  are  undetermined  by  the  asymptotic 
analyses  vary  through  the  range  from  small  scale  yielding  to  general  yield,  and  also,  therefore,  of  permitting  com¬ 
parison  between  the  theoretical  growth  equation  and  actual  experimental  results  all  the  way  up  to  and  through  the 
general  yielding  range. 

Work  on  creep  crack  analysis  was  also  in  progress  Kubo  investigated  the  nature  of  near-tip  fields  of  stationary 
cracks  in  solids  characterized  by  the  Bailey-Orowan  viscoplastic  model: 

r  =  FS"  V-a  ) 
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a  =  (Cf  /a“)-Dc.n  fMa 
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where  a  is  the  “internal  stress"  parameter,  n.  (1,  C,  D,  and  F  are  material  constants,  and  «  and  S  are  defined  as 


S  —Aa  -  a  )(a'  -  a  )  . 

This  model  can  account  for  creep  recovery,  thought  to  be  one  of  the  causes  of  high  crack  growth  rates  under 
cyclic  or  variable  loadings  at  high  temperatures.  One  important  result  of  the  work  thus  far  is  the  determination  that 
HRR  singular  fields  prevail  in  the  vicinity  of  the  crack  tip  for  any  combination  of  material  parameters.  Kubo  is  ana¬ 
lyzing  the  time  dependence  of  their  amplitude,  in  response  to  step  loading. 

Principal  Investigator:  J.  R.  Rice  (Engineering). 

Personnel:  W  J.  Drugan  (Engineering),  L.  Hermann  (Engineering).  T-L.  Sham  (Engineering),  H.  Pan  (Engineer¬ 
ing),  and  S.  Kubo  (Engineering). 

Publications:  "Comparison  of  Theory  and  Experiment  for  Elastic-Plastic  Plane  Strain  Crack  Growth,"  L.  Hermann 
and  J.  R.  Rice,  Metal  Science  14  (1980),  285-91 

“Elastic-Plastic  Crack  Growth,"  J.  R.  Rice,  in  Mechanics  of  Solids:  The  Rodney  Hill  60th  Anniversary  Volume,  ed. 
H.  G.  Hopkins  and  M.  J.  Sewell,  Pergamon  Press,  Oxford,  in  press. 

Supported  by  the  Materials  Research  Laboratory/ NSF  and  the  U.S.  Department  of  Energy. 


Elevated-Temperature  Creep  Cavitation  of  Grain  Boundaries 

Under  creep  conditions,  polycrystalline  solids  often  rupture  prematurely  due  to  the  diffusive  growth  and  coales¬ 
cence  of  grain  boundary  voids.  Driven  by  stress  gradients,  matter  diffused  from  the  void  surfaces  is  deposited 
along  the  grain  boundaries  Previous  analyses  of  the  problem  have  mainly  been  based  on  the  assumption  that 
the  grains  adjoining  the  cavitating  boundary  separate  rigidly.  Our  recent  studies  have  investigated  the  effect  of 
grain  deformability.  Needleman  and  Rice  have  analyzed  the  growth  of  cavities  along  grain  interlaces  by  the  com¬ 
bined  processes  of  grain  boundary  diffusion  and  plastic  dislocation  creep  in  the  adjoining  grains.  The  formulation 
of  the  problem  and  the  finite  element  scheme  employed  are  based  on  a  variational  principle,  which  governs  the 
combined  processes  described  above. 

The  principal  effect  is  that  creep  deformability  of  the  grains  allows  matter  diffused  from  the  cavity  walls  to  be 
accommodated  by  local  separations  of  the  adjoining  grains  in  the  vicinity  of  the  cavity,  thus  shortening  the  effec¬ 
tive  diffusion  path  length  and  allowing  high  values  of  the  local  matter  flux.  In  turn,  this  results  in  greater  rates  of 
cavity  growth  than  would  be  the  case  if  either  grain  boundary  diffusion  or  plastic  creep  flow  acted  alone. 


The  interaction  between  plastic  creep  and  diffusion  is  shown  to  be  characterized  by  a  parameter  L,  which  has 
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dimension  of  length  and  which  is  a  function  of  material  properties  (which  arise  in  describing  g.b.  diffusion  and  dis¬ 
location  creep),  temperature  and  applied  stress.  L  decreases  with  increasing  temperature  and  stress  and  (for  var¬ 
ious  pure  metals  stressed  to  10  3  x  shear  modulus  at  0  5  TM)  L  has  values  in  the  range  0  25  to  25  |im  The  con¬ 
tribution  of  dislocation  creep  to  the  cavity  growth  rate  is  shown  to  be  negligible  when  L  is  comparable  to  or  larger 
than  the  cavity  spacing,  but  significant  interactions  occur,  leading  to  growth  rates  very  much  in  excess  of  those 
predicted  on  the  basis  of  boundary  diffusion  alone,  when  L  is  comparable  to  or  smaller  than  the  cavity  size 

Needleman  and  Sham  have  begun  an  examination  of  the  effect  of  stress  triaxiality  on  cavity  growth  for  the  same 
axisymmetric  cavity  growth  mode;  work  thus  far  discussed  was  for  uniaxial  tension  Limited  results  thus  far  for 
nonlinear  viscous  materials  suggest  that  due  to  triaxiality,  the  range  where  the  interaction  between  creep  and  dif¬ 
fusion  is  important  changes  in  comparison  with  the  zero  hydrostatic  stress  case.  The  analysis  in  this  area  is  con¬ 
tinuing. 

The  work  cited  thus  far  has  assumed  that  surface  diffusion  is  rapid  enough  that  a  quasi-equilibrium  spherical-cap 
cavity  shape  is  retained.  However,  the  analysis  of  the  interaction  between  plastic  creep  and  diffusion  predicts  a 
growth  rate  so  rapid,  at  high  stress  levels,  that  this  condition  will  not  always  be  met.  Consequently.  Rice  and 
Sham  are  studying  the  effect  of  non-equilibrium  cavity  shapes  on  cavity  growth  for  the  combined  processes  of 
creep  and  diffusion.  Work  has  been  directed  toward  determining  the  cavity  profile  in  the  limiting  case  of  rapid 
growth  in  the  crack-like  mode  with  deformable  adjoining  grains.  A  preliminary  study  has  been  made  on  the  cavity 
profile  in  the  crack-like  mode,  when  the  adjoining  grains  are  moving  apart  rigidly  with  a  velocity  of  S/2.  Here  S  is 
the  effective  thickening  at  the  grain  boundary.  It  is  found  that  the  cavity  profile  is  affected  only  through  the  combi¬ 
nation  8  2v  where  v  is  the  cavity  growth  speed.  The  conditions  under  which  this  contribution  is  important  are  being 
examined. 

In  order  to  understand  elastic  relaxation  effects,  at  low  stress  conditions  for  which  there  is  negligible  plastic  flow,  a 
coupled  elasticity-diffusion  problem  has  been  formulated  by  Rice  and  Rubinstein.  The  response  of  partially  dam¬ 
aged  materials  (where  damage  is  represented  as  crack-like  cavities  on  grain  boundaries)  following  load  altera¬ 
tions  is  analyzed.  Again,  there  is  material  redistribution  from  the  void  surfaces  to  the  grain  boundaries,  and  in  this 
case  non-uniform  matter  accommodation  along  boundaries  is  accomplished  by  elastic  deformation. 

These  coupled  phenomena  are  modelled  in  terms  of  a  boundary  value  problem  for  the  elastic  stress  field,  which  is 
constrained  to  satisfy  the  grain  boundary  diffusion  and  transport  conservation  equations  on  the  grain  boundary, 
and  to  obey  stress  free  conditions  at  the  crack-like  cavity  surfaces.  On  the  basis  of  the  solution  for  the  stress  field 
produced  by  a  dislocation  suddenly  introduced  on  the  grain  boundary,  integral  equations  have  been  derived  for 
both  single  crack  and  periodically  distributed  cracks  on  the  grain  boundary.  The  integral  equations  have  been 
solved  numerically  by  Rubinstein  for  the  case  of  cyclic  loading,  and  Fourier  superposed  solutions  have  been 
developed  for  sudden  step  loadings. 

For  the  suddenly  applied  loading,  the  results  demonstrate  that  the  elastic  stress  concentration  is  effectively 
relaxed  by  diffusion  after  time  t  =  0.05t,  where  T=f3/DE,  <  is  half  of  the  distance  between  adjacent  tips  of 
neighboring  cracks  along  the  grain  boundary,  E  is  Young's  modulus,  and  D  is  the  diffusion  parameter  relating 
volumetric  flux  along  the  g.b.  to  stress  gradient.  By  t  =  0.25t  the  stress  distribution  becomes  essentially  identical 
to  that  for  rigid  grains,  i.e..  diffusion  has  become  the  dominant  process. 

Principal  Investigators:  A.  Needleman  (Engineering)  and  J.  R.  Rice  (Engineering). 

Personnel:  A.  Rubinstein  (Engineering)  and  T-L.  Sham  (Engineering). 

Publications:  “Plastic  Creep  Flow  Effects  in  the  Diffusive  Cavitation  of  Grain  Boundaries,"  A.  Needleman  and 
J.  R.  Rice,  Acta  Metallurgies,  Overview  Paper  No.  9,  28  (1980),  1315-32. 

“Energy  Variations  in  Diffusive  Cavity  Growth,”  J.  R.  Rice  and  T-j.  Chuang,  Journal  of  the  American  Ceramic 
Society,  in  press  for  1981. 

“Suppression  of  Cavity  Formation  in  Ceramics:  Prospects  for  Superplasticity,"  A.  G,  Evans,  J.  R.  Rice,  and 
J.  P.  Hirth,  Journal  of  the  American  Ceramic  Society  63  (1980),  368-75. 


Supported  by  the  U.S.  Department  of  Energy  and  the  Advanced  Research  Projects  Agency. 
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Fatigue  of  Aluminum  Alloys 

This  research  is  concerned  with  a  study  of  the  macro-  and  micro-mechanics  ol  fatigue  crack  growth  in  an  age- 
hardenable  aluminum  alloy,  e  g  .  2048  aluminum.  Special  attention  is  given  to  the  importance  of  microstructure, 
including  grain  size  and  aging  treatment  in  crack  growth,  on  the  phenomena  of  crack  tip  closure,  and  on  overload 
retardation  effects  Crack  closure  is  being  studied  in  detail  throughout  the  course  of  this  research  program  and  its 
use  in  describing  crack  growth  explored. 

This  research  has  led  to  the  further  refinement  of  experimental  techniques  used  to  measure  in-situ  crack  growth. 
This,  in  turn,  has  led  to  the  development  of  some  new  apparatus  for  measuring  small  changes  in  specimen  com¬ 
pliance  and  thus  small  increments  of  crack  growth. 

The  present  program  has  also  allowed  us  to  build  a  computer  interface  with  our  closed  loop  testing  equipment. 
This  is  now  working  as  a  data  acquisition  system. 

Principal  Investigator:  R.  J.  Asaro  (Engineering). 

Personnel:  L.  Hermann  (Engineering)  and  J.  M  Baik  (Engineering). 

Publications:  “Transitions  in  Fatigue  Crack  Closure  in  2048  Aluminum,"  R.  J.  Asaro,  L.  Hermann,  and  J.  M.  Baik, 
to  appear  in  Met.  Trans. 

“Fatigue  Crack  Growth  in  2048  Aluminum,"  J.  M.  Baik,  L.  Hermann,  and  R.  J.  Asaro,  in  preparation. 

Supported  by  the  Materials  Research  Laboratory  NSF  and  the  Air  Force  Office  of  Scientific  Research. 


Strength  and  Fracture  in  High  Strength  Titanium  Alloys 

Attention  during  the  past  year  focused  on  thermomechanical  processing  ol  Cobalt  Titanium  Alloys.  Alloys  of 
CO-5,  6,  7,  10  and  18%  Ti  were  prepared. 

These  consist  of  a  two  phase  mixture  of  HCP  CO  with  a  very  low  stacking  fault  energy  and  y,  0O3Ti,  which  is  an 
ordering  fee  alloy  of  the  Cu3Au  type.  After  annealing  at  1 100°C  and  quenching,  the  alloys  were  ductile  and  could 
be  cold  rolled  in  excess  of  50%.  Strength  determined  by  hardness  increased  rapidly  and  tended  to  a  maximum  for 
the  10%  Ti  alloy  equal  volume  fractions  of  both  phases.  Hardness  reached  Rc  36  at  8%  cold  deformation  and  Rc 
40-42  at  50%.  With  low  temperature  aging  treatments  hardness  increased  toTc  48  to  50.  From  previous  experi¬ 
ence  with  (J  titanium  alloys  we  expect  the  aging  treatment  to  become  increasingly  effective  for  heavy  deforma¬ 
tions  between  50  and  95%.  Larger  ingots  and  hot  forging  will  provide  material  for  cold  swagging  and  wiredrawing 
to  very  high  strain  levels. 

Principal  Investigator:  0.  H.  Avery  (Engineering). 

Supported  by  the  Materials  Research  LaboratoryfNSF. 
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SECTION  3 


Inorganic  Glasses 


Introduction 


The  research  ol  the  Inorganic  Glasses  Group  is  focused  on  understanding  the  structures,  relaxation  phenom¬ 
ena,  and  electronic  states  and  properties  of  the  three  principal  types  of  inorganic  glasses;  amorphous  metallic 
glasses,  amorphous  semiconductors,  and  oxide  glasses.  These  reports  of  this  research  in  1979-80  can  be  read 
using  either  the  phenomena  or  the  types  of  glass  as  the  organizing  feature.  But,  the  central  fact  that  emerges  is 
that  fundamental  issues  concerning  the  amorphous  state  in  general  are  being  addressed  by  studying  the  mater¬ 
ials. 

Both  the  properties  of  important  materials  and  discoveries  of  phenomena  are  reported.  Thus,  for  example,  the 
existence  of  the  magneto-optic  effect  in  metallic  glasses  was  demonstrated,  ultrasonic  attenuation  and  velocity 
changes  of  metallic  glasses  were  found  to  differ  by  orders  of  magnitude  from  the  predictions  of  the  two-level- 
system  tunneling  theory  (which  is  sufficient  for  insulating  glasses),  the  existence  of  a  mobility  edge  in  disordered 
electronic  systems  has  been  addressed  using  a  field  theoretic  model  of  interacting  matrices,  and  it  has  been 
shown  using  instanton  methods  that  a.c.  conductivity  from  the  localized  states  in  the  band  tail  of  a  disordered 
solid  is  given  as  o-(u>)~a)2(log  m)d  ’ 1  (d  is  the  space  dimensionality)  confirming  Anderson,  Halperin  and  Mott's 
predictions. 

Moreover,  the  existence  of  non-linear  composition  dependence  (pronounced  minimum)  of  isothermal  compres¬ 
sibility  (KJ0(T))  in  mixed  alkali  pentasilicates  was  shown,  strong  spectroscopic  evidence  for  the  existence  of 
intermediate  range  structures  in  B203  was  educed  from  isotopic  substitution  studies,  it  was  shown  that  the 
proposed  pure  metaborate  compound  does  not  exist  in  the  Mg0-B203  system,  the  distribution  of  transition 
metal  ions  in  oxide  glasses  was  studied  by  neighboring  atom  nuclear  spin  relaxation,  the  role  of  a  third  compo¬ 
nent  in  enabling  glass-formation  in  a  ternary  system  whose  binary  components  crystallize  was  studied,  and 
measurements  of  dielectric  loss  at  VHF  to  gigahertz  frequencies  on  mixed  alkali  silicate  glasses  by  time-domain 
methods  showed  that  they  exhibit  previously  unknown  composition-dependent  relaxations  in  this  range. 

In  addition,  the  relaxation  of  photoinduced  earners  in  a-Si:H  was  studied  and  shown  to  be  much  slower  than  in 
a-Si  or  chalcogenide  glasses  and  to  be  due  to  diffusion— limited  bimolecular  recombination.  The  dependence 
of  the  photoinduced  absorption  constant  for  silicon-rich  alloys  (Si1tGe J  and  a-Si:H  was  determined.  It  was 
interpreted  on  the  basis  of  a  model  in  which  the  photoinduced  absorption  band  is  due  to  transitions  of  holes 
from  trap  states  close  to  the  valence  band  into  the  valence  band  itself— a  model  which  should  permit  using  the 
picosecond  absorption  band  to  obtain  information  about  hole  traps  in  a-Si:H  and  related  materials. 
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SECTION  3 


Inorganic  Glasses 


Individual  Contributions 


Optical  and  Magneto-Optical  Properties  of  Metallic  Glass  Fe^^ 

We  measured  the  optical  constants  of  ferromagnetic  Fe^B^  glass  and  the  magneto-optic  Kerr  effect  (MOKE). 
The  optical  absorption  in  the  glass  is  compared  to  that  of  Fe  and  a  shift  of  the  structures  in  the  Fe^B^  spectrum 
towards  lower  energies  is  noted.  The  MOKE  absorption  spectrum  is  similar  in  shape  but  reduced  hy  a  factor  of 
about  five  compared  to  Fe.  This  exploratory  work  was  done  to  find  out  whether  the  magneto-optic  effect  exists  in 
metallic  glasses.  We  found  that  it  does  exist  and  is  measurable;  however,  the  interpretation  of  the  spectrum  was 
difficult  because  of  the  complexity  of  the  band  structure  of  Fe.  It  was  proposed  that  studies  of  this  effect  in  amorp¬ 
hous  materials  with  simpler  electronic  structures  may  be  more  informative  (for  example,  on  materials  containing 
Ni  instead  of  Fe;  unfortunately  these  materials  are  much  poorer  glass  formers). 

Principal  Investigator:  J.  Tauc  (Engineering  and  Physics). 

Personnel:  L.  Bieman  (Engineering)  and  S.  Ray  (Physics). 

Publication:  “Optical  and  Magneto-Optical  Properties  of  Metallic  Glass  Fe  B  ,'  S.  Ray  and  J.  Tauc,  Solid  State 
Comm.  34  (1980),  769. 

Supported  by  the  Materials  Research  Laboratory  I NSF  and  the  National  Science  Foundation. 


Ultrasonic  Studies  of  Metallic  Glasses 

The  results  of  thermal  and  ultrasonic  studies  of  insulating  glasses  at  low  temperatures  conducted  over  the  past 
few  years  are  generally  found  to  be  consistent  with  the  predictions  of  the  two-level-system  (TLS)  tunneling  theory. 
The  behavior  of  metallic  glasses  at  low  temperatures  is  less  clear.  Therefore,  measurements  of  ultrasonic  attenu¬ 
ation  and  velocity  changes  were  carried  out  on  the  metallic  glass  Pd0  ^Si,,  185Cuooe  as  a  function  of  ampli¬ 
tude,  in  the  frequency  range  10  to  90  MHz,  for  0.3  s  T  s  10  K.  A  comparison  of  the  experimental  results  with  the 
predictions  of  the  current  two-level-system  (TLS)  tunneling  theory  Is  presented.  The  amplitude-dependent  atten¬ 
uation  changes  observed  in  these  studies  are  larger  by  a  factor  of  100  to  1000  than  the  values  obtained  from  cur¬ 
rent  TLS  theory,  with  the  use  of  parameters  determined  experimentally  at  higher  frequencies  by  other  workers. 
Below  -  1  K  these  attenuation  changes  have  a  linear  (rather  than  quadratic)  dependence  on  frequency  and  very 
weak  temperature  dependence  (rather  than  1/T).  The  saturated  (amplitude-independent)  attenuation  was  found 
to  have  a  very  weak  temperature  dependence,  rather  than  T3  dependence.  A  significant  amplitude  dependence 
of  the  slope  of  the  sound  velocity  change  as  a  function  of  InT  was  also  observed. 

Principal  Investigator:  C.  Elbaum  (Physics  and  Applied  Mathematics). 

Personnel:  A.  Hikata  (Applied  Mathematics),  H.  Araki  (Materials  Research  Laboratory),  and  G.  Park  (Physics). 

Publications:  “Ultrasonic  Studies  of  Amorphous  PdSICu,"  H.  Araki,  G.  Park,  A.  Hikata,  and  C.  Elbaum,  Solid  Stale 
Comm.  32  (1979),  625. 

“Ultrasonic  Studies  of  Metallic  Glasses,*  H.  Arald,  G.  Park,  A.  Hikata,  and  C.  Elbaum,  Proceedings  of  the  Interna¬ 
tional  Conference  on  Phonon  Scattering  In  Condensed  Matter,  ed.  H.  J.  Marls  (Plenum,  1980),  p.  69. 

"Ultrasonic  Attenuation  and  Velocity  Studies  of  Amorphous  PdSiCu  at  Low  Temperatures,"  H.  Araki,  G.  Park, 
A.  Hikata,  and  C.  Elbaum,  Phys.  Rev.  B21  (1960),  4470. 
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Instantons  and  Disordered  Systems 

We  have  continued  our  work  on  disordered  systems.  It  has  been  shown  using  instanton  methods  (publication  1) 
that  a.c.  conductivity  from  the  localized  states  in  the  band  tail  of  a  disordered  solid  is  given  as  (r(u>)~u2(log 
u>)a  * 1 ,  d  is  the  space  dimensionality,  confirming  the  predictions  of  Anderson,  Halperin  and  Mott.  The  existence  of 
a  mobility  edge  in  disordered  electronic  systems  has  been  discussed  using  a  field  theoretic  model  of  interacting 
matrices  (publications  5  and  10).  It  is  shown  rigorously  that  all  states  are  localized  in  two  dimensions,  above  two 
dimensions  for  weak  disorder  there  are  mobility  edges,  but  these  merge  above  a  critical  amount  of  disorder  and 
all  states  become  localized. 

We  have  also  studied  topological  order  in  two  dimensions.  A  field  theory  of  dislocation  mediated  melting  in  two 
dimensions  has  been  developed,  a  coupled  vector  sine  Gordon  theory.  Thus  we  may  now  study  this  transition 
systematically  (publication  8).  The  effect  of  a  random  symmetry  breaking  field  on  topological  disorder  in  two 
dimensional  systems  has  been  studied  Such  a  field  would  simulate  the  interaction  of  the  idealized  two  dimen¬ 
sional  system  with  an  underlying  substrate  if  the  substrate  were  to  contain  patches  (islands)  in  which  short  range 
order  is  retained  rather  than  being  a  single  crystalline  surface,  e  g.,  charge  density  waves  in  chemisorption  and 
physisorption.  Planar  model  behavior  should  be  observable  for  hexagonal  perturbations  (publication  7).  In 
related  work  we  have  shown  that  the  ratio  of  any  two  correction  to  scaling  amplitudes  is  universal  and  have  calcu¬ 
lated  the  ratios  to  second  order  in  *  =  4-d.  Agreement  with  experiment  when  available  is  excellent  (publications  2, 
3.  and  9). 

In  the  area  of  surface  magnetism,  the  magnetization  shape  function  at  a  free  surface  has  been  calculated  explic¬ 
itly  using  renormalized  perturbation  theory.  To  leading  order 

M(z)  =  t@(6/g),2c<t>(y) 


<t>(y)  =  y°,  o  =  1  - 1! 6 


Here  3  is  the  bulk  magnetization  critical  exponent  (publication  6).  We  have  also  derived  explicitly,  using  instanton 
methods,  the  form  of  the  essential  singularity  at  the  coexistence  curve  in  an  Ising  ferromagnet  (publication  4). 

Principal  Investigator:  A.  Houghton  (Physics). 

Personnel:  M-C.  Chang  (Physics) 

Publications:  “Two  Particle  Spectral  Function  and  a.c.  Conductivity  of  an  Amorphous  System  Far  Below  the 
Mobility  Edge.  A  Calculation  of  Interacting  Instantons,"  A.  Houghton,  L.  Schafer,  and  F.  J  Wegner,  Phys.  Rev.  B 
22  (1980),  3598. 

“Critical  Amplitude  Ratio  of  the  Confluent  Singular  Term  for  the  Specific  Heat:  Calculations  to  Order  *2  for  Sys¬ 
tems  of  Continuous  Symmetry,”  A.  Houghton  and  M-C.  Change,  Phys.  Rev.  .8  21  (1980),  1881. 

“Universal  Ratios  Among  Correction  to  Scaling  Amplitudes  on  the  Coexistence  Curve,”  A.  Houghton  and 
M-C.  Chang,  Phys.  Rev.  Lett.  44  (1980),  785. 

The  Metastable  Ising  Magnet  in  a  Negative  Field,"  A.  Houghton  and  T.  C.  Lubensky,  Phys.  Lett.  77 A,  (1980),  479. 

“Noncompact  a  Models  and  the  Existence  of  a  Mobility  Edge  in  Disordered  Electronic  Systems  Near  Two  Dimen¬ 
sions,"  A.  Houghton,  A.  Jevicki,  R.  D.  Kenway,  and  A.  M.  Pruisken,  Phys.  Rev.  Lett.  45,  (1980),  394. 

“Singularities  of  the  Surface  Magnetisation  Profile  near  the  Critical  Point,”  A.  Houghton,  D.  R.  Grempel,  and 
S.  C.  Ying,  Phys.  Lett.  A  78,  (1980),  295. 

The  Effects  of  a  Random  Symmetry  Breaking  Field  on  Topological  Order  in  Two  Dimensions,"  A.  Houghton, 
R.  D.  Kenway,  and  S.  C.  Ying,  Phys.  Rev.  B  23,  (1981),  289. 
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“A  Field  Theory  of  Dislocation  Mediated  Melting  in  Two  Dimensions.''  A  Houghton  and  M.  C  Ogilvte.  J.  Phys  A 
L499-L454  (1980) 

“Correction  to  Scaling  Amplitude  Ratios  for  the  Superfluid  Transition  in  ''He."  A  Houghton  and  M-C  Chang. 
Phys  Rev  B  to  be  published. 

“Noncompact  <r  Models  and  the  Existence  of  a  Mobility  Edge  in  Disordered  Electronic  Systems  near  Two  Dimen¬ 
sions."  A  Houghton,  A.  Jevicki.  R  D.  Kenway.  and  A  M  Pruisken.  Proceedings  of  the  International  Conference 
on  Statistical  Mechanics,  Edmonton.  Alberta,  Canada,  1980 

Supported  by  the  Materials  Research  Laboratory  NSF  and  the  National  Science  Foundation. 


Time  Resolved  Spectroscopy  of  Amorphous  Semiconductors 

We  continued  to  work  on  the  relaxation  of  photoinduced  absorption  in  amorphous  semiconductors.  In  the  previ¬ 
ous  study  we  found  that  in  non-hydrogenated  amorphous  silicon  and  chalcogenide  glasses  the  photogenerated 
carriers  decay  very  fast  (on  a  picosecond  scale)  but  in  hydrogenated  amorphous  silicon  (a-Si:H)  the  decay  was 
much  longer  (exceeding  the  maximum  time  of  150  psec  measurable  with  our  picosecond  equipment).  Recently 
we  studied  the  decay  in  a-Si:H  in  the  time-range  of  0.5  usee  to  10  msec  and  found  that  the  decay  can  be 
described  as  due  to  diffusion  limited  bimolecular  recombination.  The  diffusion  constant  depends  on  time  as 
1 11";  the  parameter  a  (<1 )  was  found  to  be  between  0.6  to  0.7  in  sputtered  a-Si:H  while  in  glow-discharge  a-Si:H 
it  vaned  with  temperature.  With  our  model  for  the  photoinduced  absorption  band  obtained  from  steady  state  stud¬ 
ies.  the  decay  data  can  be  used  to  get  information  about  the  transport  of  electrons  (more  mobile  carriers)  which 
recombine  with  trapped  holes  The  time  dependence  of  the  diffusion  constant  points  out  the  dispersive  character 
of  the  transport,  as  observed  also  in  other  experiments  on  these  materials,  and  is  interpreted  on  the  basis  of  the 
continuous-random-walk  theory 

Principal  Investigator  J  Tauc  (Engineering  and  Physics). 

Personnel:  Z  Vardeny  (Engineering)  and  S.  Ray  (Physics). 

Publications:  “Picosecond  Relaxations  in  Amorphous  Semiconductors,"  D.  E.  Ackley.  J.  Tauc,  and  W.  Paul,  J. 
Noncrystalline  Solids  35  &  36,  (1980),  957. 

"Optical  Studies  of  Excess  Carrier  Recombination  in  a-Si:H:  Evidence  for  Dispersive  Diffusion,"  Z.  Vardeny, 
P  O'Connor,  S  Ray,  and  J  Tauc,  Phys.  Rev.  Lett.  44,  (1980),  1267 
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Photoinduced  Optical  Absorption  in  Amorphous  Hydrogenated  Silicon 

We  observed  an  infrared  absorption  band  below  the  absorption  edge  induced  by  illumination  in  the  fundamental 
absorption  band  in  hydrogenated  amorphous  tetrahedrally  bonded  semiconductors,  as  reported  previously.  In 
a-Ge:H  it  was  plausible  to  ascribe  this  band  to  photoinduced  polaron  hopping  but  in  a-Si:H  the  spectrum  could  not 
be  interpreted  in  this  way.  Recently  we  found  that  in  a-Si:H  and  silicon  rich  a-Si,  xGeK  alloys  the  photoinduced 
absorption  constant  Aa~(fiu>  -  E0)'  2/fiu>  suggesting  transitions  from  a  state  with  a  relatively  sharply  defined 
energy  to  a  band  with  the  square  root  density  of  electronic  states.  We  proposed  a  model  in  which  the  photoin¬ 
duced  absorption  band  is  ascribed  to  transitions  of  holes  from  trap  states  dose  to  the  valence  band  into  the 
valence  band.  Using  this  model,  study  of  the  picosecond  absorption  band  can  be  used  to  obtain  information  about 
the  hole  traps  in  a-Si.H  and  related  materials. 

Principal  Investigator:  J.  Tauc  (Engineering  and  Physics). 

Personnel:  P.  O’Connor  (Physics). 
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Publications:  “Photoinduced  Absorption  in  Hydrogenated  Amorphous  Tetrahedrally-Bonded  Semiconductors,” 
P.  O'Connor  and  J.  Tauc,  J.  Noncrystalline  Solids  35  &  36,  (1980),  699, 

"Spectrum  of  Photoinduced  Optical  Absorption  in  a-Si:H,”  P,  O'Connor  and  J.  Tauc,  Solid  State  Comm.  36, 
(1980),  947 

“Investigation  of  Defect  Electronic  States  and  Recombination  in  Amorphous  Tetrahedrally  Bonded  Semiconduc¬ 
tors  by  Photoinduced  Absorption,"  P.  O'Connor  and  J.  Tauc,  to  be  published. 
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Composition  Dependence  of  Isothermal  Compressibility  in  Mixed  Alkali  Glass  Form¬ 
ing  Systems 

In  connection  with  our  experimental  and  theoretical  efforts  to  understand  the  mixed  alkali  effects  on  a  range  of 
properties,  especially  as  they  appear  in  dynamic  mechanical  properties,  we  have  sought  to  relate  the  composition 
dependence  of  the  properties  of  the  glass-forming  systems  at  their  state  at  the  fictive  temperature,  Tf,  to  those  of 
the  glasses  themselves.  Glasses  in  the  series  xCs20(1-x)Na205Si02  were  studied  and  a  novel  and  significant 
result  was  obtained  in  the  variation  of  the  low  frequency  bulk  modulus  at  T,  with  composition,  i.e.,  Kto(T,)(x).  It 
was  found  that  Kt0(T()  has  a  pronounced  minimum  with  respect  to  x.  This  conclusion  was  made  possible 
because  our  measurements  also  show  that  n2  varies  monotonically  with  x,  so  that  (de/dC)PT  changes  slowly 
with  x.  and  that  T  vanes  slowly  with  x.  so  the  isothermal  compressibility  can  be  calculated  assuming  that  density 
variations  primarily  affect  light  scattering  intensity  ratios.  The  fact  that  KT  0(Tf)  varies  non-linearly  with  x  allows 
us  to  pose  the  question  of  whether  this  property  of  the  materials  at  T,  is  relevant  to  those  of  the  glasses  which, 
although  they  are  in  a  different  state  below  T(,  are  the  material  actually  measured  and  must  retain  the  physical 
evidence  of  the  T(  properties  in  some  form.  The  hypothesis  is  that  they  do  in  the  form  of  spatial  variations  in  strain 
state,  which  manifest  themselves  in  the  response  of  the  glasses  to  mechanical  stresses  in  dynamic  mechanical 
experiments  and  the  mixed  alkali  effect  on  low  frequency  mechanical  loss. 

Principal  Investigator:  W.  M.  Risen.  Jr.  (Chemistry). 

Personnel:  E.  I.  Kamitsos  (Chemistry). 

Publication:  "Mixed  Alkali  Glasses:  xCs20(1-x)Na20'5Si02,“  W.  M.  Risen,  Jr„  G.  B.  Rouse,  and  E.  I.  Kamitsos, 
J.  Noncryst.  Solids,  1981,  in  press. 
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Isotopic  Study  of  Structure  in  Borate  Glasses 

The  basic  question  about  the  structure  of  oxide  glasses  concerns  the  randomness  of  the  systems.  This  question 
has  been  joined  at  the  level  of  intermediate  range  order,  since  short  range  order  follows  chemical  bonding  princi¬ 
ples  in  them  and  they  are  disordered  over  long  ranges  in  the  absence  of  orienting  stresses.  While  special  inter¬ 
mediate  range  structures  have  been  identified  spectroscopically  in  our  work  on  arsenatophosphates,  the  continu¬ 
ous  random  network  (CRN)  approach  holds  that  the  spectra  and  properties  of  bulk-quenched  glasses  are 
consistent  with  random  connections  and  orientations  of  all  possible  types.  This  issue  is  central  to  the  structure  of 
B203(gl),  since  the  CRN  approach  has  been  proposed  as  the  best  way  to  treat  the  X-ray  distribution  function 
while  all  previous  spectral  work  had  led  to  the  proposition  that  intermediate  range  structures — boroxyl  rings — are 
the  dominant  features. 

In  this  work  the  B203  glass  system  was  prepared  in  the  isotopic  forms  Pure  B20’*  (gl),  pure  B20’3  (gl),  and, 
by  several  routes,  0.5  B20’*  0.5  B20’®  (gl),  and  their  Raman  spectra  obtained  and  analyzed.  Since  the 
principal  features  in  such  spectra  arise  from  symmetric  modes  based  on  B-0  stretches,  it  is  possible  to  distin¬ 
guish  between  a  boroxyl  ring-based  structure  and  a  CRN  structure.  The  CRN  predicts  that  the  principal  mode 
arises  from  the  density  of  states  as  a  collective  mode  which  would  shift  as  one  band  with  a  frequency  given  by  the 
average  oxygen  mass  in  the  network.  The  boroxyl  ring-based  structure  would  lead  to  four  bands,  in  the  intensity 
ratio  1:33:1,  for  the  50%  B20’®  50%  B20'®  glass,  with  the  bands  giving  rise  to  the  relative  intensity  of  1 
being  at  the  positions  seen  for  the  pure  isotopic  glasses.  That  was  observed.  The  glass  is  composed  nearly  com¬ 
pletely  of  boroxyl  rings,  and  any  CRN  character  must  apply  only  to  the  connections  between  these  large  struc¬ 
tures. 
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average  oxygen  mass  in  ttie  network.  The  boroxyl  ring-based  structure  would  lead  to  tour  bands,  in  the  intensity 
ratio  1:33:1,  for  the  50%  B20'®  50%  B20'3  glass,  with  the  bands  giving  rise  to  the  relative  intensity  of  1 
being  at  the  positions  seen  for  the  pure  isotopic  glasses.  That  was  observed.  The  glass  is  composed  nearly  com¬ 
pletely  of  boroxyl  rings,  and  any  CRN  character  must  apply  only  to  the  connections  between  these  large  struc¬ 
tures. 

Principal  Investigator:  W.  M.  Risen,  Jr.  (Chemistry). 

Personnel:  C.  H.  Windisch  (Chemistry). 

Publications:  “Raman  Study  of  Oxygen-labelled  B203  Glasses,"  W.  M.  Risen,  Jr.  and  C.  H  Windisch,  J.  Noncrys¬ 
talline  Solids.  1961,  in  press. 

“Spectra  and  Structure  of  Isotopically  Labelled  Alkali  Borate  Glasses,'  W.  M.  Risen,  Jr.  and  C.  H  Windisch,  J. 
Noncrystalline  Solids,  to  be  submitted. 
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Nuclear  Magnetic  Resonance  Studies  of  iron  Lead  Borate  Glasses 

Oxide  glasses  containing  greater  than  3  Mol  %  Fe203  are  generally  superparamagnetic  or  ferromagnetic  making 
the  observation  of  nuclear  magnetic  resonance  (NMR)  for  the  nuclei  of  non-magnetic  atoms  in  these  glasses 
impossible.  However,  glasses  containing  as  much  as  40  Mol  %  Fe203  which  remain  paramagnetic  can  be  pro¬ 
duced  in  the  system  xFe203  -  (1-x)[PbO  SBjOjj.  Hence,  this  system  is  ideal  for  studying  B11  and  Pb207  NMR  in 
an  iron-rich  glass. 

The  linewidth  of  the  narrow  response  of  the  B"  NMR  spectrum  has  been  obtained  as  a  function  of  composition 
and  applied  field  These  data  are  in  good  agreement  with  a  model  in  which  the  iron-ions  are  randomly  distributed 
in  the  glass  matrix.  In  addition,  computer  simulations  of  the  B"  and  Pb207  NMR  iineshape  have  been  compared 
with  the  experimental  spectra  to  derive  values  of  B"  quadrupole  interaction  parameters,  ratios  of  4-coordinated 
to  3-coordinated  atoms,  and  Pb207  chemical  shift  interaction  parameters  as  a  function  of  composition. 

Principal  Investigator:  P.  J.  Bray  (Physics). 

Personnel:  See  below  (Organic  Solids). 

Publications:  See  below  (Organic  Solids) 

Supported  by  the  Materials  Research  Laboratory/NSF  and  the  National  Science  Foundation. 


O17  and  H1  NMR  Studias  of  Si02 

SiO,  enriched  to  50%  O’7  has  been  fabricated  in  conjunction  with  Professor  William  Risen  and  his  students  in 
the  Department  of  Chemistry.  O'7  and  H1  NMR  spectra  of  the  resulting  amorphous  powder  have  been  obtained. 
The  O17  spectrum  exhibits  sharply  defined  structure.  The  H'  spectrum  indicates  protons  are  present  in  a  concen¬ 
tration  of  approximately  102  ppm,  and  they  are  not  present  as  Interstitial  water  molecules  but  probably  exist  as 
hydroxyl  groups  in  the  glass. 

Principal  Investigator:  P.  J.  Bray  (Physics). 

Personnel:  See  below  (Organic  Solids). 

Publications:  See  below  (Organic  Solids). 

Supported  by  the  Materials  Research  Laboratory/NSF  and  the  National  Science  Foundation. 
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NMR  Studies  of  Glasses  and  Polycrystaillne  Materials  in  the  ZnO  •  B203  System 

The  B”  NMR  spectra  lor  zinc  borate  compounds  were  secured  and  analyzed  to  identity  the  boron-oxygen  struc¬ 
tural  groupings  present  in  these  materials.  Results  for  polycrystalline  materials  at  the  3:1  and  4:3  compositions 
(zinc  oxide  to  boron  oxide)  were  found  to  be  indicative  of  boron-oxygen  groupings  which  correlate  with  structures 
proposed  in  the  literature.  A  third  compound  at  the  diborate  (1:2)  composition  is  currently  under  investigation. 

Zinc  borate  glasses  were  made  and  studied.  The  NMR  spectra  were  found  to  be  very  similar  throughout  the  range 
of  glass  formation.  This  indicates  that  a  similarity  in  structure  exists  between  all  the  glasses  in  this  system.  The 
fraction  of  boron  in  four-coordination.  N4,  was  determined  for  each  of  the  glasses.  A  consistent  difference  in  N4 
was  found  for  glasses  having  the  same  compositions  with  different  thermal  histories. 

Principal  Investigator:  P.  J.  Bray  (Physics). 

Personnel:  See  below  (Organic  Solids). 

Publications:  See  below  (Organic  Solids). 

Supported  by  the  Materials  Research  Laboratory/NSF  and  the  National  Science  Foundation. 


The  Structures  of  Compounds  and  Glasses  In  the  System  MgO  -  B203  Using  B11  NMR 

Structural  groupings  in  the  system  MgO  -  B203  were  studied  by  the  technique  of  computer  fitting  B"  NMR 
experimental  lineshapes.  The  results  are  strongly  dependent  on  sample  preparation  and  indicate  that  the  pure 
metaborate  compound  does  not  exist.  Previous  models  of  the  structure  of  glasses  in  this  system  must  be  reexa¬ 
mined  in  light  of  these  findings. 

Principal  Investigator:  P  J.  Bray  (Physics). 

Personnel:  See  below  (Organic  Solids). 

Publications:  See  below  (Organic  Solids). 
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Tellurium  Borate  Glaaaes 

Tellurium  borate  glasses  containing  either  lithium  oxide  or  aluminum  oxide  as  third  components  have  been  made 
and  studied  using  NMR.  The  binary  tellurium  borate  system  was  found  to  have  a  large  region  of  immiscibility. 
However,  homogeneous  glasses  could  be  made  through  the  addition  of  lithium  oxide  or  alumina.  The  B1'  NMR 
spectra  for  glasses  in  the  above  systems  were  secured  and  used  to  study  the  behavior  of  N4,  the  fraction  of  boron 
atoms  in  four-coordination.  The  N4  results  indicate  that  the  Te02  in  the  glasses  contributes  to  the  formation  of 
four-coordinated  boron  atoms.  There  is  an  apparent  change  in  the  coordination  of  tellurium  atoms  from  four  to 
three  with  the  subsequent  release  of  oxygen  atoms.  The  Te1*9  spectra  in  these  glasses  are  under  investigation 
and  will  aid  in  understanding  both  the  behavior  of  tellurium  and  boron  in  these  materials. 

Principal  Investigator:  P.  J,  Bray  (Physics). 

Personnel:  See  below  (Organic  Solids). 

Publications:  See  below  (Organic  Solids). 

Supported  by  the  Materials  Research  Laboratory/NSF  and  the  National  Science  Foundation. 
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Time  Domain  Spectroscopy  and  Dielectric  Properties  of  Glassy  Materials 

Preliminary  measurements  of  several  mixed  alkali  glasses  made  by  W.  M.  Risen's  research  group  have  indicated 
small  but  definite  losses  at  VHP  to  gigahertz  frequencies  which  are  composition  dependent  and  in  excess  of  ordi¬ 
nary  conduction  effects  Work  in  progress  is  directed  at  better  definition  of  these  losses  after  further  refinement  of 
time  domain  difference  methods  previously  developed  in  this  laboratory. 

A  program  in  collaboration  with  G.  T  Evans  (Oregon  State)  was  started  to  test  his  Brownian  dynamics  calcula¬ 
tions  of  kinetics  of  polymer  chain  configurational  changes.  These  can  be  probed  by  dielectric  relaxation  of  electric 
dipole  orientations,  and  there  are  intriguing  indications  that  a  basis  for  molecular  understanding  of  the 
“skewed-arc"  relaxation  function  found  for  many  solid  materials  may  be  at  hand,  at  least  for  linear  chain  polymers. 

Related  work,  with  support  primarily  from  other  sources,  was  a  series  of  measurements  of  ion-dipole  interaction 
effects  in  a  series  of  electrolyte  solutions  as  a  further  test  of  Hubbard-Onsager  and  Wolynes-Hubbard  theories  of 
kinetic  depolarization.  This  until  recently  unrecognized  effect  was  confirmed  by  our  measurements  of  sulfuric 
acid,  and  is  now  believed  to  have  a  significant  bearing  on  ionic  mobilities  in  a  variety  of  systems. 

Considerable  time  was  spent  preparing  two  review  manuscripts  for  publication:  one  on  dipole  correlation  function 
theories  of  dielectric  relaxation  in  solids  (for  a  British  Institute  of  Physics  symposium),  and  the  other  on  Fourier 
transform  methods  for  study  of  dielectric  polarization  (as  part  of  an  ACS  symposium  volume). 

Principal  Investigator:  R.  H.  Cole  (Chemistry). 

Publication:  "Evaluation  of  Dielectric  Behavior  by  Time  Domain  Spectroscopy.  3.  Precision  Difference  Methods,” 
R.  H.  Cole,  S.  Mashimo,  P.  Winsor,  IV,  J.  Phys.  Chem.  84  (1980),  786. 
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Organic  Solids 

Both  N'4  and  Cl35  NQR  data  have  been  obtained  for  a  series  of  compounds  containing  nitrogen  directly  bonded 
to  chlorine.  The  existence  of  an  inverse  correlation  between  the  N’4  and  Cl35  quadrupole  coupling  constants  is 
interpreted  in  terms  of  a  simple  model  dealing  with  charge  excesses  and  deficiencies  at  the  respective  nuclei.  The 
relationship  between  the  ionic  character  of  the  N-CI  bond  and  the  NQR  data  is  discussed,  and  the  possibility  of 
obtaining  a  value  for  the  atomic  quadrupole  coupling  constant  of  N14  (e?qpQ)  is  noted.  The  Townes  and  Dailey 
method  of  computing  electron  charge  densities  from  the  NQR  frequencies  is  shown  to  be  satisfactory  in  giving 
consistent  results,  provided  that  the  hybridization  of  the  nitrogen  orbitals  is  known. 

The  ,4N  NQR  spectra  of  22  monosubstituted  anilines  have  been  obtained  and  analyzed  in  the  framework  of  the 
Townes  and  Dailey  theory.  Satisfactory  correlation  of  the  NQR  data  with  both  the  Hammett  a  parameters  and  the 
in  vitro  biological  activities  of  the  corresponding  sulfanilimides  have  been  found.  The  results  obtained  suggest 
that  the  nitrogen  lone-pair  orbital  is  more  sensitive  than  the  nitrogen-carbon  sigma  orbital  is  to  substituent  effects. 

,4N  NQR  spectra  have  been  obtained  at  77°K  for  several  organic  dye  cations  and  one  free  base.  The  cations  are 
weakly  paramagnetic  at  77°K  exhibiting  broad  and  featureless  ESR  spectra  characteristic  of  weak  charge-trans¬ 
fer  complexes.  The  positive  charge  in  the  cations  is  largely  delocalized.  The  existence  of  a  simple  relationship 
between  the  ,4N  NQR  data  and  optical  properties  for  three  cyanine  laser  dyes  has  been  studied. 

,4N  NQR  spectra  have  been  obtained  for  a  series  of  monomeric  carbazole  compounds  including  N-vinyl  carba- 
zole,  which  in  the  polymeric  form  of  this  compound  is  one  of  the  most  sensitive  photoconductlve  organic  poly¬ 
mers.  Electron  distributions  determined  from  the  NQR  data  show  clearly  and  quantitatively  the  delocalization  of  pi 
electrons  in  carbazole  as  compared  to  pyrrole.  Effects  of  substituents  such  as  aldehyde  and  vinyl  groups  on  the 
molecular  charge  distributions  have  been  studied. 

Principal  Investigator:  P.  J.  Bray  (Physics). 

Personnel:  F.  Buchoitz  (Physics),  W.  J.  Deil  (Physics),  A.  E.  Gelssburger  (Physics),  I.  A.  Harris  (Physics), 
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S.  G.  Greenbaum  (Physics),  M.  Lui  (Physics),  and  M.  L.  Buess  (Physics). 

Publications:  “B”  NMR  Studies  of  the  Structure  of  Borate  Glasses,"  P  J.  Bray,  S.  A.  Feller,  G.  E.  Jellison,  Jr., 
and  Y.  H.  Yun,  J.  Non-Crystalline  Solids,  38  &  39  (I960),  93-98. 

“Si2®  NMR  Spectra  of  Glasses  and  Polycrystalline  Compounds  in  the  System  K20-Si02,"  I.  A.  Harris,  Jr.  and 
P.  J.  Bray,  Phys.  and  Chem.  of  Glasses  21,  (1980),  156. 

“B"  Nuclear  Magnetic  Resonance  Studies  of  Li20-B203  Glasses  at  High  Li20  Content,”  Y.  H.  Yun  and 
P.  J.  Bray,  submitted  to  J.  Non-Crystalline  Solids. 

The  Use  of  Magnetic  Fields  for  ,4N  NQR  Spin-Echo  Signal  Enhancement,”  M.  L.  Buess  and  P.  J.  Bray,  accepted 
for  publication  in  the  J.  Mag.  Res. 

"N’4  Nuclear  Quadrupole  Resonance  in  Carbazoles,”  N.  S.  Kim  and  P.  J.  Bray,  submitted  to  J.  Organic  Magnetic 
Resonance 

“B"  NMR  Studies  of  Glasses  in  the  System  xFe203(1-x)[PbO-3B203],"  F.  Bucholtz  and  P.  J.  Bray,  Am.  Cer. 
Soc.  Bull.  59,  (1980),  863. 

“Investigation  of  Nitrogen-Chlorine  Bonds  by  14N  and  “Cl  NQR,”  S.  G.  Greenbaum  and  P.  J.  Bray,  Bull.  Am. 
Phys.  Soc.,  25(1980),  272. 

“,4N  Nuclear  Quadrupole  Resonance  Study  of  Substituted  Anilines  and  Prediction  of  in  vitro  Activity  of  Sulfani- 
lamides,"  S.  G.  Greenbaum,  S.  N.  Subbarao,  P.  J.  Bray,  and  T.  Oja,  J.  Organic  Magnetic  Resonance  14(1980), 
379. 

”'4N  Nuclear  Quadrupole  Resonance  in  Weakly  Paramagnetic  Organic  Dye  Cations,"  S.  G.  Greenbaum  and 
P.  J.  Bray,  submitted  to  J.  Magnetic  Resonance. 

“The  Determination  of  the  Structures  of  Compounds  and  Glasses  in  the  System  Mg0-B203  Using  B"  NMR,” 
M.  J.  Park,  K.  S.  Kim,  and  P.  J.  Bray,  Phys.  Chem.  Glasses  20,(1979),  31. 
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SECTION  4 

Chemisorption  and  Related  Surface  Interactions 


Introduction 


The  work  in  this  thrust  area  is  concerned  with  the  basic  properties  of  solid  surfaces  and  with  the  fundamental 
interactions  occurring  in  chemisorption  processes  and  surface  reactions.  The  experimental  studies  cover  a 
range  of  model  substrates,  from  well-defined  single-crystal  surfaces  to  metal  clusters  on  certain  polymers.  The 
theoretical  work  extends  from  a  simple  jellium  picture  to  more  realistic  models  of  alkali  and  transition  metals. 

Static  surface  properties,  such  as  crystallographic  and  electronic  structure,  have  been  investigated  experimen¬ 
tally  by  Estrup  and  Stiles,  using  low  energy  electron  diffraction  (LEED),  surface  reflectance  spectroscopy 
(SRS),  photoelectron  spectroscopy  (UPS,  XPS),  and  other  methods.  At  the  same  time,  Ying  has  carried  out  the¬ 
oretical  work  on  the  electronic  and  magnetic  properties  of  metal  surfaces,  and  has  considered  the  effects  of  sur¬ 
face  defects.  Several  studies  are  underway  concerning  phase  diagrams  and  surface  phase  transitions;  this 
work  is  likely  to  receive  increasing  emphasis  and  will  involve  a  particularly  close  collaboration  between  experi¬ 
mentalists  and  theorists. 

Studies  of  dynamic  surface  properties  include  investigations,  by  Greene,  of  the  scattering  and  ionization  of 
atomic  beams;  and  the  laser-spectroscopic  investigations,  by  Diebold,  of  collisional  energy  transfer  to  diatomic 
molecules.  In  addition,  Estrup  has  studied  the  kinetics  of  dissociation  and  desorption  of  simple  molecules  on 
metal  surfaces;  and  Risen  has  investigated  the  reaction  of  similar  species  on  metals  supported  in  ionomers. 
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SECTION  4 

Chemisorption  and  Related  Surface  Interactions 


Individual  Contributions 


Experimental  Studies  of  Chemisorption 

Studies  of  surface  st'jcture  and  surface  kinetics  have  been  carried  out  on  well-defined  systems  consisting  of 
chemisorbed  layers  on  single  crystals  of  molybdenum  and  tungsten. 

The  nature  of  the  reconstructed  Mo(001)  and  W(001)  surfaces  has  been  investigated  by  LEED  and  associated 
techniques.  On  the  clean  substrates,  at  low  temperature,  the  surface  atoms  are  found  to  be  displaced  along  the 
■  1 1  >  directions  but  adsorption  of  hydrogen  (and  probably  other  adsorbates)  causes  a  shift  to  the  ■10>  directions. 
Evidence  has  been  obtained  indicating  that  the  surface  phases  result  from  a  periodically  varying  adatom  occupa¬ 
tion  probability  coupled  to  the  displacement  waves.  Data  for  hydrogen  desorption  suggest  that  the  substrate  rear¬ 
rangement  is  the  major  cause  of  the  observed  coverage  dependence  of  the  adsorbate  properties.  Similarly,  stud¬ 
ies  by  photoelectron  spectroscopy  of  the  kinetics  of  CO  dissociation  on  Mo(001)  show  that  the  rate  behavior  can 
be  explained  by  an  activation  energy  which  increases  at  the  same  coverage  at  which  the  substrate  changes  struc¬ 
ture. 

In  collaboration  with  P.  J  Stiles  (see  below),  surface  reflectance  studies  have  been  made  of  chemisorbed  layers 
on  the  W(1 10)  surface.  Information  about  adsorbate  bonding  and  about  surface  phase  transition  has  been 
obtained  by  this  technique. 

Principal  Investigator;  P  J.  Estrup  (Physics  and  Chemistry). 

Personnel .  R.  A.  Barker  (Chemistry),  R.  Robertshaw  (Physics),  and  S.  Semancik  (Physics). 

Publications:  “Simple  DC  Method  for  the  Continuous  Measurement  of  Work  Function  Changes,”  P.  J.  Estrup, 
R.  A.  Barker,  T  E.  Felter,  and  S.  Semancik,  J.  Vac.  Sci.  Technol.  17  (1980),  755. 

"Surface-Reflectance-Spectroscopy  Studies  of  H  on  W(110):  Surface  Band  Structure  and  Adsorbate  Geometry.” 
P.  J.  Estrup,  G.  Blanchet,  and  P.  J.  Stiles,  Phys.  Rev.  Lett.  44  (1980),  171. 

“Atomic  Displacements  on  the  Mo(001)  and  H/Mo(001)  Surfaces,"  P.  J.  Estrup,  R.  A.  Barker,  and  S.  Semancik, 
Surlace  Sci.  94  (1980),  L162. 

“Photoemission  Method  for  the  Study  of  Surface  Kinetics:  CO  on  Mo(100),"  P.  J.  Estrup  and  S.  Semancik.  J.  Vac. 
Sci.  Technol.  17  (1980),  233. 

“Substrate  Rearrangement  and  2D  Ordering,"  P.  J.  Estrup,  in  Ordering  in  Two  Dimensions,  S.  Sinha,  ed. 
(North-Holland)  39  (1980). 

Supported  by  the  Materials  Research  LaboratorylNSF. 


Surface  Reflectance  Spectroscopy 

It  has  been  found  that  surface  reflectance  spectroscopy  (SRS)  is  a  powerful  tool  for  the  study  of  surface  electronic 
structure  and  bonding.  By  including  reflectance  measurements  at  different  angles  of  polarization,  we  have 
obtained  information  not  only  about  the  band  structure  but  also  about  the  geometry  of  the  adsorbate  complex. 
Recent  work  has  shown,  furthermore,  that  surface  phase  transitions  give  rise  to  characteristic  reflectance 
changes,  and  SRS  has  been  employed  to  map  the  two-dimensional  phase  diagram  (in  the  temperature-coverage 
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plane)  for  the  O'W(IIO)  chemisorption  system. 

Principal  Investigator  P.  J.  Stiles  (Physics). 

Personnel:  G.  Blanchet  (Physics)  and  B  Parker  (Physics). 

Publications:  “Chemisorption-Induced  Changes  of  the  Surface  Dielectric  Function  for  W(110),’’  P  J.  Stiles  and 
G  Blanchet,  Phys.  Rev.  B21  (1980),  3273. 

"Surface-Reflectance-Spectroscopy  Studies  of  H  on  W(110);  Surface  Band  Structure  and  Adsorbate  Geometry," 
P  J  Stiles,  G.  Blanchet,  and  P  J.  Estrup,  Phys.  Rev.  Lett.  44  (1980),  171. 

Supported  by  the  Materials  Research  Laboratory INSF. 


Phase  Transitions  in  Chemisorbed  Systems 

We  have  developed  a  statistical  mechanical  model  for  the  structural  phase  transitions  observed  on  many  metal 
surfaces.  The  effect  of  adsorbates  on  this  transition  is  studied  in  detail.  Depending  on  the  coverage  and  tempera¬ 
ture.  the  adsorbate  can  either  suppress,  enhance  or  change  the  transition  from  a  commensurate  order  into  an 
incommensurate  one.  Thus  the  chemisorbed  metallic  surface  provides  a  unique  and  rich  system  for  the  study  of 
structural  phase  transitions.  In  the  case  of  incommensurate  order,  one  can  obtain  a  physical  realization  of  the 
X-V  model  The  effect  of  random  impurities  on  the  topological  order  in  this  case  has  been  investigated  via  field 
theoretical  methods. 

Principal  Investigator:  S.  C.  Ying  (Physics). 

Personnel:  K.  H.  Lau  (Physics). 

Publications:  "Effect  of  H-Adsorption  on  the  Displacive  Transition  of  W(001)  Surface,"  K.  H.  Lau  and  S.  C.  Ying, 
Phys.  Rev.  Lett.  44  (1980),  1222. 

“Hydrogen  Chemisorption  and  the  Displacive  Transition  of  W(100)  Surface,"  K.  H.  Lau  and  S.  C.  Ying,  in  Ordering 
in  Two  Dimensions,  North  Holland  (1980). 

"Effect  of  Substrate  Distortion  on  the  Ordered  Phases  in  Chemisorption  Systems."  S.  C.  Ying,  invited  talk  deliv¬ 
ered  at  the  27th  National  Symposium  of  A.V.S.  (1980)  and  published  in  Journal  of  Vacuum  Science  and  Technol¬ 
ogy 

“Effect  of  a  Random  Symmetry-Breaking  Field  on  the  Topological  Order  in  Two  Dimensions,”  A.  Houghton, 
R.  D.  Kenway,  and  S.  C.  Ying,  Phys.  Rev.  B  23  (1981),  289. 

Supported  by  the  Materials  Research  LaboratoryiNSF  and  the  National  Science  Foundation. 


Electronic  Structure  of  Defects  and  Surfaces  in  Simple  Metals 

The  density  matrix  variational  method  has  been  applied  to  the  study  of  low  index  surfaces  of  alkali  metals.  The 
three-dimensional  electronic  density  and  surface  potential  are  obtained  for  the  first  time  for  the  entire  series  of 
alkali  metals.  They  differ  significantly  from  the  corresponding  results  of  the  widely  used  jeilium  model. 

A  corresponding  study  has  been  made  for  a  monovacancv  in  Na.  The  discrete  lattice  nature  is  found  to  change 
the  Friedel  oscillations  around  the  vacancy  significantly  from  the  jeilium  picture. 

Principal  Investigator:  S.  C.  Ying  (Physics). 

Personnel:  K.  P.  Bohnen  (Physics)  and  M.  Lega  (Physics). 
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Publications:  “Electronic  Structure  of  K(001 )  Surface  by  the  Density  Matrix  Method."  K.  P.  Bohnen  and  S.  C.  Ying, 
Solid  State  Comm  30  (1979).  301 

‘Self-Consistent  Study  if  Surface  of  Simple  Metals  by  the  Density  Matrix  Method:  (100)  and  (110)  Surfaces  of 
Na,  K.  Rb  and  Cs,"  K  P  Bohnen  and  S.  C.  Ying,  Phys.  Rev.  B  22  (1980),  1806. 

Supported  by  the  Materials  Research  Laboratory/ NSF  and  the  National  Science  Foundation. 


Magnetic  Properties  of  Surfaces 

We  have  developed  a  new  method  for  calculating  the  Green's  functions  for  a  tight  binding  Hamiltonian.  This 
method  was  applied  to  the  study  of  the  magnetic  properties  of  surfaces  of  Ni.  V  and  Cr.  Substantial  surface 
enhancement  is  found  for  both  V  and  Cr.  in  agreement  with  recent  experimental  findings.  By  including  the  con¬ 
tribution  of  spin  fluctuations,  we  have  also  accounted  satisfactorily  for  the  transition  temperature  and  the  Curie- 
Weiss  like  susceptibility  of  these  systems  at  finite  temperatures. 

In  addition,  we  have  investigated  the  surface  magnetization  profile  of  an  Ising  type  magnet  near  the  critical  temp¬ 
erature. 

Principal  Investigator:  S.  C.  Ying  (Physics). 

Personnel:  D.  R.  Grempel  (Physics). 

Publications:  “Spin  Fluctuation  at  the  Surface  of  Vanadium,"  D.  R.  Grempel  and  S.  C.  Ying,  Phys.  Rev.  Lett.  45 
(1980).  1018. 

“Singularities  in  Surface  Magnetization  Profile,"  D.  R.  Grempel,  A.  Houghton,  and  S.  C.  Ying,  Phys  Letters  A  78 
(1980),  295. 

Supported  by  the  Materials  Research  Laboratory! NSF  and  the  National  Science  Foundation. 


The  Scattering  and  Ionization  of  Atoms  at  Surfaces 

The  first  part  of  the  study  of  the  efficiency  with  which  atoms  are  ionized  on  surfaces  is  nearly  complete.  Atomic 
beams  of  K,  Na,  Li,  and  Tl  have  been  directed  at  surfaces  of  Pt  and  Si,  and  the  variation  of  the  yields  of  scattered 
ions  with  temperature  has  been  measured.  The  results  show  that  the  statistical  ratio  of  ions  to  atoms  leaving  the 
surface  is  found  except  for  some  special  cases.  These  deviations  lead  to  conclusions  about  the  state  of  the  sur¬ 
face.  Analysis  of  the  composition  of  the  surface  by  Auger  spectroscopy  can  be  done  between  the  ionization 
measurements  to  check  for  the  presence  of  impurities.  Preliminary  measurements  on  an  insulator,  Al203,  sug¬ 
gest  that  the  statistical  ratio  is  no  longer  produced,  so  the  analysis  of  results  will  be  able  to  give  information  on 
kinetic  as  well  as  equilibrium  processes. 

The  apparatus  originally  chosen  for  measurements  of  the  elastic  scattering  of  rare  gases  from  the  surfaces  of 
crystals  of  graphite  has  been  shown  to  be  unsuitable  for  the  purpose.  The  lack  of  enough  differential  pumping  and 
the  relatively  high  ultimate  pressure  were  the  main  problems.  A  new  apparatus  has  been  designed,  the  compo¬ 
nents  have  been  built,  and  the  assembly  is  now  being  finished.  The  angular  distributions  of  elastically  scattered 
atoms,  that  the  new  apparatus  should  yield,  promise  to  give  useful  information  about  the  structure  of  the  surface. 

Principal  Investigator:  E.  F.  Greene  (Chemistry). 

Personnel:  A.  LeGrand  (Physics),  J.  Keeley  (Chemistry),  and  M.  Pickering  (Chemistry). 

Publication:  ‘Activation  of  Chemical  Reaction  by  Impact  of  Molecules  on  a  Surface,"  M.  S.  Connolly, 
E.  F,  Greene,  C.  Gupta,  P.  Marzuk,  T.  M.  Morton,  C.  Parks,  G.  Staker,  J.  Phys.  Chem.  85  (1981).  In  press. 
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Inelastic  Scattering  of  Diatomic  Molecules  from  Solid  Surfaces 

As  noted  in  our  last  report,  the  presence  of  scattered  molecules  in  the  target  chamber  has  increased  the  problem 
of  background  interference  so  that  we  are  unable  to  distinguish  between  background  absorption  and  absorption 
from  the  scattered  beam.  We  are  at  present  developing  pulsed  beam  sources  so  that  we  can  eventually  eliminate 
this  problem.  In  the  interim  we  have  been  investigating  a  new  tool  for  scattering  studies.  We  have  investigated  a 
new  universal  detector  that  can  serve  as  a  beam  detector  in  studies  of  gas-solid  scattering.  This  device  should  be 
particularly  suitable  for  use  with  the  pulsed  valve  now  on  order,  although  it  will  function  when  an  amplitude  modu¬ 
lated.  continuous  beam  is  used. 

Principal  Investigator:  G.  J.  Diebold  (Chemistry). 

Personnel:  J.  Hayden  (Chemistry). 

Publication:  "Direct  Detection  of  Momentum  Flux  in  Atomic  and  Molecular  Beams,”  J.  G.  Choi,  J.  Hayden, 
M.  O'Connor,  and  G.  J.  Diebold,  to  be  published. 

Supported  by  the  Materials  Research  Laboratory  I NSF  and  the  Petroleum  Research  Fund. 


Surface  Complex  Models  on  Rh-PFSA  lonomers 

The  studies  of  reactions  of  small  molecules  with  metal  ions  supported  in  PFSA  ionomers  are  designed  to  help 
elucidate  the  structure  and  bonding  of  chemisorbed  species  on  metals  and  supported  metals.  Our  recent  work 
with  Rh-PFSA  has  established  a  clear  connection  between  these  systems  and  those  known  to  be  catalytically 
active.  Chemisorption  of  CO  on  Rh  supported  by  Al203  or  Si02,  for  example,  is  a  widely  studied  step  in  the  pro¬ 
cesses  catalyzed  by  supported  Rh.  A  range  of  vibrational  bands  for  such  sorbed  CO  molecules  typically  is 
observed  and  assigned  to  CO  on  metallic  Rh.  Our  work  on  Rh-PFSA  reacted  with  CO  or  CO/H2  shows,  however, 
that  while  all  of  the  features  seen  on  the  active  catalysts  are  observed  in  this  system  (as  well  as  another  for 
weakly  chemisorbed  CO  on  Rh(lll)),  most  of  the  previous  interpretations  must  be  changed.  The  two  strong  fea¬ 
tures  commonly  seen  with  CO  on  Rh(AI203)  are  not  due  to  CO  on  Rh(O)  at  all,  but  are  due  to  oxidized  rhodium 
(Rh(l)),  and  two  higher  frequency  features  are  due  to  CO  on  Rh(lll)  sites.  The  lower  frequency  (ca  1850-1930 
cm  ')  features  for  CO  on  dispersed-Rh  catalysts  are  not  due  to  a  distribution  of  bridging  CO  structures  alone,  but 
to  a  distribution  of  them  and  of  Rh  oxidation  states.  Finally,  those  features  which  are  characteristic  of  CO  on 
Rh(O)  centers  are  observed  in  Rh-PFSA.  Moreover,  both  the  high  temperature  (ca  300°K)  and  low  temperature 
(ca  100-150°K)  forms  of  CO  on  Rh(0)  can  be  formed  through  different  chemical  routes  with  Rh-PFSA. 

The  work  on  this  and  related  M-PFSA  systems  continues  in  order  to  model  effectively  and  study  chemisorption 
species  on  a  variety  of  metals.  Preparation  In  this  way  permits  studies  which  are  not  possible  on  the  catalysts 
themselves. 

Principal  Investigator:  W.  M.  Risen,  Jr.  (Chemistry). 

Personnel:  D.  M.  Barnes  (Chemistry). 

Publication:  ‘Reactions  of  Supported  Rhodium  Species  with  CO:  Comparison  of  the  Rh-PFSA  System  with  Rh  on 
AljOj  or  SIOj,”  W.  M.  Risen,  Jr.  and  D.  M.  Barnes,  J.  Catalysis,  in  submission,  1981. 
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SECTION  5 


Lower  Dimensional  Systems 


Pseudo  One  Dimensional  Conductors 

The  effect  of  impurities  and  the  coupling  of  several  chains  via  phonons  for  quasi-one-dimensional  systems  has 
been  considered.  For  the  study  of  electron  localization  and  delocalization  (1 ),  we  considered  the  large  momentum 
(2k()  scattering  of  the  electrons  by  impurities  in  addition  to  the  usual  backward  scattering  interactions.  This  prob¬ 
lem  bears  a  very  strong  mathematical  similarity  to  the  case  of  coupling  of  several  chains  via  phonons  (2,3) 
(modulo  the  randomness  of  the  impurities).  For  the  latter  problem,  we  consider  BS  model  on  each  chain  and  the 
large  momentum  part  of  the  electron-phonon  coupling.  Mandelstam  operators  (generalized  to  the  case  of  several 
chains)  are  used  to  bosonize  the  two  problems.  The  resulting  theories  are  formally  related  to  the  SU(2)  Massive 
Thirring  model  (on  each  chain)  with  the  proper  renormalizations.  The  mass  gaps,  the  ranges  in  which  they  exist 
and  thus  the  ranges  of  localization  and  delocalization  are  obtained  by  using  a  method  developed  earlier  (4).  We 
found  that  these  ranges  depend  both  on  2V-Un  and  Un  in  contrast  to  some  other  results  (5).  We  also  con¬ 
cluded  that  the  transition  (from  the  massless  to  the  massive  phases)  has  an  XY-like  behavior  and  that  the  prob¬ 
lem  is  equivalent  to  two  coupled  Coulomb  gases.  Methods  of  calculation  of  the  correlation  and  the  free  energy  are 
discussed. 

We  found  that  the  expression  for  the  periodic  boundary  conditions  (PBC)  obtained  in  the  diagonalization  of  ID 
theories  by  the  Bethe  Ansatz  method  contains  all  the  physics  of  the  problem  and  that  they  can  be  used  to  relate 
and  classify  various  one-dimensional  and  two-dimensional  models  (5).  In  addition  to  the  exact  spectrum,  we 
obtain  an  infinite  number  of  conserved  charges,  the  S-Matrix,  and  various  correlation  functions.  We  are  then  able 
to  take  the  continuum  limit  of  lattice  theories  without  any  ambiguities.  It  is  found  that  two  different  limits  of  the 
spin- 1/2  XYZ  chain  leads  the  Bose-Gas  and  the  Massive  Thirring  Model  (6).  In  view  of  the  results  obtained  in  this 
investigation  involving  1 D  and  2D  models,  it  is  concluded  that  for  a  given  internal  symmetry  there  is  only  one  fun¬ 
damental  model  and  only  one  Bethe  Ansatz.  The  others  are  obtained  by  various  limits  or  transformations.  We 
find  for  example  that  the  “gauge"  equivalence  of  some  1 D  models  is  realized  as  a  mapping  of  the  kernel  and/or 
the  inhomogeneous  parts  of  these  integral  equations.  These  considerations  allow  the  classification  of  one-di¬ 
mensional  or  two-dimensional  models. 
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Surface  Excitations  in  Conducting  Solids 

In  conducting  solids  in  which  the  electron  density  does  not  drop  abruptly  from  its  bulk  value  to  zero  at  the  surface, 
a  new  class  of  surface  modes  called  “higher  multipole  excitations"  can  occur.  We  have  investigated  these  modes 
for  the  case  of  a  degenerate  polar  semiconductor  with  a  surface  depletion  layer.  The  higher  multipole  excitations 
interact  with  the  polar  phonons  giving  rise  to  “multipolariton  modes"  of  the  surface  depletion  layer.  The  dispersion 
of  these  modes  has  been  studied  for  a  simple  two-step  model  in  the  long  wavelength  non-retarded  limit. 

Principal  Investigator:  J.  J.  Quinn  (Physics). 

Personnel:  K.  S.  Yi  (Physics). 

Publications:  “A  Simplified  Treatment  of  Exchange  and  Correlation  in  Semiconducting  Surface  Inversion  Layers," 
J.  J.  Quinn,  R.  K.  Kalia,  G.  Kawamoto,  and  S.  C.  Ying,  Solid  State  Comm.  34  (1980),  423. 

“Surface  Helicon  Waves  in  the  Faraday  Geometry;  Do  They  Exist,”  J.  J.  Quinn  and  P.  Halevi,  Solid  State  Comm. 
33  (1980),  467 

'Theory  of  the  Effect  of  a  dc  Current  on  Surface  Polaritons  in  n-Type  Silicon,”  J.  J.  Quinn,  B.  G.  Martin  and 
R.  F.  Wallis,  to  appear  Surf.  Science. 

Theory  of  the  Effect  of  a  dc  Current,  Including  the  Electron  Thermal  Pressure  Gradient  on  Surface  Polaritons  in 
n-Type  Silicon,”  J.  J.  Quinn,  B.  G.  Martin,  and  R.  F.  Wallis,  Proc.  Int.  Conf.  on  Surfaces  and  Thin  Films,  Cannes, 
1980. 

“Exchange-Correlation  Effects  in  Silicon  Inversion  Layer:  Valley  Occupancy  Phase  Transitions,"  J.  J.  Quinn, 
W.  L.  Bloss,  S.  C.  Ying,  T.  Cole,  and  B.  D.  McCombe,  Proc.  Int.  Conf.  on  Semiconductor  Physics,  Kyoto,  1980. 

“Exchange  Correlation  Effects  in  Silicon  (111)  Inversion  Layers:  Stress  Enhanced  Valley  Occupancy  Phase 
Transition,"  J.  J.  Quinn,  W.  L.  Bloss,  and  S.  C.  Ying,  to  appear  Phys.  Rev. 

“Magneto-Plasma  Surface  Waves  of  a  Degenerate  Semiconductor  in  the  Faraday  Geometry  -  Effect  of  the  Pres¬ 
ence  of  a  Metallic  Screen,"  J.  J.  Quinn,  K.  S.  Yi,  and  P.  Halevi,  Phys.  Rev.  22  (1980),  6247. 

‘Surface  Modes  of  a  Degenerate  Polar  Semiconductor  with  a  Surface  Depletion  Layer,"  J.  J.  Quinn,  S.  Das 
Sarma.  and  A  Eguiluz,  to  appear  Solid  State  Comm. 

‘Properties  of  a  Surface  Inversion  Layer  in  a  Polar  Semiconductor,"  J.  J.  Quinn,  G.  Kawamoto,  and  W.  L.  Bloss,  to 
appear  Phys.  Rev. 

"Multipolaritons  of  a  Polar  Semiconductor  Surface  Depletion  Layer,"  J.  J.  Quinn,  S.  Das  Sarma,  and  A.  Eguiluz, 
Proc.  of  Int.  Conf.  on  Phys.  of  Compound  Semiconductors,  Williamsburg,  1981. 
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“Evidence  for  Valley  Occupancy  Phase  Transition  in  Si-(IOO)  Inversion  Layers,"  J.  J.  Quinn,  T.  Cole, 
B.  D.  McCombe,  and  R.  K.  Kalia,  submined  to  Phys.  Rev.  Leners. 

Supported  by  the  Materials  Research  Laboratory/NSF,  the  National  Science  Foundation,  and  the  Office  of 
Naval  Research. 


Low-Dimensional  Quantum  Field  Theories  and  Their  Application  to  Disordered  Sys¬ 
tems 

We  studied  the  zero  temperature  electrical  conductivity  of  a  disordered  solid  in  2  +  e  dimensions.  This  is  related  to 
the  critical  behavior  of  a  non-linear  a-model  on  the  non-compact  space  0(M,M)/0(M)x0(M)  in  the  limit  M^O.  We 
calculated  the  0- function  to  one  loop  order  for  general  Einstein  spaces,  deducing  from  it  that  all  electronic  states 
are  localized  in  2  dimensions.  Above  2  dimensions,  for  weak  disorder  there  are  mobility  edges,  but  these  merge 
above  a  critical  amount  of  disorder  and  all  states  become  localized. 

Guerin  and  Kenway  formulated  a  strong-coupling  expansion  in  the  inverse  fermion  propagator  for  the  N-species 
Gross-Neveu  model,  and  summed  it  to  all  orders  in  the  large-N  limit.  We  recovered  all  the  exact  results  previously 
derived  from  the  weak-coupling  expansion  and  were  able  to  compare  different  regularization  schemes. 

Houghton,  Ying,  and  Kenway  looked  at  the  effect  of  quenched  random  symmetry-breaking  fields  on  topological 
order  in  the  2  dimensional  XY  model.  This  was  to  understand  whether  substrate  impurities  interfere  with  a  Koster- 
litz-Thouless  phase  transition  in  the  absorbate.  We  used  a  renormalization  group  analysis  of  the  equivalent 
sine-Gordon  field  theory  to  show  that  the  transition  survives  provided  the  random  field  has  locally  three-fold  sym¬ 
metry  or  higher. 

Principal  Investigator:  R.  D.  Kenway  (Physics). 

Personnel:  A.  Houghton  (Physics),  S.  C.  Ying  (Physics),  F.  Guerin  (Physics),  A.  Jevicki  (Physics),  and 
A.  M.  Pruisken  (Physics). 

Publications:  “Non-Compact  a-Models  and  the  Existence  of  a  Mobility  Edge  in  Disordered  Electronic  Systems 
Near  Two  Dimensions."  A.  Houghton,  A.  Jevicki,  R.  D.  Kenway,  and  A.  M.  Pruisken,  Phys.  Rev.  Lett.  45  (1980), 
394. 

"A  Strong-Coupling  Expansion  for  Fermion  Field  Theories  and  the  Large-N  Limit  of  the  Gross-Neveu  Model," 
F.  Guerin  and  R.  D.  Kenway,  Nucl.  Phys.  0776  (1980),  168. 

“The  Effects  of  a  Random  Symmetry-Breaking  Field  on  Topological  Order  in  Two  Dimensions,”  A.  Houghton, 
R.  D.  Kenway,  and  S.  C.  Ying,  Phys,  Rev.  823  (1981),  298. 

Supported  by  the  Materials  Research  Laboratory/NSF. 


Resonance  Raman  Spectra  of  Highly  Conducting  Poly-p-phenylene  Complexes 

The  discovery  of  the  highly  conducting  charge-transfer  complexes  of  poly-p-phenylene  has  attracted  a  great  deal 
of  attention.  Through  chemical  doping  with  electron-acceptors,  such  as  AsFs,  or  donors,  such  as  Na,  the  electri¬ 
cal  conductivity  of  these  conjugated  ir-eiectron  systems  increases  from  insulating  to  metallic  values 
(1x10s  S  m1).  The  major  problem  in  studying  such  systems  as  these,  potyacetytenes,  or  others  is  in  discovering 
the  nature  of  the  polymer  after  doping— l.e.,  how  its  structure,  electron  states,  etc.  have  been  changed  by  doping. 
This  is  a  problem  because  vibrational  spectroscopic  methods  have  not  been  successful  previously  in  obtaining 
this  information  on  the  highly  conducting  forms. 

In  this  work  the  resonance  Raman  spectra  of  poty-p-phenyiene  (PPP)  and  its  highly  conductive  AsFg-  and  Na- 
doped  charge  transfer  complexes  have  been  obtained  and  shown  to  reveal  the  effects  on  the  polymer  caused  by 
doping.  There  is  selective  enhancement  of  different  vibrational  bands  associated  with  the  in-ring  and  bridge 
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bonds  of  the  different  samples.  One  key  result  is  that  doping  with  either  type  of  dopant  causes  the  C-C  bond  join¬ 
ing  the  phenyl  rings  to  become  stronger  and  at  least  one  of  the  in-ring  bonds  to  become  weaker,  so  all  C-C  bonds 
along  the  backbone  of  the  ir-electron  system  become  more  alike  than  they  are  in  the  undoped  PPP.  Another  is 
that  the  resonance  enhancements  show  the  relationship  between  changes  in  molecular  configuration  caused  by 
optical  transitions  in  the  donor-  or  acceptor-doped  PPP  and  the  vibrational  modes  of  the  polymer.  Finally,  the  data 
are  analyzed  in  terms  of  the  proposed  electronic  structures  of  pure  and  doped  PPP. 

Principal  Investigator:  W  M.  Risen,  Jr.  (Chemistry). 

Personnel:  C.  H.  Tzinis  (Chemistry)  and  E.  I.  Kamitsos  (Chemistry). 

Publication:  “Raman  Spectra  of  Highly  Conducting  Poly-p-phenylene  Charge  Transfer  Complexes,”  W.  M.  Risen, 
Jr.,  C.  H.  Tzinis,  and  R.  H.  Baughman,  Solid  State  Comm.,  in  submission. 

Supported  by  the  Materials  Research  Laboratory  I NSF  and  the  Office  of  Naval  Research. 


Lower  Dimensional  Studies 

This  period  involved  the  initiation  of  many  new  experiments.  We  have  begun  measurements  to  determine  if  the 
predicted  2-1  valley  transition  of  Bloss  and  Quinn  exists.  We  have  determined  new  information  on  valley  degen¬ 
eracies  in  (110)  and  (111)  surface  MOSFET  structures  in  Si.  We  see  a  transition  from  6  valleys  to  2  valleys  in  new 
(111)  MOSFETs  as  well  as  just  2  valleys  on  other  (111)  MOSFET  structures  (prepared  differently).  We  see  4  val¬ 
leys  on  (1 10)  surfaces  on  our  newer  samples.  At  sufficiently  high  magnetic  fields,  we  see  a  transition  to  1  valley 
structure. 

We  have  indications  that  high  pressure  oxidation  will  iead  to  a  better  Ge:GeO(  structure  to  study  two-dimensional 
electron  gases.  Work  continued  on  electrons  on  helium,  but  was  unsuccessful.  We  are  continuing  that  study  and 
the  plasma  oxidation  of  GaAs. 

Principal  Investigator:  P.  J.  Stiles  (Physics). 

Personnel:  E.  Crisman  (Physics),  C.  Huber  (Physics),  J.  Lee  (Physics),  K.  C.  Woo  (Physics),  and  R.  Zeller  (Engi¬ 
neering). 

Publication:  “Quasi-Two-Dimensional  Electrons  In  Surface  Space  Charge  Layers,"  P.  J.  Stiles,  to  be  published. 

Supported  by  the  Materials  Research  Laboratory/NSF,  the  National  Science  Foundation,  and  the  Office  of 
Naval  Research. 


SECTION  6 


Low  Temperature  Properties  of  Materials 


Study  of  Rotons  in  Liquid  Helium-Four 

The  theory  ot  roton  second  sound  in  superfluid  helium-4  is  considered.  We  first  derive  the  hydrodynamic  equa¬ 
tions  which  determine  the  motion  of  a  pure  roton  gas.  The  collisions  between  rotons  are  assumed  to  conserve 
energy,  momentum,  and  the  number  of  rotons.  The  hydrodynamic  equations  without  dissipative  terms  have  the 
same  form  as  the  corresponding  equations  for  an  ordinary  fluid.  When  lowest  order  dissipative  effects  are 
included,  it  is  found  that  there  are  some  extra  dissipative  processes  in  addition  to  heat  conduction  and  viscosity. 
We  next  consider  modifications  of  the  hydrodynamic  equations  which  occur  when  the  effects  of  collisions 
between  rotons  and  phonon  are  included,  and  when  there  are  some  roton-roton  collisions  in  which  the  roton 
number  changes.  It  is  found  that  at  the  temperatures  and  pressures  where  roton  second  sound  has  so  far  been 
observed,  the  wave  motion  is  best  described  as  isothermal  second  sound.  We  show  that  current  experimental 
data  indicate  that  in  nearly  all  collisions  between  rotons  ('99.97%)  the  total  number  of  rotons  is  conserved. 

Principal  Investigator:  H  J.  Maris  (Physics). 

Personnel:  R.  W.  Cline  (Physics). 

Publications:  “Velocity  of  Roton  Second  Sound,"  R.  W.  Cline  and  H.  J.  Maris,  in  Proceedings  of  the  Third  Interna¬ 
tional  Conference  on  Phonon  Scattering  in  Condensed  Matter,  ed.  by  H.  J.  Maris  (Plenum,  New  York,  1980),  p. 
161. 

“Roton  Second  Sound  and  Roton  Scattering,"  R.  W.  Cline  and  H.  J.  Maris,  to  appear  in  Phys.  Rev.  B„  April  1981 . 
Supported  by  the  Materials  Research  Laboratory/ NSF  and  the  National  Science  Foundation. 


Kapitza  Resistance 

We  have  carried  out  an  experiment  to  measure  the  Kapitza  resistance  between  liquid  and  solid  helium,  and  have 
obtained  some  preliminary  data.  This  experiment  is  important  because  the  surface  of  solid  helium  is  highly  perfect 
and  clean,  and  so  one  has  the  opportunity  to  study  the  Kapitza  resistance  in  the  absence  of  surface  damage  or 
dirt.  We  hope  to  make  more  measurements  and  analyze  the  data  during  the  next  few  months. 

Principal  Investigator:  H.  J.  Maris  (Physics). 

Personnel:  T.  E.  Huber  (Physics). 

Supported  by  the  Materials  Research  Laboratory/NSF  and  the  National  Science  Foundation. 


Energy-Selective  Phonon  Detector 

We  have  continued  our  experimental  study  of  phonon  detection  using  superconducting  zinc  films.  The  character¬ 
istics  of  the  detect or  have  been  explored  by  means  of  several  experiments  with  phonons  and  rotons  in  superfluid 
helium.  Measurements  of  the  threshold  energy  for  phonon  detection,  and  its  correlation  to  changes  in  the  i-V 
characteristics  with  magnetic  field  have  been  made.  A  tentative  explanation  of  the  detection  mechanism  in  terms 
of  phonon-activated  ftuxoid  motion  has  been  proposed. 


Principal  Investigator:  H.  J.  Maris  (Physics). 
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Personnel:  fl.  W.  Cline  (Physics). 

Publications:  “Energy-Selective  Detection  of  Phonons  by  Superconducting  Zinc  Films,"  R.  W.  Cline  and 
H.  J.  Maris,  to  be  published. 

“Observation  of  a  New  Type  of  Energy-Selective  Phonon  Detector,”  R.  W.  Cline  and  H.  J.  Maris,  in  Proceedings 
ot  the  Third  International  Conference  on  Phonon  Scattering  in  Condensed  Matter,  ed.  by  H.  J.  Maris  (Plenum, 
New  York,  1980),  p.  317. 

Supported  by  the  Materials  Research  Laboratory/NSF  and  the  National  Science  Foundation. 


ESR  of  Atomic  Hydrogen  at  4.2°K 

ESR  measurements  at  9GHz  have  been  performed  on  atomic  hydrogen  in  a  cavity  the  walls  of  which  are  coated 
with  solid  molecular  hydrogen  at  low  temperatures.  At  4.2°K  and  below,  a  broad  line  spectrum  is  observed  and 
ascribed  to  atoms  trapped  on  the  surface  of  solid  H2  at  fixed  sites.  The  broad  line  spectrum,  which  appears  in 
addition  to  the  narrow  line  spectrum  associated  H  atoms  in  the  gas,  has  several  notable  characteristics — a  g  shift, 
shift  in  hyperfine  interaction  and  line  width  that  are  compatible  with  atoms  in  surface  states.  The  absorption  spec¬ 
trum  persists  indefinitely  below  4°K  when  the  source  of  H  atoms  is  turned  off.  And  as  the  temperature  is  lowered, 
the  broad  spectrum  increases  in  intensity  as  the  narrow  spectrum  decreases.  Furthermore,  on  turning  off  the 
source  of  H  atoms  and  warming  to  4.2°K,  the  broad  spectrum  is  stronger  than  with  the  source  on  and  the  narrow 
gas  resonance  present. 

We  interpret  these  observations  as  indicating  that  the  broad  line  spectrum  originates  from  H  atoms  trapped  at 
fixed  sites  on  the  surface  of  molecular  H2.  Support  for  this  contention  comes  from  a  comparison  of  the  measured 
parameters  of  this  spectrum  with  those  found  by  other  workers  for  H  atoms  incorporated  in  solid  H2,  from  the 
difficulty  of  others  in  creating  atoms  in  solid  H2  by  deposition,  and  from  the  dependence  of  the  intensity  of  the 
broad  line  on  the  narrow  line.  Our  measurements  require  a  density  of  immobile  surface  states  of  between  1010 
and  lO"/cmJ,  suggesting  that  the  localized  sites  are  limited  to  defects  and  grain  boundaries  on  the  surface. 

Principal  Investigator:  G.  M.  Seidel  (Physics). 

Personnel:  R.  Mayer  (Physics)  and  J.  Cameron  (Physics). 

Supported  by  the  Materials  Research  Laboratory/NSF  and  the  National  Science  Foundation. 


Recovery  of  Solid  Helium-4  After  Plastic  Deformation 

Recovery  of  both  bcc  and  hep  phases  of  solid  helium-4,  after  being  subjected  to  plastic  deformation,  has  been 
investigated  by  measuring  attenuation  changes  as  a  function  of  time.  Bcc  and  hep  He4  show  distinctly  different 
behavior.  In  bcc  He4,  the  attenuation  changes  somewhat  upon  deformation  and  relaxes  within  a  time  of  order  of 
one  minute.  In  hep  He4,  attenuation  changes  considerably  upon  deformation  and  it  takes  a  time  of  order  of  10 
minutes  to  relax. 

Principal  Investigator:  C.  Elbaum  (Physics  and  Applied  Mathematics). 

Personnel:  H.  Kwun  (Materials  Research  Laboratory). 

Supported  by  the  Materials  Research  Laboratory/NSF. 
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Finite  Amplitude  Weve  Propagation  in  Solid  and  Liquid  Helium-4 

Amplitude  response  of  a  longitudinal  ultrasonic  wave  (10MHz)  was  investigated  in  solid  helium  (bcc  and  hep 
phases)  as  well  as  in  liquid  helium  (normal  and  superfluid  phases).  The  amplitude  response  curves  show  satura¬ 
tion  behavior  which  is  well  described  by  an  equation  for  amplitude  loss  based  on  second-harmonic  generation. 

Principal  Investigator:  C.  Elbaum  (Physics  and  Applied  Mathematics). 

Personnel:  H.  Kwun  (Materials  Research  Laboratory). 

Publication:  ‘Finite  Wave  Amplitude  Propagation  in  Solid  and  Liquid  Helium  4,"  A.  Hikata,  H.  Kwun,  and 
C.  Elbaum,  Phys.  Rev.  B21  (1980),  3932. 

Supported  by  the  Materials  Research  Laboratory/NSF. 


Nonlinear  Wave  Interaction  In  a  Dispersionless  Liquid 

in  addition,  nonlinear  interactions  of  two  ultrasonic  waves,  propagating  collinearly  in  a  dispersionless  liquid 
medium  were  investigated. 

Principal  Investigator:  C.  Elbaum  (Physics  and  Applied  Mathematics). 

Personnel:  H  Kwun  (Matenals  Research  Laboratory). 

Publication:  “Nonlinear  Interactions  of  Two  Ultrasonics  in  a  Dispersionless  Liquid  Medium,”  H.  Kwun,  A.  Hikata, 
and  C.  Elbaum,  J.  Appl.  Phys.  51  (1980),  4645. 

Supported  by  the  Materials  Research  Laboratory/NSF. 


Anharmonic  Properties  of  Liquid  and  Solid  Helium 

This  study  was  undertaken  in  order  to  determine  such  features  as  higher  (third)  order  elastic  constants,  wave 
coupling,  and  the  response  of  dislocations  (in  solids)  to  high  amplitude  ultrasonic  waves.  To  this  end  the  ampli¬ 
tude  response  (i.e.,  relation  between  input  and  output  amplitudes  of  a  wave  transmitted  through  a  medium)  of 
10  MHz  longitudinal  wave  is  investigated  in  hep  and  bcc  solid  4He  as  well  as  normal  fluid  and  superfluid  4He.  In 
all  four  phases  the  acoustic  saturation  in  the  amplitude  response  curve  was  observed.  The  results  are  well 
described  by  an  equation  for  amplitude  loss  based  on  second-harmonic  generation.  The  nonlinear  parameter  (j  of 
these  phases  was  evaluated.  It  is  found  that  the  magnitude  3  in  solid  4He  is  comparable  to  that  of  classical  solids. 
Dislocation  contributions  to  the  amplitude  response  were  also  investigated  and  were  ruled  out  as  a  significant 
source  of  the  observed  effects.  It  is  concluded,  therefore,  that  the  large  nonlinear  amplitude  response,  which 
leads  to  saturation  observed  in  solid  4He  should  be  attributed  to  the  unusually  large  compressibility,  i.e.,  small 
elastic  constants  (second  order),  which,  in  turn,  induces  in  solid  4He  much  larger  strain  amplitudes  for  a  given 
applied  stress  than  in  ordinary  solids. 

Principal  Investigator:  C.  Elbaum  (Physics  and  Applied  Mathematics). 

Personnel:  A.  Hikata  (Applied  Mathematics)  and  H.  Kwun  (Materials  Research  Laboratory). 

Publication:  ‘Finite  Amplitude  Wave  Propagation  in  Solid  and  Liquid  4He,”  A.  Hikata,  H.  Kwun,  and  C.  Elbaum, 
Phys.  Rev.  B1  (1980),  3932. 


Supported  by  the  Materials  Research  Laboratory/NSF. 
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Anharmonic  Effects  in  Ultrasonic  Propagation  Near  Second-Order  Phase  Transitions 

Ultrasonic  studies  near  second-order  phase  transitions  have  been  extended  here  to  sufficiently  high  amplitudes 
of  the  waves,  so  that  nonlinear  (anharmonic)  behavior  can  be  observed.  In  particular,  anharmonic  effects  were 
studied  by  measuring  the  temperature  dependence  of  second  harmonics  of  a  104MHz  ultrasonic  wave  generated 
near  the  x  transition  of  pressurized  liquid  4He.  The  anharmonic  coupling  coefficient  C  between  the  fundamental 
wave  and  the  second  harmonic  is  found  to  diverge  as  ~iti  '(t  =  (T/Tk)-i]  near  Tk.  A  phenomenological  relation 
was  developed  expressing  C  in  terms  of  known  relevant  thermodynamic  quantities.  This  relation  accounts  well  for 
the  experimental  results. 

Principal  Investigator:  C.  Elbaum  (Physics  and  Applied  Mathematics). 

Personnel:  A.  Hikata  (Applied  Mathematics)  and  H.  Kwun  (Materials  Research  Laboratory). 

Publication:  “Anharmonic  Effects  in  Ultrasonic  Propagation  Near  the  Lambda  Transition  of  4He,"  H.  Kwun, 
A.  Hikata,  and  C.  Elbaum,  J.  Low  Temperature  Physics  39  (1980),  497. 

Supported  by  the  Materials  Research  LaboratoryINSF. 
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SECTION  7 


Other  Research 


Polymers 

Recent  experiments  by  Liao  and  Morawetz'  show,  contrary  to  physical  intuition,  that  the  activation  energy  tor 
rotational  barrier  crossings  in  a  long-chain  molecule  in  solution  is  the  same  as  tor  an  analog  molecule  with  only  a 
single  rotational  barrier.  We  have  studied  this  process  by  computer  simulation,  using  the  technique  of  Brownian 
dynamics  applied  to  a  model  of  a  polymer  chain  regarded  as  a  system  of  linked  rigid  bodies.  The  computer  simu¬ 
lation  results  were  found  to  agree  qualitatively  with  the  experimental  results  and  provided  insight  into  the  reasons 
underlying  the  activation  energy  insensitivity. 


’T.-P.  Liao  and  H.  Morawetz.  Macromolecules  13  (1980),  1228. 

Principal  Investigator:  J.  H.  Weiner  (Engineering). 

Personnel:  D.  Perchak  (Physics). 

Publication:  "Activation  Energy  Insensitivity  to  Barrier-Crossing  Correlations  in  Long-Chain  Molecules,"  D.  Perc¬ 
hak  and  J.  H.  Weiner,  submitted  to  Macromolecules. 

Supported  by  Materials  Research  Laboratory/GRI. 


Dynamic  Mechanical  Response  of  Fiber-Reinforced  Materials 

A  new  program  has  been  started  to  study  the  dynamic  mechanical  response  of  fiber-reinforced  materials.  Their 
response  both  when  the  deformations  are  elastic  and  plastic  is  being  investigated,  and  at  the  same  time,  their 
fracture  strength  and  the  nature  of  crack  propagation  in  such  materials  is  being  observed.  This  program  is  both 
experimental  and  theoretical  and  is  being  carried  out  in  collaboration  with  Professor  A.  C.  Pipkin. 

The  specimens  used  in  these  studies  are  of  two  types.  First,  specimens  of  Scotchply  2001  (an  epoxy-resin  rein¬ 
forced  with  glass  fibers)  are  being  tested.  Two  types  of  Scotchply  sheet  have  been  obtained,  in  one  the  glass 
fibers  are  all  oriented  parallel  to  each  other,  in  a  single  direction;  in  the  other  they  are  in  the  form  of  a  square 
mesh,  one  half  being  in  one  direction  and  the  other  half  in  a  direction  perpendicular  to  it. 

The  other  type  of  specimen  which  is  being  studied  are  beams  of  lead  and  tin,  which  have  been  reinforced  along 
the  beam  direction  either  by  steel  piano  wires  or  by  wires  of  phospher-bronze.  Three  types  of  experimental  test 
have  been  performed.  These  are  (a)  small  elastic  deformation,  (b)  large  plastic  deformation,  and  (c)  crack  propa¬ 
gation  studies. 

The  Scotchply  epoxy  resin  has  been  used  in  studies  (a)  and  (c)  and  the  metal  specimens  in  studies  (a)  and  (b). 
The  experimental  set-ups  for  the  three  types  of  investigation  have  been  as  follows:  for  (a)  beams  have  been  set 
into  longitudinal,  flexural  and  torsional  resonance  and  from  the  resonant  frequencies,  the  elastic  moduli  and  the 
mode  of  deformation  has  been  inferred;  for  (b)  small  charges  of  high  explosive  and  high  speed  Impacts  in  the 
"High  G"  testing  machine  have  been  carried  out;  while  for  (c)  precracks  have  been  inserted  In  the  specimens 
which  have  been  subjected  to  impact  loads. 

Interesting  results  have  been  obtained  in  all  these  investigations  although  they  are  all  in  the  early  stages  of  study. 
In  particular  it  has  been  shown  that  when  beam  specimens  with  their  fibers  oriented  parallel  to  the  axis  of  the 
beam  are  set  into  flexural  oscillation  and  the  anisotropy  of  the  material  is  very  high  (l.e.,  the  ratio  E/G  is  very  large; 
E  is  the  effective  longitudinal  modulus  of  the  fibers  and  G  is  the  shear  modulus  of  the  matrix)  the  dynamic 


mechanical  response  of  the  material  corresponds  to  Q,  the  shear  modulus  of  the  material,  rather  than  E,  the 
extensions!  modulus  ot  the  fibers.  (This  is  in  agreement  with  the  prediction  of  M.  Sayir  who  very  recently  pub¬ 
lished  a  paper  on  this  effect  Ing.  Archiv.  49  (1980),  309-30.) 

High  speed  cine-photographs  of  transverse  dynamic  impacts  on  fiber  reinforced  beams  are  being  used  to  test  the 
accuracy  of  the  theoretical  predictions  of  A.  J.  M.  Spencer  and  his  group  on  the  response  of  such  materials. 
Finally,  measurements  have  been  started  on  the  dynamic  stress  intensity  factors  which  occur  when  brittle  frac¬ 
tures  are  propagated  through  sheets  of  fiber  reinforced  materials. 

Work  has  also  continued  in  a  program,  started  earlier,  on  the  anomalous  torsional  response  of  tubes  of  rubber 
subjected  to  large  extensional  strains.  By  studying  the  change  in  anomalous  response  of  the  behavior  of  speci¬ 
mens  of  a  number  of  different  rubbers  at  different  temperatures,  information  has  been  obtained  on  the  relation  of 
the  crystallization  of  the  rubber  to  its  biaxial  mechanical  response. 

Principal  Investigator:  H.  Kolsky  \ Applied  Mathematics). 

Personnel:  S.  Guofang  (Applied  Mathematics),  J.  Mosquera  (Engineering),  and  P.  Russo  (Applied  Mathematics). 

Publications:  "The  Role  of  Stress  Waves  in  Penetration  Processes,”  H.  Kolsky,  chapter  in  Ballistic  Materials  and 
Penetration  Mechanics,  ed.  by  R.  Laible  (Elsevier  Scientific  Publishing  Co.,  1980),  pp.  185-223. 

“The  Viscoelastic  Response  of  Highly  Stretched  Rubber  Specimens  to  Superposed  Torsional  Deformation," 
H.  Kolsky  and  A.  C.  Pipkin.  Ing.  Archiv.  49  (1980),  337-45. 

“The  Anomalous  Mechanical  Response  of  Highly  Stretched  Tubes  of  Natural  Rubber  under  Conditions  of  Biaxial 
Loading,”  J.  Y.  Liu,  ARPA  Technical  Report  AM  49  (1979). 
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Photoelectrolysis 

Preparation  and  Comparison  of  the  Photoelectronic  Properties  of  Sr.Nb,0T  and 
B«0  s  Sr0  5  NbjOs 

The  two  alkaline  earth  niobates  Sr2Nb207  and  Ba05  Sr05  Nb206  have  been  prepared,  their  electronic  proper¬ 
ties  measured,  and  their  photoresponses  compared.  The  indirect  band  gap  in  Sr2Nb207  is  3.86  eV  compared 
with  3.38  eV  for  Ba05  SrQ5  Nb2Og  Hence,  photoanodes  composed  of  Sr2Nb207  respond  to  much  less  of  the 
"white"  light  spectrum  than  those  made  from  Ba05  SrQ5  Nb2Oe.  Nevertheless,  their  electrical  outputs  at  an 
anode  potential  of  0.8  eV  with  respect  to  SCE  in  0.2  M  soldium  acetate  under  “white"  xenon  arc  irradiation  of 
1.25  W/cm?  are  comparable 

Preparation  and  Photoelectronic  Properties  of  FeNb04 

The  orthorhombic  a-Pt>02  phase  of  FeNbO,  was  prepared  and  its  photoelectronic  properties  measured:  Sin¬ 
tered  disks  were  shown  to  be  n  type  and  gave  a  resistivity  of  40  fl-cm.  Measurements  of  the  photoresponse  gave 
a  flat-band  potential  between  0.1  and  0.4  V  vs  SCE  at  a  pH  of  8.5  and  an  optical  band  gap  of  2.08(2)  eV.  Several 
high-energy  band  gaps  at  2.68(2),  3.24(2)  and  4.38(2)  eV  were  also  determined.  There  appears  to  be  an 
enhancement  of  the  quantum  efficiency  due  to  the  presence  of  [FeCy  active  centers  while  retaining  the  funda¬ 
mental  characteristics  of  the  [NbOJ  octahedra. 

Preparation  and  Photoelectronic  Properties  of  the  System  Cd2Qef.)ISi)[04 

Members  of  the  system  CdjGe^Si^O,  where  0sxs0.4  have  been  prepared.  These  compounds  were 
observed  to  crystallize  with  the  olivine  structure,  space  group  Pbnm.  The  resistivity,  Hall  mobility,  flat-band  poten¬ 
tial,  band  gaps,  and  stability  were  determined  as  functions  of  composition.  The  variation  of  these  photoelectronic 
properties  can  be  attributed  to  the  reduction  ot  the  cell  parameters  with  increasing  silicon  substitution.  Substitu¬ 
tion  of  silicon  for  Germanium  reduces  the  loss  of  photocurrent  from  25%  after  one  hour  for  *  =  0.0  to  only  6%  after 
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twenty-two  hours  for  x  =  0.4. 

Principal  Investigator:  A.  Wold  (Chemistry). 

Personnel:  See  below  (Hydrodesulfurization  Catalysts). 
Publications:  See  below  (Hydrodesulfurization  Catalysts). 


Hydrodesulfurization  Catalysts 

The  Preparation  and  Properties  of  Poorly  Crystallized  CoS2  and  RuS2 

Amorphous  or  poorly  crystallized  samples  of  RuS2  and  CoS.,  were  made  by  sulfurizing  (NH4)RuCI6  and 
(NH3)CoCI2  at  low  temperatures  with  H2S.  The  amorphous  products  were  annealed  under  various  conditions 
and  the  resulting  degree  of  crystallinity  was  determined.  For  RuS2,  the  average  crystallite  sizes  ranged  from  26A 
for  the  poorlv  crystallized  material,  to  greater  than  2000A  for  the  well-crystallized  material.  CoS2  average  crystal¬ 
lite  sizes  ranged  from  130A  to  greater  than  2000A.  The  stoichiometry,  appearance  of  the  crystalline  pyrite 
phases,  cell  parameters,  and  magnetic  properties  were  studied  as  a  function  of  annealing  temperature. 

Principal  Investigator:  A.  Wold  (Chemistry). 

Personnel:  A.  Wold  (Chemistry),  K.  Dwight  (Chemistry),  W.  J.  Croft  (Chemistry),  J.  Hormadaly  (Chemistry), 
S.  N.  Subbarao  (Chemistry),  R.  Kershaw  (Chemistry),  B.  Khazai  (Chemistry),  J.  D.  Passaretti  (Chemistry), 
J.  Koenitzer  (Chemistry),  R.  R.  Chianelli  (Exxon  Research  &  Development),  and  V.  S.  Nguyen  (Chemistry). 

Publications:  “High  Pressure  Phase  Transformation  of  Platinum  Sulfide,"  A.  Wold,  R.  Collins,  R.  Kaner,  P.  Russo, 
and  D.  A  vignant,  I norg.  Chem.  18  (1979),  727. 

“The  Electrical,  Optical  and  Photoconducting  Properties  of  Fe2xCrx03  (0  x  0.47),"  A.  Wold,  P.  Merchant, 
R.  Collins,  R.  Kershaw,  and  K.  Dwight,  J.  Solid  State  Chem.  27  (1979),  307. 

“Interband  Transitions  of  Semiconducting  Oxides  Determined  from  Photoelectrolysis  Spectra,”  A.  Wold, 
F.  P  Koffyberg,  and  K.  Dwight,  Solid  State  Comm  30  (1979),  433. 

“The  Preparation  of  Amorphous  RuS2  at  Ambient  Temperature,”  A.  Wold,  J.  D.  Passaretti,  and  R.  C.  Collins, 
Materials  Res.  Bull.  14  (1979),  1167 

‘Preparation  and  Photoelectric  Properties  of  n-Type  Cd2Ge04,”  A.  Wold,  V.  S.  Nguyen,  S.  N.  Subbarao, 
R.  Kershaw,  and  K.  Dwight,  Materials  Res.  Bull.  14  (1979),  1535. 

‘Electrical  and  Optical  Properties  of  High-Purity  p-Type  Single  Crystals  of  GeFe  O.,"  A.  Wold,  P.  Strobel,  and 
F.  P.  Koffyberg,  J.  Solid  State  Chem.  31  (1980),  209. 

“Preparation  and  Properties  of  a  Copper  Niobium  Oxyfluoride,"  A.  Wold,  B.  Khazai,  K.  Dwight,  and  E.  Kostiner, 
Inorg.  Chem.  19  (1980),  1670. 

“Preparation  and  Comparison  of  the  Photoelectric  Properties  of  Sr2Nb207  and  Ba05Sr05Nb2O8,”  A.  Wold, 
J.  Hormadaly,  S.  N.  Subbarao,  R.  Kershaw,  and  K.  Dwight,  J.  Solid  State  Chem.  33  (1980),  27. 
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Investigations  of  Solar  Cells  Based  on  Chalcopyrite  Semiconductors  and  Their  Alloys 

In  the  previous  MRL  Annual  Technical  Report,  we  summarized  the  reasons  for  the  interest  in  the  chalcopyrite 
semiconductors,  of  which  CulnSe2  is  a  prototype,  for  solar  cell  applications.  Briefly,  it  has  already  been  shown 
that  photovoltaic  devices  that  can  convert  somewhat  more  than  10%  of  the  sunlight  incident  on  them  into  electric¬ 
ity  can  be  fabricated  from  single  crystals  of  CulnSe2  and  that  cells  having  efficiencies  up  to  6%  can  be  fabricated 
from  thin  films  of  CulnSe2.  Our  group  at  Brown  has  focused  its  attention  on  alloys  of  these  ternary  chalcopyrites 
like  CulnyGa  Se^Te^,  in  which  one  can  prepare  semiconductors  having  a  wide  range  of  bandgaps.  Such 
materials  can  be  combined  with  a  wide  bandgap  semiconductor  like  CdS  to  produce  heterojunction  photovoltaic 
cells  in  which  there  is  no  lattice  mismatch  and  in  which  intrinsic  interface  stale  concentrations  should  have  a  neg¬ 
ligible  concentration.  Solar  cells  made  from  such  pairs  of  lattice  matched  semiconductors  should  attain  higher 
efficiencies  than  cells  like  the  n-CdS/p-CulnSe2  cells  in  which  lattice  mismatch  is  about  1%. 

During  the  period  1  July  1979  through  30  June  1980,  we  have  continued  our  work  on  these  materials  and  on  solar 
cells  fabricated  from  them.  We  have  continued  to  synthesize  powders  of  four  and  five  element  alloys  of 
Cu<Ag(1  >(lnS2iSe2(,.2)  and  of  CulnyGa^,  Se2zTe2  .  Targets  for  rf-sputtering  were  prepared  from  these 
synthesized  powders.  Thin  films  of  the  Cu-ln-Ga-Se-Te  system  were  deposited  by  sputtering  from  these  targets. 
The  alloys  chosen  for  sputtering  had  forbidden  energy  bandgaps  of  about  1 .25  eV,  1 .35  eV  and  1 .45  eV,  all  with 
the  same  lattice  constant  as  CdS.  Optical  transmission  measurements  on  these  films  indicated  that  the  bandgaps 
were  close  to  their  anticipated  values.  X-ray  diffraction  studies  showed  that  they  had  either  the  sphalerite  or  the 
chalcopyrite  structure.  Analysis  of  EDAX  spectra  indicated  that  the  films  were  essentially  stoichiometric.  SEM 
studies  showed  that  grain  sizes  of  1  um  were  achievable;  such  grain  sizes  are  compatible  with  acceptable 
efficiency  photovoltaic  cells.  Devices  made  from  such  films  by  covering  them  with  evaporated  films  of  CdS  had 
open  circuit  voltages  of  0.5  V  but  their  other  photovoltaic  properties  were  poor  so  that  their  overall  solar  energy 
conversion  efficiency  was  low.  Identification  of  the  factors  limiting  the  photovoltaic  characteristics  of  these  films  is 
the  target  of  continuing  research  in  this  area. 

Films  of  the  ternaries  and  ternary  alloys  were  also  deposited  by  chemial  spray  pyrolysis,  a  deposition  technique 
which  is  even  more  attractive  from  an  economic  point  of  view  than  rf-sputtering.  Briefly,  to  deposit  films  by  this 
method  an  aqueous  solution  containing  a  stoichiometric  proportion  of  the  ions  in  the  desired  material  is  sprayed 
as  a  fine  mist  over  a  glass  substrate  which  is  heated  to  a  temperature  between  200°C  and  400°C.  When  the  mist 
strikes  the  substrate,  the  desired  ions  deposit  on  the  glass  and  form  the  compound  and  the  other  products  of  the 
reaction  volatilize.  We  have  prepared  thin  films  of  CulnS2,  CulnSe2,  AglnS2,  and  of  four  and  five  element  alloys 
of  the  Cu-Ag-ln-S-Se  system  by  this  method  Photovoltaic  cells  have  been  made  from  these  films  by  evaporating 
CdS  over  them 

We  have  also  continued  working  on  Semiconductor-Insulator-Semiconductor  (SIS)  solar  cells.  Large  area  (4 
cm2)  devices,  having  solar  energy  conversion  efficiencies  up  to  1 4%,  have  been  made  by  depositing  a  thin  layer 
of  indium-tin-oxide  (ITO)  over  single  crystal  silicon.  The  ITO  was  deposited  by  argon-beam  sputtering.  The  per¬ 
formance  of  these  SIS  devices  is  comparable  to  that  of  single  crystal  solar  cells  in  which  the  p/n  junction  is  formed 
by  diffusion. 

We  have  also  completed  theoretical  calculations  on  the  ultimate  efficiency  which  can  be  expected  from  new 
design  solar  cells  which  employ  minority  carrier  mirrors  and  optical  mirrors.  A  minority-carrier-mirror  (MCM)  is  an 
electrostatic  barrier,  internal  to  the  photovoltaic  cell,  at  which  minority  carriers  are  reflected  back  toward  the 
charge-separating  p/n  junction.  Such  MCM’s  can  increase  the  open  circuit  voltage  of  PV  cells  if  the  separation 
between  the  p/n  junction  and  the  MCM  is  less  than  the  minority  carrier  diffusion  length.  Our  calculations  show  that 
the  efficiency  of  thin  film  direct  gap  photovoltaic  cells,  specifically  CulnSe2  cells,  can  in  principle  attain  a  value  of 
26%  in  cells  whose  total  thickness  lies  between  one  and  two  microns.  The  insight  provided  by  these  calculations 
is  guiding  our  experiments;  we  are  trying  to  realize  the  benefits  of  MCM's  in  ternary  semiconductor  solar  cells 
fabricated  by  rf-sputtering. 

Finally,  we  have  been  exploring  application  of  these  theoretical  concepts  to  germanium  cells  intended  for  thermo- 
photovoltaic  energy  converters  (TPV)  which  are  of  Interest  to  the  Naval  Research  Laboratories.  As  the  experi¬ 
mental  part  of  this  effort  we  have  measured  the  minority  carrier  recombination  velocity  at  the  interface  between 
germanium  crystals  and  a  native  oxide,  produced  by  high  pressure  oxidation,  by  the  transient  capacitance 
method. 
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Nonequilibrium  Plasmas  In  Semiconductors 

Recently  developed  picosecond  infrared  techniques  have  been  applied  to  the  study,  for  the  first  time,  of  the  real 
time  dynamics  of  dense,  highly  degenerate,  transient  electron  hole  systems  in  PbTe,  inSb,  and  Hg  Cd„Te.  In 
these  experiments  an  intense  ultrashort  pulse  generated  a  high  density  electron-hole  gas  (>10”cm  *)  through 
interband  absorption  Aw>>  E  .  The  optical  properties  of  the  excited  system  were  then  examined  by  picosecond 
infrared  radiation  at  10.6  pin  and  5.3  pm  wavelengths  corresponding  to  probing  of  intraband  and  near  band  edge 
interband  transitions,  respectively.  From  time  resolved  decay  measurements  we  have  determined  concentration 


dependent  recombination  rates,  dominated  most  likely  by  rapid  Auger  processes.  For  example,  a  characteristic 
decay  time  of 5  40  psec  was  measured  in  PbTe  at  room  temperature  at  an  excess  density  of  10'9cm3  Meas¬ 
urements  for  InSb  and  Hg,  .Cdje  have  also  yielded  data  for  picosecond  relaxation  rates  at  the  high  excess 
densities  and  at  varying  temperatures  In  interpreting  the  results,  initial  attempts  have  been  made  to  include  in  the 
dielectnc  modeling  the  contnbution  by  many-body  effects,  which  appear  to  be  quite  significant  in  our  experimental 
conditions  In  this  we  have  been  aided  by  experimental  measurements  of  the  “dynamic"  plasma  edge  (o>  =  a>p(t)) 
and  time  varying  interband  transition  rates  (absorption  and  gain)  at  the  probe  beam  frequencies. 

Elsewhere  we  have  theoretically  examined  the  guiding  of  millimeter  waves  as  interface  electromagnetic  waves 
(surface  plasmons)  in  moderately  doped  semiconductors  of  high  material  quality.  Specifically,  we  have  calculated 
the  propagation  characteristics  and  modulation  of  surface  plasmons  in  n-GaAs.  Our  results  suggest  that  useful 
guiding  and  control  of  millimeter  wave  signals  by  these  excitations  may  be  possible. 

Principal  Investigator:  A.  V.  Nurmikko  (Engineering). 

Personnel:  W.  Lemley  (Engineering)  and  J.  Ajo  (Engineering). 

Publications:  "Study  of  Surface  Recombination  in  GaAs  and  InP  by  Picosecond  Optical  Techniques,” 
A.  V.  Nurmikko.  C.  A.  Hoffman,  and  H.  J.  Gerritsen,  J.  Appl.  Phys.  51  (1980),  1603. 

"Generation  of  Ultrashort  Pulses  in  Synchronous  Pumping  of  Near-Millimeter  Wave  Lasers,”  A.  V.  Nurmikko  and 
W.  Lemley,  J  Infrared  and  MM-waves  1  (1980),  85. 

“Guiding  and  Control  of  Millimeter  Waves  by  Surface  Plasmon  Phenomena  in  Semi-conductors,"  A.  V.  Nurmikko, 
D  M  Bolle.  and  S.  Talisa,  J.  Infrared  and  MM-waves  1  (1980),  3. 

"Nonlinear  Transmission  of  Picosecond  10.6  urn  Pulses  in  InSb,”  B.  D.  Schwartz  and  A.  V.  Nurmikko,  Optics  Let¬ 
ters  5  (1980),  5. 

"Picosecond  Infrared  Spectroscopy  in  Narrow-Gap  Semiconductors,"  B.  D.  Schwartz  and  A.  V.  Nurmikko,  in 
Picosecond  Phenomena  II,  Springer-Verlag  Series  in  Chemical  Physics  14  (1980),  p.  303. 

Supported  by  the  Materials  Research  Laboratory / NSF,  the  National  Science  Foundation,  the  Air  Force  Office  of 
Scientific  Research,  and  the  U  S.  Department  of  Energy/Solar  Energy  Research  Institute. 


Theory  of  the  Hall  Effect  In  Amorphous  Semiconductors 

Charge  transport  in  amorphous  semiconductors  is  generally  assumed  to  occur  in  two  regimes:  extended  state 
motion  for  carriers  with  energies  above  the  mobility  edge,  and  transport  among  Anderson  localized  states  for  car¬ 
riers  with  energies  below  the  edge.  In  the  former  regime,  the  scattering  is  strong  and  the  mean  free  path  of  the 
order  of  the  interatomic  spacing,  too  small  for  conventioanl  Bloch-Wilson  band  theory  to  apply.  The  random 
phase  model  has  been  used  here,  where  phase  coherence  of  the  one  electron  wave  function  is  taken  to  be  ran¬ 
dom  from  one  site  to  the  next.  In  one  part  of  our  study,  it  was  shown  that  this  model  leads  to  difficulties  which  can 
be  avoided  by  taking  the  phases  to  be  correlated  over  the  mean  free  path  (phase  correlation  distance)  but  taking 
phase  differences  of  distinct  pairs  of  sites  to  be  uncorrelated.  In  particular,  it  was  shown  that  the  previously 
derived  expression  for  the  Hall  mobility  uH  is  recovered — an  expression  which  has  been  successfully  applied  to 
transport  data  in  a  number  of  disordered  systems  and  liquid  semiconductors. 

In  fitting  the  temperature  dependence  of  a  number  of  transport  properties  of  amorphous  semiconductors,  many 
workers  (notably  N.  F.  Mott  and  his  school)  assume  a  two  band  model  in  which  the  partial  Hall  effect  due  to  the 
localized  states  is  negligible  in  comparison  to  that  of  the  extended  states  described  above.  In  order  to  investigate 
the  validity  of  this  assumption,  we  have  studied  the  Hall  effect  due  to  hopping  transport  in  systems  with  both  ran¬ 
dom  site  position  and  random  local  site  energies  (earlier  studies  incorporated  the  former  but  neglected  the  latter). 
We  have  carried  out  a  percolation-theoretic  calculation  of  the  Hall  effect  for  this  case  and  find  that  for  statistical 
reasons  the  Hall  mobility  is  indeed  very  small  as  had  been  assumed. 


Principal  Investigator:  L.  R.  Friedman  (Physics). 
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Publications:  "Hall  Effect  in  the  Random  Phase  Model  II,"  L.  Friedman,  Phil.  Mag.  41  (1980),  347. 

“The  Hall  Effect  in  the  Variable  Range  Hopping  Regime.’'  L.  Friedman  and  M.  Poliak,  to  be  submitted  to  Phil.  Mag. 
Supported  by  the  Materials  Research  Laboratory/ NSF. 


X-Ray  Topography 

Studies  of  oxidation  in  zinc  by  x-ray  topography  and  correlation  with  surface  features  observed  by  scanning  elec¬ 
tron  microscopy,  microprobe  and  optical  and  Nomarski  microscopy  have  continued.  Electron  microprobe  meas¬ 
urements  of  the  oxygen  and  zinc  x-ray  intensities,  in  regions  of  crystal  surface  selected  on  the  basis  of  optical  and 
x-ray  topographic  observations,  indicate  that  a  correlation  exists  between  the  surface  defects,  oxide  growth  and 
lattice  defects.  Berg  Barrett  x-ray  topographs  show  the  dislocation  structure  within  the  crystal  and  its  alteration 
during  oxide  growth  on  the  surface.  Loop  and  spiral  dislocations  are  observed  in  the  presence  of  large  numbers  of 
slip  dislocations  of  all  three  basal  Burger's  vectors.  Dislocation  spirals  are  often  composed  of  bundles  of  these 
dislocations.  Dislocation  loop  growth  and  climb  of  dislocation  spirals,  associated  with  vacancy  injection  into  the 
lattice  during  oxidation,  are  observed  during  oxidation  at  room  temperature.  Hexagonal  crystallographic  pits  are 
observed  optically  on  the  crystal  surface  which  give  a  correlatable  contrast  in  Berg  Barrett  topographs.  The 
results  indicate  that  the  oxidation  process  is  closely  related  to  and  perhaps  controlled  by  the  lattice  defect  density 
and  dislocation  arrangements  near  the  surface. 

In  addition,  a  small  effort  has  been  devoted  to  activities  based  on  work  carried  out  in  the  past  on  x-ray  topographic 
observations  of  ferromagnetic  domains  and  associated  contrast  effects  and  on  small  angle  grain  boundaries  in 
zinc  crystals  These  activities  were  in  connection  with  the  U.S.-France  Seminar  on  Applications  of  x-ray  topogra¬ 
phy  to  Materials  Science  held  March  12-14,  1980  in  Paris,  France.  As  an  outgrowth  of  these  activities  we  have 
also  begun  preparations  for  some  work  using  synchrotron  radiation  for  topography  in  cooperation  with  Dr.  Masao 
Kuriyama  of  the  National  Bureau  of  Standards  in  Washington,  DC. 

Principal  Investigator  B  Roessler  (Engineering). 

Personnel:  K  C  Yoo  (Engineering)  and  W.  Oates  (Engineering). 

Supported  by  the  Matenals  Research  Laboratory/ NSF. 


Nonlinear  Wave  Interactions  in  Condensed  Matter 

The  study  of  nonlinear  (anharmonic)  properties  of  solids  and  liquids  has  a  long  history.  Many  questions  in  this 
area,  however,  remain  unanswered.  One  of  these  questions  concerns  parametric  amplification  under  special 
conditions.  We  conducted,  therefore,  a  study  of  this  problem. 

Nonlinear  interactions  of  two  ultrasonic  waves,  propagating  collinearty  in  a  dispersionless  liquid  medium,  are 
investigated  experimentally  for  both  degenerate  and  nondegenerate  cases.  Results  are  described  of  the  first 
experimental  observation  of  parametric  amplification  of  the  signal  wave  for  the  degenerate  case.  For  the  nonde¬ 
generate  case,  no  such  amplification  of  the  amplitude  of  the  signal  wave  is  observed.  These  experimental  results 
are  found  to  be  in  good  agreement  with  the  results  of  the  theoretical  analysis  presented  in  the  paper. 

Principal  Investigator:  C.  Elbaum  (Physics  and  Applied  Mathematics). 

Personnel:  A.  Hikata  (Applied  Mathematics)  and  H.  Kwun  (Materials  Research  Laboratory). 

Publication:  ‘Nonlinear  Interactions  of  Two  Ultrasonic  Waves  In  a  Dispersionless  Liquid  Medium,"  H.  Kwun, 
A.  Hikata,  and  C.  Elbaum,  J.  Appl.  Phys.  51  (1980),  4645, 

Supported  by  the  Materials  Research  Laboratory/NSF. 


mmmmrnmmmmmi,, , 


Lecturers  in  the  Materials  Research  Laboratory  Program,  1979-80 

September  20,  1979 

Solid  State  Seminar:  Physics  Department,  Dr.  Erio  Tosatti,  international  Center  for  Theoretical  Physics 
and  the  University  of  Trieste.  Italy.  INSTABILITIES,  PHASE  TRANSITIONS  AND  MEETING  OF 
CRYSTAL  SURFACES 

September  24,  1 979 

Physics  Colloquium:  Professor  F.  M.  Pipkin,  Department  of  Physics,  Harvard  University, 
MEASUREMENTS  OF  THE  FINE  STRUCTURE  OF  HYDROGENIC  ATOMS. 

September  26,  1979 

Materials  Science  Seminar:  Division  of  Engineering,  Professor  John  Hirth,  Department  of  Metallurgical 
Engineering,  Ohio  State  University,  MICROSTRUCTURAL  MODELS  FOR  FRICTION  AND  WEAR. 

September  27,  1979 

Joint  Matenals  Science  and  Fracture  Group  Seminar:  Division  of  Engineering,  Dr.  John  Hancock  of  the 
Department  of  Mechanical  Engineering,  University  of  Glasgow,  DUCTILE  FRACTURE. 

Solid  State  Seminar:  Physics  Department,  Dr.  T.  N.  Taylor,  Los  Alamos  Scientific  Lab,  SURFACE 
STUDIES  OF  U02  BY  LOW  ENERGY  ELECTRON  DIFFRACTION  AND  ION  SCATTERING 
SPECTROSCOPY. 

October  3,  1979 

Materials  Science  Seminar:  Division  of  Engineering,  Dr.  A.  F.  Szewczyk,  Brown  University, 
MICROSTRUCTURAL  OBSERVATIONS  OF  MONOTONICALLY  AND  CYCLICALLY  DEFORMED 
DUAL-PHASE  STEEL. 

October  4,  1979 

Theoretical  Seminar:  Physics  Department,  Dr.  Goran  Senjanovic,  Department  of  Physics,  University  of 
Maryland,  BROKEN  SYMMETRIES  AT  HIGH  TEMPERATURE. 

October  11,  1979 

Solid  State  Seminar:  Physics  Department,  Professor  Brooke  Gregory,  Trinity  College. 
MAGNETIZATION  OF  MULTIPLY-CONNECTED  SUPERCONDUCTORS. 

October  15,  1979 

Solid  Mechanics  Seminar:  Division  of  Engineering,  Professor  K.  B.  Broberg,  Lund  Institute  of  Technol¬ 
ogy,  CRACK  PROPAGATION  BEYOND  THE  RAYLEIGH  WAVE  VELOCITY. 

October  18,  1979 

Theoretical  Seminar:  Physics  Department,  Dr.  Tom  Clark,  Los  Alamos,  AXIAL  CURRENTS, 
SUPERCURRENTS  AND  ANOMALIES. 

Solid  State  Seminar:  Physics  Department,  Dr.  Walter  Bloss,  Physics  Department,  Brown  University, 
EXCHANGE-CORRELATION  INDUCED  TRANSITION  IN  SILICON  INVERSION  LAYERS. 
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October  22,  1979 

Solid  Mechanics  Seminar:  Division  of  Engineering,  Dr.  Viggo  Tvergaard,  Department  of  Solid  Mechan¬ 
ics,  Technical  University  of  Denmark,  INFLUENCE  OF  VOIDS  ON  SHEAR  BAND  INSTABILITIES 
UNDER  PLANE  STRAIN  CONDITIONS. 

October  23,  1979 

Materials  and  Electrical  Sciences  Seminar:  Division  of  Engineering,  Professor  Julian  Aulytner,  Central 
Laboratory  of  X-Ray  and  Electron  Microscopy,  Polish  Academy  of  Sciences,  Warsaw,  Poland,  X-RAY 
AND  ELECTRON  MICROSCOPIC  INVESTIGATION  OF  ION-IMPLANTED  CRYSTALS. 

October  25,  1979 

Solid  State  Seminar:  Physics  Department,  Professor  L.  Friedman,  St.  Andrews  University,  Scotland, 
and  Brown  University,  THE  HALL  EFFECT  IN  AMORPHOUS  AND  LIQUID  SEMICONDUCTORS. 

October  26,  1979 

Chemistry  Colloquium:  Dr.  F.  Stillinger,  Bell  Laboratory,  BRIDGING  STRUCTURAL  CHEMISTRY  AND 
STATISTICAL  MECHANICS  WITH  THE  POLARIZATION  MODEL. 

November  1,  1979 

Solid  State  Seminar:  Physics  Department,  Dr.  U.  Strom,  Naval  Research  Laboratory,  DISORDER 
PHENOMENA  IN  SUPERIONIC  CONDUCTORS. 

Solid  Mechanics  Seminar:  Division  of  Applied  Mathematics,  Dr.  Paolo  Podio-Guidugli,  Department  of 
Mathematics,  Carnegie-Mellon  University,  SOME  QUALITATIVE  ASPECTS  IN  STRUCTURAL 
STABILITY. 

November  5,  1979 

Physics  Colloquium:  Dr.  L.  Nosanow,  National  Science  Foundation,  SPIN  ALIGNED  HYDROGEN;  A 
POSSIBLE  NEW  QUANTUM  SYSTEM. 

November  7,  1979 

Joint  Applied  Mathematics  and  Solid  Mechanics  Seminar:  Divisions  of  Applied  Mathematics  and  Engi¬ 
neering,  Professor  William  Prager,  Savognin,  Switzerland,  SOME  UNEXPECTED  RESULTS  IN 
STRUCTURAL  OPTIMIZATION. 

November  15,  1979 

Solid  State  Seminar:  Physics  Department,  Professor  E.  Abrahams,  Rutgers  University, 
LOCALIZATION  AND  NON-OHMIC  EFFECTS  IN  TWO  DIMENSIONS. 

November  16. 1979 

Chemistry  Colloquium:  Dr.  M.  Saunders,  Yale  University,  NON-CLASSICAL  DELOCALIZATION  IN 
CATIONS  BY  NMR. 

November  19, 1979 

Solid  Mechanics  Seminar:  Division  of  Engineering,  Dr.  C.  T.  Chon,  Ford  Scientific  Research  Labs, 
ANALYSIS  OF  COMPOSITE  DRIVE  SHAFT  AND  TUBULAR  LAP  JOINT. 
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November  20,  1979 

Solid  State  Seminar:  Physics  Department,  Dr.  James  Smith,  Los  Alamos  Scientific  Lab, 
SUPERCONDUCTIVITY  IN  F-ELECTRON  SYSTEMS. 

Chemistry  Colloquium:  Dr.  R.  Silbey,  Massachusetts  Institute  of  Technology,  ELECTRONIC  ENERGY 
TRANSFER  IN  DISORDERED  SOLIDS. 

November  26,  1979 

Physics  Colloquium:  Professor  N.  Ashcroft,  Department  of  Physics,  Cornell  University,  VERY  DENSE 
PHASES  OF  HYDROGEN. 

November  28,  1979 

Materials  Science  Seminar:  Division  of  Engineering,  Dr.  Alan  T.  English,  Supervisor  Reliability  Studies, 
Bell  Laboratories.  Murray  Hill,  New  Jersey,  PHOTOLITHOGRAPHY  UTILIZING  INORGANIC 
RESISTOR  MATERIALS. 

November  29,  1979 

Solid  State  Seminar:  Physics  Department,  Dr.  David  Zehner,  Oak  Ridge  National  Laboratory, 
GEOMETRICS:  SURFACE  STRUCTURE  VIA  HIGH  ENERGY  ION  SCATTERING. 

December  3,  1979 

Physics  Colloquium:  Professor  Nihat  Barker,  Department  of  Physics,  Massachusetts  Institute  of  Tech¬ 
nology,  RENORMALIZATION  GROUP  THEORY:  PHASE  TRANSITIONS  IN  ABSORBED  SYSTEMS. 

December  5.  1979 

Materials  Science  Seminar:  Mr.  Henry  Cialone,  Materials  Science  Division,  Brown  University,  THE 
ROLE  OF  HYDROGEN  IN  THE  DUCTILE  FRACTURE  OF  PLAIN  CARBON  STEELS. 

December  6,  1979 

Solid  State  Seminar:  Physics  Department,  Dr.  Norton  Lang,  IBM,  Yorktown  Heights,  CORE  HOLES 
AND  CHEMISORPTION. 

December  7, 1979 

Solid  Mechanics  Seminar:  Division  of  Engineering,  Professor  Sol  Bodner,  Division  of  Mechanical  Engi¬ 
neering,  Technion,  Haifa,  Israel,  DEVELOPMENT  AND  APPLICATION  OF  MULTI  DIMENSIONAL 
HARDENING  LAWS  FOR  ELASTIC-VISCOPLASTIC  MATERIALS. 

December  10,  1979 

Solid  Mechanics  Seminar:  Division  of  Engineering,  Dr.  John  Dempsey,  Department  of  Civil  Engineer¬ 
ing,  Northwestern  University,  SOME  REMARKS  ON  SINGULAR  SOLUTIONS  IN  PLANE  ELASTICITY. 

December  12,  1979 

Solid  Mechanics  Seminar:  Division  of  Engineering,  Professor  Eric  Belssner,  University  of  California  at 
San  Diego,  NEW  THOUGHTS  ON  THE  PROBLEM  OF  BENDING  AND  TWISTING  OF  BEAMS. 

Oecember  13, 1979 

Solid  State  Seminar  Physics  Department,  Dr.  H.  K.  Heniah,  Department  of  Physics,  Pennsylvania  State 
University,  RELAXATION  SEMICONDUCTORS. 
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December  14,  1979 

Physics  Colloquium:  Professor  T,  W.  Ruljgrok,  Institute  voor  Theoretische  Fysica,  The  Netherlands,  ON 
THE  HAMILTONIAN  AND  LAGRANGIAN  FORMULATION  OF  CLASSICAL  DYNAMICS  FOR 
PARTICLES  WITH  SPIN. 

January  7.  1980 

Solid  Mechanics  Seminar  Division  of  Engineering,  Professor  Michael  P.  Cleary,  Department  of 
Mechanical  Engineering,  Massachusetts  Institute  of  Technology,  SOME  MECHANICS  ASPECTS  OF 
UNDERGROUND  FRACTURING. 

January  9,  1980 

Materials  Science  Seminar:  Division  of  Engineering,  Dr.  M.  P.  Puls,  Whiteshell  Nuclear  Research 
Establishment,  Atomic  Energy  of  Canada  Ltd.,  Pinawa,  Manitoba,  STRAIN  ENERGY  EFFECTS  AND 
HYDRIDE  INDUCED  CRACK  GROWTH. 

January  11,1 980 

Chemistry  Colloquium:  K.  Eisenthal,  Columbia  University,  PICOSECOND  LASER  STUDIES  OF 
PROTON  TRANSFER  PROCESSES. 

Solid  Mechanics  Seminar:  Division  of  Engineering,  Professor  Atle  Gjelsvik,  Department  of  Civil  Engi¬ 
neering,  Columbia  University,  PLASTIC  ANALYSIS  IN  OFFSHORE  STRUCTURES. 

February  1,  1980 

Chemistry  Colloquium:  Or.  T.  Madey,  National  Bureau  of  Standards,  USE  OF  ELECTRON 
STIMULATED  DESORPTION  FOR  DETERMINING  THE  STRUCTURES  OF  ABSORBED 
MOLECULES. 

February  7,  1980 

Solid  State  Seminar:  Physics  Department,  Dr.  David  Emin,  Sandia  Laboratory,  Albuquerque,  New  Mex¬ 
ico,  POLARON  VERSUS  DEFECT  PICTURE  IN  AMORPHOUS  SEMICONDUCTORS. 

February  8,  1980 

Chemistry  Colloquium:  Dr.  Gary  Pruentt,  University  of  Pennsylvania,  STATE  SELECTED  REACTIONS 
UNDER  MOLECULAR  BEAM  CONDITIONS. 

February  14, 1980 

Division  of  Applied  Mathematics  Seminar:  Dr.  Jerry  Nedelman,  Cornell  University  Graduate  School  of 
Medical  Sciences.  Biomathematics  Division,  FACILITATED  DIFFUSION  OF  OXYGEN  AND  CARBON 
MONOXIDE  IN  THE  HIGH  AFFINITY  REGIME. 

Energy  Lecture:  Physics  Department,  Professor  E.  A.  Farber,  University  of  Florida,  Gainesville,  Florida, 
SOLAR  ENERGY:  ITS  CONSERVATION  AND  UTILIZATION. 

February  15, 1980 

Solid  State  Seminar:  Physics  Department,  Dr.  Jorge  Hirsch,  University  of  Chicago  and  University  of  Cal¬ 
ifornia.  San  Diego,  SINGULAR  THERMODYNAMIC  PROPERTIES  OF  RANDOM  MAGNETIC  CHAINS. 


February  20,  1980 

Materials  Science  Seminar:  Division  of  Engineering,  John  Fisher,  Brown  University,  VOID 
NUCLEATION  IN  SPHEROIDIZED  STEELS  DURING  TENSILE  DEFORMATION. 

February  21,  1980 

Solid  State  Seminar:  Physics  Department,  Dr.  James  W.  Davenport,  Brookhaven  National  Laboratory, 
Long  Island,  CHEMISORPTION  AND  SOLUTION  OF  HYDROGEN  IN  TRANSITION  METALS. 

February  22,  1980 

Chemistry  Colloquium:  Dr.  Mark  M.  Green,  New  York  Polytechnical  Institute,  MASS  SPECTROMETRY 
IS  NOT  SPECTROMETRY. 

February  29.  1980 

Solid  State  Seminar:  Physics  Department,  Professor  E.  Burstein,  University  of  Pennsylvania,  GIANT 
RAMAN  SCATTERING  BY  MOLECULES  ON  METAL  SURFACES. 

March  6,  1980 

Solid  State  Seminar:  Physics  Department,  Dr.  D.  R.  Hamann,  Bell  Telephone  Laboratories, 
CHEMISORPTIVE  BONDS  ON  SEMICONDUCTORS:  STRUCTURE  AND  SPECTRA. 

March  7.  1980 

Chemistry  Colloquium:  Dr.  J.  R.  Dorfman,  University  of  Maryland,  ELEMENTARY  EXCITATIONS  IN 
CLASSICAL  LIQUIDS. 

Division  of  Applied  Mathematics  Colloquium:  Professor  N.  J.  Zabusky,  University  of  Pittsburgh, 
VORTEX  DYNAMICS:  STABILITY,  INTERACTIONS  AND  BREAKING  OF  COHERENT 
STRUCTURES. 

March  10,  1980 

Physics  Colloquium:  Professor  James  S.  Langer,  Department  of  Physics,  Camegie-Mellon  University, 
INSTABILITIES  AND  PATTERN  FORMATION  IN  SOLIDIFICATION. 

March  12,  1980 

Solid  Mechanics  Seminar:  Division  of  Engineering,  Professor  Alan  Needleman,  Brown  University, 
FLOW-LOCALIZATION  IN  THE  PLANE  TENSILE  TEST. 

March  13,  1980 

Solid  State  Seminar:  Physics  Department,  Professor  Leo  P.  Kadanoff,  University  of  Chicago,  THE 
WORK  OF  SATO,  MIWA  AND  JIMBO:  THE  ISING  MODEL  AND  THE  DIRAC  EQUATION. 

March  14.  1980 

Division  of  Applied  Mathematics  Seminar:  Dr.  Nicholas  Alikakos,  Purdue  University,  ON  THE 
EQUATION  OF  POROUS  MEDIA. 

March  17. 1980 

Physics  Colloquium:  Or.  Arthur  Aahkln,  Bell  Telephone  Laboratories.  TRAPPING  OF  PARTICLES  BY 
LASER  LIGHT. 
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March  20,  I960 

Solid  Mechanics  Seminar:  Division  ol  Engineering,  Dr.  R.  D.  James,  Department  of  Aeronautical  Engi¬ 
neering  and  Mechanics,  University  of  Minnesota,  CO-EXISTENT  PHASES  IN  ELASTICITY. 

Solid  State  Seminar:  Physics  Department,  Professor  Samuel  C.  Fain,  Jr.,  University  of  Washington, 
Seattle,  PHASE  TRANSITIONS  IN  PHYSISORBID  LAYERS. 

March  24.  1980 

Solid  Mechanics  Seminar:  Division  of  Engineering,  Professor  A.  J.  Holzer,  Department  ol  Mechanical 
Engineering,  Camegie-Mellon  University,  DYNAMIC  MECHANICAL  PROPERTIES  USING  SHORT 
LOAD  CELLS. 

March  27,  1980 

Energy  Lecture:  Department  of  Physics,  Dr.  Richard  T.  Kennedy,  Commissioner,  United  States  Nuclear 
Regulatory  Commission,  AVAILABILITY  OF  PRIMARY  ENERGY  RESOURCES  AND  TECHNOLOGY. 

April  11,  1980 

Chemistry  Colloquium:  Dr.  Hans  Nowotny,  University  of  Connecticut,  STRUCTURAL  PRINCIPLES  OF 
ALLOY  PHASES. 

April  14,  1980 

Expository  Talk:  Division  of  Applied  Mathematics,  Mr.  Paul  Massatt,  Brown  University,  STABILITY  AND 
FIXED  POINTS  OF  DISSIPATIVE  SYSTEMS. 

April  17,  1980 

Solid  State  Seminar:  Department  of  Physics,  Dr.  Douglas  Osheroff,  Bell  Telephone  Laboratories, 
NUCLEAR  ANTIFERROMAGNETIC  RESONANCE  AND  SOLID  He3. 

April  21,  1980 

Solid  Mechanics  Seminar:  Division  of  Engineering,  Dr.  Shiro  Kubo,  Osaka  University,  Japan,  A 
PHENOMENOLOGICAL  ANALYSIS  OF  CREEP  CRACK  PROPAGATION  IN  DUCTILE  MATERIALS. 

Physics  Colloquium:  Professor  Paul  Meijer,  Catholic  University  of  America.  Washington,  DC. 
SUPERCOOLING  OF  WATER  AND  OTHER  LIQUIDS. 

Materials  Science  Seminar:  Division  of  Engineering,  Dr.  -Ing  J.  Heinrich,  Deutsche  Forschungs-und 
Versuchsanstalt  fur  Luft-und  Raumfahrt,  E.V.,  PROCESSING  OF  STRUCTURAL  CERAMICS. 

April  28.  1980 

Solid  Mechanics  Seminar:  Or.  J.  W.  Nunziato,  Sandia  Laboratories,  A  THEORY  OF  MULTI-PHASE 
MIXTURES  WITH  APPLICATIONS  TO  TWO-PHASE  FLOWS. 

May  1,  I960 

Solid  State  Seminar:  Department  of  Physics,  Dr.  D.  Greywall.  Bell  Telephone  Laboratories,  THE  He3 
QUASIPARTICLE  SPECTRUM  FOR  DILUTE  SOLUTIONS  OF  He*  IN  SUPERFLUID  He4. 

May  7.  1960 


Materials  Science  Seminar:  Professor  R.  O.  Heckroodt,  Department  of  Materials  Science,  University  of 
Capetown,  South  Africa,  TELLURITE  GLASSES  AND  RARE  EARTH  PHOTOCHROMISM. 
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May  8.  1980 

Solid  State  Seminar:  Department  of  Physics,  Dr.  R.  Penn,  National  Bureau  of  Standards,  EFFECT  OF 
d-BAND  HOLES  ON  THE  PHOTOEMISSION  SPECTRUM  OF  Ni. 

Solid  Mechanics  Seminar:  Division  of  Engineering,  Professor  William  N.  Findley,  Brown  University, 
THE  COMBINED  TENSION  AND  TORSION  CREEP  AND  RELAXATION  OF  2618  A  AT  200  C  UNDER 
VARIABLE  NONLINEAR  STRESSES 

May  9,  1980 

Materials  Science  Seminar:  Division  of  Engineering,  Professor  I.  M.  Bernstein,  Department  of  Metal¬ 
lurgy  and  Materials  Science,  Carnegie-Mellon  University,  MICROMECHANICS  OF  CLEAVAGE 
FRACTURE  IN  EUTECTOID  STEELS. 

May  12.  1980 

Solid  Mechanics  Seminar:  Division  of  Engineering,  Professor  Janusz  Klepaczko,  Institute  of  Funda¬ 
mental  &  Technological  Research,  Warsaw,  Poland  and  the  Division  of  Engineering,  Brown  University, 
LOADING  RATE  SPECTRA  IN  FRACTURE  DYNAMICS. 

May  13,  1980 

Seminar  on  Identification  and  Control:  Division  of  Applied  Mathematics,  Professor  B.  Jakubczyk,  Polish 
Academy  of  Science,  Warsaw,  Poland  and  the  Division  of  Applied  Mathematics,  Brown  University, 
REALIZATIONS  OF  NONLINEAR  SYSTEMS. 

May  19,  1980 

Solid  Mechanics  Seminar:  Division  of  Engineering,  Professor  William  A.  Nash,  Department  of  Civil 
Engineering,  University  of  Massachusetts,  SEISMIC  EFFECTS  ON  LIQUID  STORAGE  TANKS. 

June  20,  1980 

Materials  Science  Seminar:  Division  of  Engineering,  Dr.  Jean  Louis  Chermant,  Laboratory  of  Crystal¬ 
lography  and  Chemistry  of  Solids,  University  of  Caen,  Caen,  France,  SINTERING:  THEORY  AND 
TOPOLOGY. 
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List  of  Staff  in  tha  Mata  rial*  Raaaarch  Laboratory  Program 

APPLIED  MATHEMATICS 

NAME 

TITLE 

C.  Elbaum 

Professor 

H.  Kolsky 

(Physics) 

Professor 

(Engineering) 

A.  C.  Pipkin 

(Physics) 

Professor 

A.  Hikata 

Professor 

(Research) 

L.  Mannion 

Research  Assistant 

J.  Mosquera 

Research  Assistant 

J.  Marsella 

Technical  Assistant 

P.  Russo 

Technical  Assistant 
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List  of  Staff  in  the  Materials  Research  Laboratory  Program 

CHEMISTRY 


NAME 

TITLE 

R.  H.  Cole 

P.  J.  Estrup 

Professor 

Professor 

E.  F.  Greene 

W.  M.  Risen,  Jr. 

A.  Wold 

(Physics) 

Professor 

Professor 

Professor 

(Engineering) 

K.  Dwight 

Associate  Professor 

G.  J.  Diebold 

Assistant  Professor 

J.  Hormadly 

S.  N.  Subbarao 

Research  Associate 

Research  Associate 

R.  A.  Barker 

D.  M.  Barnes 

J.  Hayden 

E.  1.  Kamitsos 

J.  Keeley 

B.  Khazai 

V.  S.  Nguyen 

J.  D.  Passaretti 

M.  Pickering 

C.  H.  Tzinis 

C.  H.  Wfndisch 

Research  Assistant 

Research  Assistant 

Research  Assistant 

Research  Assistant 

Research  Assistant 

Research  Assistant 

Research  Assistant 

Research  Assistant 

Research  Assistant 

Research  Assistant 

Research  Assistant 

J.  Koenitzer 

Research  Undergraduate  Student 
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List  of  Staff  in  the  Materials  Research  Laboratory  Program 


ENGINEERING 


NAME 


TITLE 


R.  J.  Asaro 
D.  H.  Avery 
R.  J.  Clifton 
J.  Duffy 

L.  B.  Freund 

M.  Glicksman 
J.  Gurland 

G.  S.  Heller 

H.  Kolsky 


J.  J.  Loferski 
J.  R.  Rice 
M.  H.  Rich  man 
B.  Roessler 
J.  Tauc 

J.  H.  Weiner 

A.  Wold 


A.  Needleman 
A.  V.  Nurmikko 


Professor 

Professor 

Professor 

Professor 

Professor 

Professor 

Professor 

Professor 

Professor 

(Applied  Mathematics) 
(Physics) 

Professor 

Professor 

Professor 

Professor 

Professor 

(Physics) 

Professor 

(Physics) 

Professor 

(Chemistry) 

Associate  Professor 
Associate  Professor 


K.  K.  Lo 


Assistant  Professor 


J.  Klepaczko  Visiting  Professor,  Research 

J.  Shewchun  Visiting  Professor 


V.  Tvergaard 


Visiting  Associate  Professor 


P.  Kumar 
X.  Marken scoff 


Visiting  Assistant  Professor 
Visiting  Assistant  Professor 


J.  Piekoszewski 


Senior  Research  Associate 


S.  Kubo 


Postdoctoral  Research  Fellow 


L.  Bieman 
R.  Dmowska 
Z.  Vardeny 


Research  Associate 

Research  Associate,  Geological  Sciences 
Research  Associate 


Abou-Sayed 

Ajo 

M.  Baik 

Burgers 

Burk 

Case 

Cfaione 


Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
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R.  Deshmukh 
A.  Douglas 
W.  J.  Drugan 
Y.  E  rcil 
J  R.  Fisher 
A.  Gilat 
O  Gregory 
J.  Hong 
M.  Kachanov 
F.  K.  Lehner 
W.  Lemley 

C.  H.  Li 
V.  C.  Li 

M.  Magida 
A.  Morrone 

J.  Mosquera 
J  Pan 

D.  Perchak 
D  Pierce 
A  Rosakis 

A.  Rubenstein 
A.  Ruina 
T.  L.  Sham 
A.  Szewczyk 

N.  Triantafyllidis 

K.  C.  Yoo 
R.  Zeller 


Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 


J.  C.  Lambropolis 

T.  Simon 

Undergraduate  Student 
Undergraduate  Student 

E.  E.  Crisman 

R.  H.  Hawley 

L.  Hermann 

Senior  Research  Engineer 
Senior  Research  Engineer 
Senior  Research  Engineer 

B  Agbara 

J  Beall 

R.  Kershaw 

T.  Kirst 

Research  Engineer 
Research  Engineer 
Research  Engineer 
Research  Engineer 

R.  Beaulieu 

J.  W.  Fogarty 

G.  J.  LaBonte 

W  Oates 

P  Rush 

H.  Stanton 

Technical  Assistant 
Technical  Assistant 
Technical  Assistant 
Technical  Assistant 

Tev.:  teal  Assistant 
Technical  Assistant 
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List  of  Staff  In  the  Materials  Research  Laboratory  Program 
MATERIALS  RESEARCH  LABORATORY 


NAME  TITLE 

L  Simpson  Senior  Research  Associate 


R.  C.  Abeyaratne 
H.  Araki 
E.  Berkcan 
R.  D.  Kenway 
H.  Kwun 


Research  Associate 
Research  Associate 
Research  Associate 
Research  Associate 
Research  Associate 
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List  of  Staff  in  ths  Materials  Research  Laboratory  Program 

PHYSICS 


NAME 

TITLE 

P.  J.  Bray 

Professor 

L.  N  Cooper 

Professor 

C.  Elbaum 

Professor 

(Applied  Mathematics) 

P  J  Estrup 

Professor 

(Chemistry) 

A.  Houghton 

Professor 

H.  Kolsky 

Professor 

(Applied  Mathematics) 
(Engineering) 

H.  J.  Maris 

Professor 

J.  J.  Quinn 

Professor 

G.  M.  Seidel 

Professor 

P.  J.  Stiles 

Professor 

J.  Tauc 

Professor 

(Engineering) 

S.  C.  Ying 

Associate  Professor 

L.  Friedman 

Visiting  Professor 

K.  H.  Lau 

Research  Associate 

K.  P.  Bohnen 

Visiting  Research  Associate 

G.  Blanchet 

Research  Assistant 

F.  Bucholtz 

Research  Assistant 

M.  L.  Buess 

Research  Assistant 

J.  Cameron 

Research  Assistant 

M.  C.  Chang 

Research  Assistant 

R.  W.  Cline 

Research  Assistant 

W.  J.  Dell 

Research  Assistant 

A.  E.  Geissburger 

Research  Assistant 

S.  G.  Greenbaum 

Research  Assistant 

D.  Grempel 

Research  Assistant 

M.  D.  Han 

Research  Assistant 

1.  A.  Harris 

Research  Assistant 

C.  Huber 

Research  Assistant 

T.  E.  Huber 

Research  Assistant 

J.  S.  Kotol 

Research  Assistant 

J.  Lee 

Research  Assistant 

M.  Lega 

Research  Assistant 

A.  LeGrand 

Research  Assistant 

M.  Lui 

Research  Assistant 

R.  Mayer 

Research  Assistant 

P.  O’Connor 

Research  Assistant 

G.  Park 

Research  Assistant 

B.  Parker 

Research  Assistant 

S.  Ray 

Research  Assistant 

G.  Robertshaw 

Research  Assistant 
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S.  Semancik 
M.  Spitzer 
E.  Vera 
K  C  Woo 
K.  S.  Yi 


Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
Research  Assistant 
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Index,  Principal  Investigators 


A  Abeyaratne,  R.  C.,  14 
Asaro,  R.  J.,  14,  17,  23 
Avery,  D.  H.,  23 

B  Bray,  P.  J.,  29-31 

C  Clifton,  R.  J.,  2,  4-5 
Cole,  R.  H.,  31 
Cooper,  L.  N,  39 

D  Diebold.  G.  J„  37 
Duffy.  J.,  3,  10 

E  Elbaum,  C.,  2,  25,  43-45,  52 
Estrup,  P.  J.,  34 

P  Freund,  L,  B„  5,  8-9 
Friedman,  L.  R.,  51 

G  Greene.  E.  F.,  36 
Gurland,  J.,  12-13 

H  Houghton,  A.,  26 

K  Kenway,  R  D „  40 
Kolsky,  H  ,  47 

L  Lo,  K.  K.,  10 
Loferski,  J.,  50 

M  Maris,  H.  J.,  42 

N  Needleman,  A.,  14-17,  22 

Nurmikko,  A.  V.,  51 

P  Pipkin,  A.  C.,  1 1 

Q  Quinn,  J.  J.,  39 

R  Rice,  J.  R.,  15,  19,  21-22 
Richman,  M.  H.,  18 
Risen,  Jr.,  W.  M.,  28-29,  37,  41 
Roessler,  B.,  50,  52 

S  Seidel,  G.  M„  43 

Shewchun,  J.,  50 
Simpson,  L.  A.,  11-12 
Stiles,  P.  J.,  35,  41 

T  Tauc,  J.,  25,  27 

W  Weiner,  J.  H„  4,  46 

Wold,  A.,  48 


Y  Ying,  S.  C„  35-38 


END 
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